* E02RBF Example Program Text * Mark 14 Revised. NAG Copyright 1989. * .. Parameters .. INTEGER L, M, IA, IB, IC, IW PARAMETER (L=4,M=4,IA=L+1,IB=M+1,IC=IA+IB-1,IW=IB*(2*IB+3)) INTEGER NOUT PARAMETER (NOUT=6) * .. Local Scalars .. DOUBLE PRECISION ANS, TVAL, X INTEGER I, IFAIL * .. Local Arrays .. DOUBLE PRECISION A(IA), B(IB), CC(IC), W(IW) * .. External Subroutines .. EXTERNAL E02RAF, E02RBF * .. Intrinsic Functions .. INTRINSIC DBLE, EXP * .. Executable Statements .. WRITE (NOUT,*) 'E02RBF Example Program Results' CC(1) = 1.0D0 DO 20 I = 1, IC - 1 CC(I+1) = CC(I)/DBLE(I) 20 CONTINUE IFAIL = 1 * CALL E02RAF(IA,IB,CC,IC,A,B,W,IW,IFAIL) * WRITE (NOUT,*) IF (IFAIL.EQ.0) THEN WRITE (NOUT,*) ' X Pade True' DO 40 I = 1, 10 X = DBLE(I)/10.0D0 IFAIL = 0 * CALL E02RBF(A,IA,B,IB,X,ANS,IFAIL) * TVAL = EXP(X) WRITE (NOUT,99999) X, ANS, TVAL 40 CONTINUE ELSE WRITE (NOUT,99998) ' ** E02RAF returned with IFAIL = ', IFAIL END IF * 99999 FORMAT (1X,F6.1,3E15.5) 99998 FORMAT (1X,A,I5) END