60 IF ( a.GT.0.0 .AND. r.GT.0.0 ) go
to 10
61 WRITE (*,*)
'In GENGAM - Either (1) Location param A <= 0.0 or'
62 WRITE (*,*)
'(2) Shape param R <= 0.0 - ABORT!'
63 WRITE (*,*)
'A value: ',a,
'R value: ',r
65 + (
'Location or shape param out of range in GENGAM - ABORT!')
real function gengam(a, r)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to