*     G05LZF Example Program Text
*     Mark 20 Release. NAG Copyright 2001.
*     .. Parameters ..
      INTEGER          NOUT, M, NR
      PARAMETER        (NOUT=6,M=2,NR=((M+1)*(M+2))/2)
*     .. Local Scalars ..
      INTEGER          IFAIL, IGEN
*     .. Local Arrays ..
      DOUBLE PRECISION C(M,M), R(NR), X(M), XMU(M)
      INTEGER          ISEED(4)
*     .. External Subroutines ..
      EXTERNAL         G05KBF, G05LZF
*     .. Executable Statements ..
      WRITE (NOUT,*) 'G05LZF Example Program Results'
      WRITE (NOUT,*)
*     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)
*
      C(1,1) = 2.0D0
      C(2,1) = 1.0D0
      C(1,2) = 1.0D0
      C(2,2) = 3.0D0
*
      XMU(1) = 1.0D0
      XMU(2) = 2.0D0

      IFAIL = 0

*     Set up reference vector and generate numbers
      CALL G05LZF(0,M,XMU,C,M,X,IGEN,ISEED,R,NR,IFAIL)
      WRITE (NOUT,99999) X

*     Generate numbers
      CALL G05LZF(2,M,XMU,C,M,X,IGEN,ISEED,R,NR,IFAIL)
      WRITE (NOUT,99999) X

      STOP
*
99999 FORMAT (1X,F10.4)
      END