* G01EFF Example Program Text * Mark 14 Release. NAG Copyright 1989. * .. Parameters .. INTEGER NIN, NOUT PARAMETER (NIN=5,NOUT=6) * .. Local Scalars .. DOUBLE PRECISION A, B, G, P INTEGER IFAIL * .. External Functions .. DOUBLE PRECISION G01EFF EXTERNAL G01EFF * .. Executable Statements .. WRITE (NOUT,*) 'G01EFF Example Program Results' WRITE (NOUT,*) WRITE (NOUT,*) + 'Gamma deviate Alpha Beta Lower tail prob.' WRITE (NOUT,*) * Skip heading in data file READ (NIN,*) 20 READ (NIN,*,END=40) G, A, B IFAIL = 1 * P = G01EFF('Lower',G,A,B,IFAIL) * IF (IFAIL.EQ.0) THEN WRITE (NOUT,99998) G, A, B, P GO TO 20 ELSE WRITE (NOUT,99999) ' ** G01EFF returned with IFAIL = ', IFAIL END IF 40 CONTINUE * 99999 FORMAT (1X,A,I5) 99998 FORMAT (1X,F9.2,F13.2,F9.2,F14.4) END