* F01ECF Example Program Text * Mark 22 Release. NAG Copyright 2008. * .. Parameters .. INTEGER NIN, NOUT PARAMETER (NIN=5,NOUT=6) INTEGER NMAX PARAMETER (NMAX=10) INTEGER LDA PARAMETER (LDA=NMAX) * .. Local Scalars .. INTEGER I, IERR, IFAIL, J, N * .. Local Arrays .. DOUBLE PRECISION A(LDA,NMAX) * .. External Subroutines .. EXTERNAL F01ECF, X04CAF * .. Executable Statements .. WRITE (NOUT,*) 'F01ECF Example Program Results' WRITE (NOUT,*) * Skip heading in data file READ (NIN,*) READ (NIN,*) N IF (N.LE.NMAX) THEN * * Read A from data file * READ (NIN,*) ((A(I,J),J=1,N),I=1,N) * * Find exp( A ) * IFAIL = -1 CALL F01ECF(N,A,LDA,IFAIL) * IF (IFAIL.EQ.0) THEN * * Print solution * IERR = 0 CALL X04CAF('General',' ',N,N,A,LDA,'Exp(A)',IERR) * ELSE * WRITE (NOUT,*) WRITE (NOUT,*) 'Failure in F01ECF' WRITE (NOUT,99999) IFAIL END IF ELSE WRITE (NOUT,*) 'NMAX too small' END IF * 99999 FORMAT (6X,'IFAIL = ',I5) END