* G01DCF Example Program Text * Mark 14 Revised. NAG Copyright 1989. * .. Parameters .. INTEGER N, IW PARAMETER (N=6,IW=3*N/2) INTEGER NOUT PARAMETER (NOUT=6) * .. Local Scalars .. DOUBLE PRECISION ERREST, ETOL, EXP1, EXP2, SUMSSQ INTEGER I, IFAIL, J, K * .. Local Arrays .. DOUBLE PRECISION PP(N), VEC(N*(N+1)/2), WORK(IW) * .. External Subroutines .. EXTERNAL G01DAF, G01DCF * .. Executable Statements .. WRITE (NOUT,*) 'G01DCF Example Program Results' ETOL = 0.0001D0 IFAIL = 1 * CALL G01DAF(N,PP,ETOL,ERREST,WORK,IW,IFAIL) * IF (IFAIL.EQ.0) THEN EXP1 = PP(N) EXP2 = PP(N-1) SUMSSQ = 0.0D0 DO 20 I = 1, N SUMSSQ = SUMSSQ + PP(I)*PP(I) 20 CONTINUE IFAIL = 1 * CALL G01DCF(N,EXP1,EXP2,SUMSSQ,VEC,IFAIL) * IF (IFAIL.EQ.0) THEN WRITE (NOUT,*) WRITE (NOUT,99999) 'Sample size = ', N WRITE (NOUT,*) 'Variance-covariance matrix' K = 1 DO 40 J = 1, N WRITE (NOUT,99998) (VEC(I),I=K,K+J-1) K = K + J 40 CONTINUE ELSE WRITE (NOUT,*) WRITE (NOUT,99997) ' ** G01DCF returned with IFAIL = ', + IFAIL END IF ELSE WRITE (NOUT,*) WRITE (NOUT,99997) ' ** G01DAF returned with IFAIL = ', IFAIL END IF * 99999 FORMAT (1X,A,I2) 99998 FORMAT (1X,6F8.4) 99997 FORMAT (1X,A,I5) END