Listing of file='29A' on disk='vmedia/glbr_22.wvd.zip'
# Sector 644, program filename = '29A' 0005 DIM T(8) 0007 PRINT "TO END PROGRAM INPUT 0." 0010 PRINT "INPUT 'INTEGRATION LIMIT X'" 0020 INPUT X2 0030 IF X2<0 THEN 10 0040 IF X2=0 THEN 999 0060 X1=0 0080 FOR N1=1 TO 7 0090 L=X2-X1 0100 X=X1 0110 GOSUB 900 0120 Y1=Y 0130 X=X2 0140 GOSUB 900 0150 T(1)=(Y1+Y)/2 0160 N=1 0170 FOR H=1 TO N1 0180 U=0 0190 M=L/(2*N) 0200 FOR J=1 TO 2*N-1 STEP 2 0210 X=X1+J*M 0220 GOSUB 900 0230 U=U+Y 0240 NEXT J 0250 T(H+1)=(U/N+T(H))/2 0260 F=1 0270 FOR J=H TO 1 STEP -1 0280 F=4*F 0290 T(J)=T(J+1)+(T(J+1)-T(J))/(F-1) 0300 NEXT J 0310 N=2*N 0320 NEXT H 0330 I2=T(1)*L 0335 IF N1=1 THEN 350 0340 IF ABS(I1-I2)<=1E-4*ABS(I2) THEN 700 0350 I1=I2 0360 NEXT N1 0370 PRINT "UNABLE TO COMPUTE THE INTEGRAL TO 4 SIGNIFICANT DIGITS." 0380 PRINT "CLOSEST APPROXIMATION IS:" 0700 PRINT 0710 PRINT "INTEGRAL=";1.1283792*I2 0720 PRINT 0730 GOTO 10 0900 Y=EXP(-(X^2)) 0910 RETURN 0999 END