5 FUNCTION GAUSS3(F,A,B,EPS)
9 DATA W/0.1012285,.2223810,.3137067,.3623838,.0271525,
10 & .0622535,0.0951585,.1246290,.1495960,.1691565,
12 DATA X/0.9602899,.7966665,.5255324,.1834346,.9894009,
13 & .9445750,0.8656312,.7554044,.6178762,.4580168,
20 IF(ABS(Y).LE.DELTA) RETURN
28 1 S8=S8+W(I)*(F(C1+U)+F(C1-U))
31 3 S16=S16+W(I)*(F(C1+U)+F(C1-U))
34 IF(ABS(S16-S8).GT.EPS*(1.+ABS(S16))) GOTO 4
39 IF(ABS(Y).GT.DELTA) GOTO 2
43 7 FORMAT(1X,'GAUSS3....TOO HIGH ACURACY REQUIRED')