! 正規分布による2項分布の近似 ! XがB(n,p)に従うとき累積確率分布P(X≦k)を計算する OPTION ARITHMETIC NATIVE FOR n=1000 TO 40000 STEP 1000 PRINT n, LET p=1/6 LET m=n*p LET s2=n*p*(1-p) LET t=0 LET k=0 DO LET t=t+1/SQR(2*PI*s2)*EXP(-(k-m)^2/(2*s2)) IF t>0.05 THEN EXIT DO LET k=k+1 LOOP PRINT k-1, LET p=0.16 LET m=n*p LET s2=n*p*(1-p) LET t=0 LET k=n DO LET t=t+1/SQR(2*PI*s2)*EXP(-(k-m)^2/(2*s2)) IF t>0.05 THEN EXIT DO LET k=k-1 LOOP PRINT k+1 NEXT n END