* G05NBF Example Program Text * Mark 20 Release. NAG Copyright 2001. * .. Parameters .. INTEGER NOUT, N PARAMETER (NOUT=6,N=8) * .. Local Scalars .. INTEGER I, IFAIL, IGEN, K, M * .. Local Arrays .. INTEGER IPOP(N), ISAMPL(N), ISEED(4) * .. External Subroutines .. EXTERNAL G05KBF, G05NBF * .. Executable Statements .. WRITE (NOUT,*) 'G05NBF Example Program Results' WRITE (NOUT,*) * Initialise the seed to a repeatable sequence ISEED(1) = 1762543 ISEED(2) = 9324783 ISEED(3) = 1542344 ISEED(4) = 742355 * IGEN identifies the stream. IGEN = 1 CALL G05KBF(IGEN,ISEED) WRITE (NOUT,99999) ' Samples from the first ', N, ' integers' WRITE (NOUT,*) WRITE (NOUT,*) ' Sample size Values' DO 20 I = 1, N IPOP(I) = I 20 CONTINUE DO 40 M = 1, N IFAIL = 0 CALL G05NBF(IPOP,N,ISAMPL,M,IGEN,ISEED,IFAIL) * WRITE (NOUT,99998) M, (ISAMPL(K),K=1,M) 40 CONTINUE STOP * 99999 FORMAT (1X,A,I1,A) 99998 FORMAT (1X,I6,10X,8I3) END