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,
19 IF(abs(
y).LE.delta)
RETURN
30 3 s16=s16+
w(i)*(
f(
c1+u)+
f(
c1-u))
33 IF(abs(s16-s8).GT.
eps*(1.+abs(s16))) goto 4
38 IF(abs(
y).GT.delta) goto 2
42 7
FORMAT(1
x,
'GAUSS4....TOO HIGH ACURACY REQUIRED')