*     G13AJF Example Program Text
*     Mark 14 Revised. NAG Copyright 1989.
*     .. Parameters ..
      INTEGER          NPMAX, ISTMAX, IFVMAX, NXMAX, IW
      PARAMETER        (NPMAX=10,ISTMAX=10,IFVMAX=10,NXMAX=30,IW=250)
      INTEGER          NIN, NOUT
      PARAMETER        (NIN=5,NOUT=6)
*     .. Local Scalars ..
      DOUBLE PRECISION C, RMS
      INTEGER          I, IFAIL, IFV, IST, KFC, NFV, NPAR, NST, NX
*     .. Local Arrays ..
      DOUBLE PRECISION FSD(IFVMAX), FVA(IFVMAX), PAR(NPMAX), ST(ISTMAX),
     +                 W(IW), X(NXMAX)
      INTEGER          ISF(4), MR(7)
*     .. External Subroutines ..
      EXTERNAL         G13AJF
*     .. Intrinsic Functions ..
      INTRINSIC        MAX
*     .. Executable Statements ..
      WRITE (NOUT,*) 'G13AJF Example Program Results'
*     Skip heading in data file
      READ (NIN,*)
      READ (NIN,*) NX, (MR(I),I=1,7), NFV
      NPAR = MR(1) + MR(3) + MR(4) + MR(6)
      IF (NX.GT.0 .AND. NX.LE.NXMAX .AND. NPAR.GT.0 .AND. NPAR.LE.NPMAX)
     +    THEN
         READ (NIN,*) (X(I),I=1,NX)
         READ (NIN,*) (PAR(I),I=1,NPAR)
         READ (NIN,*) KFC, C
         IST = MR(4) + MR(7) + MR(2) + MR(5) + MR(3) + MAX(MR(1),MR(6)
     +         *MR(7))
         IFV = MAX(1,NFV)
         IFAIL = 0
*
         CALL G13AJF(MR,PAR,NPAR,C,KFC,X,NX,RMS,ST,IST,NST,NFV,FVA,FSD,
     +               IFV,ISF,W,IW,IFAIL)
*
         WRITE (NOUT,*)
         WRITE (NOUT,99999) 'The residual mean square is ', RMS
         WRITE (NOUT,*)
         WRITE (NOUT,99998) 'The state set consists of ', NST, ' values'
         WRITE (NOUT,99997) (ST(I),I=1,NST)
         WRITE (NOUT,*)
         WRITE (NOUT,99996) 'The ', NFV,
     +     ' forecast values and standard errors are -'
         WRITE (NOUT,99995) (FVA(I),FSD(I),I=1,NFV)
      END IF
*
99999 FORMAT (1X,A,F9.2)
99998 FORMAT (1X,A,I1,A)
99997 FORMAT (1X,4F11.4)
99996 FORMAT (1X,A,I2,A)
99995 FORMAT (10X,2F10.2)
      END