5 C*********GAUSSIAN ONE-DIMENSIONAL INTEGRATION PROGRAM*************
7 FUNCTION GAUSS1(F,A,B,EPS)
11 DATA W/0.1012285,.2223810,.3137067,.3623838,.0271525,
12 & .0622535,0.0951585,.1246290,.1495960,.1691565,
14 DATA X/0.9602899,.7966665,.5255324,.1834346,.9894009,
15 & .9445750,0.8656312,.7554044,.6178762,.4580168,
22 IF(ABS(Y).LE.DELTA) RETURN
30 1 S8=S8+W(I)*(F(C1+U)+F(C1-U))
33 3 S16=S16+W(I)*(F(C1+U)+F(C1-U))
36 IF(ABS(S16-S8).GT.EPS*(1.+ABS(S16))) GOTO 4
41 IF(ABS(Y).GT.DELTA) GOTO 2
45 7 FORMAT(1X,'GAUSS1....TOO HIGH ACURACY REQUIRED')