* G01ASF Example Program Text * Mark 14 Release. NAG Copyright 1989. * .. Parameters .. INTEGER NIN, NOUT PARAMETER (NIN=5,NOUT=6) INTEGER NMAX, LDX PARAMETER (NMAX=5,LDX=10) * .. Local Scalars .. INTEGER I, IFAIL, J, LDPLOT, M, NSTEPX, NSTEPY, OUTCHN * .. Local Arrays .. DOUBLE PRECISION WORK(5*NMAX), X(LDX,NMAX) INTEGER IWORK(LDX), N(NMAX) CHARACTER PLOT(60,132) * .. External Subroutines .. EXTERNAL G01ASF, X04ABF * .. Executable Statements .. WRITE (NOUT,*) 'G01ASF Example Program Results' OUTCHN = NOUT * Skip heading in data file READ (NIN,*) * Set advisory message unit for plot output to OUTCHN CALL X04ABF(1,OUTCHN) READ (NIN,*) M, (N(I),I=1,M) READ (NIN,*) NSTEPX, NSTEPY DO 20 J = 1, M READ (NIN,*) (X(I,J),I=1,N(J)) 20 CONTINUE LDPLOT = NSTEPY IFAIL = 1 WRITE (NOUT,*) * CALL G01ASF('Print',M,N,X,LDX,NSTEPX,NSTEPY,PLOT,LDPLOT,WORK, + IWORK,IFAIL) * IF (IFAIL.NE.0) THEN WRITE (NOUT,*) WRITE (NOUT,99999) ' ** G01ASF returned with IFAIL = ', IFAIL END IF * 99999 FORMAT (1X,A,I5) END