*     G05MRF Example Program Text
*     Mark 20 Release. NAG Copyright 2001.
*     .. Parameters ..
      INTEGER          K, NOUT, N, NR
      PARAMETER        (K=4,NOUT=6,N=20,NR=6007)
*     .. Local Scalars ..
      INTEGER          I, IFAIL, IGEN, J, M
*     .. Local Arrays ..
      DOUBLE PRECISION P(K), R(NR)
      INTEGER          ISEED(4), X(N,K)
*     .. External Subroutines ..
      EXTERNAL         G05KBF, G05MRF
*     .. Executable Statements ..
      WRITE (NOUT,*) 'G05MRF Example Program Results'
      WRITE (NOUT,*)
*     Set the distribution parameters P and M
      P(1) = 0.08D0
      P(2) = 0.1D0
      P(3) = 0.8D0
      P(4) = 0.02D0
      M = 6000
*     Initialise the seed to a repeatable sequence
      ISEED(1) = 1762543
      ISEED(2) = 9324783
      ISEED(3) = 42344
      ISEED(4) = 742355
*     IGEN identifies the stream.
      IGEN = 1
      CALL G05KBF(IGEN,ISEED)
*     Choose MODE = 2
      IFAIL = 0
      CALL G05MRF(2,M,K,P,N,X,N,IGEN,ISEED,R,NR,IFAIL)
*
      DO 20 I = 1, N
         WRITE (NOUT,99999) (X(I,J),J=1,K)
   20 CONTINUE
      STOP
*
99999 FORMAT (1X,10(1X,I12))
      END