6 FUNCTION GAUSS4(F,A,B,EPS)
10 DATA W/0.1012285,.2223810,.3137067,.3623838,.0271525,
11 & .0622535,0.0951585,.1246290,.1495960,.1691565,
13 DATA X/0.9602899,.7966665,.5255324,.1834346,.9894009,
14 & .9445750,0.8656312,.7554044,.6178762,.4580168,
21 IF(ABS(Y).LE.DELTA) RETURN
29 1 S8=S8+W(I)*(F(C1+U)+F(C1-U))
32 3 S16=S16+W(I)*(F(C1+U)+F(C1-U))
35 IF(ABS(S16-S8).GT.EPS*(1.+ABS(S16))) GOTO 4
40 IF(ABS(Y).GT.DELTA) GOTO 2
44 7 FORMAT(1X,'GAUSS4....TOO HIGH ACURACY REQUIRED')