Mark Wickens

10 GOSUB400:IFx=<0 THENp=1:GOTO80
20 IFn>100 THENz=x:x=((x/n)^(1/3)-(1-2/9/n))/SQR(2/9/n):GOSUB"LIB0:S6210":GOSUB400:p=y:x=z:GOTO80
30 IFx>400 THENp=0:GOTO80
40 a=EXP(-x/2):p=a:y=2
50 IFn MOD 2<>0 THENz=x:x=SQRx:GOSUB"LIB0:S6210":GOSUB400:p=2*y:a=SQR(2/PI)*a/x:x=z:y=1
60 FORi=y TOn-2 STEP2:a=a*x/i:p=p+a:NEXT
80 y=p:RETURN
400 ONERRORGOTO500:RETURN
500 IFERR=1 THENCLS:ONERRORGOTO0
510 sr=1:RESUME80