Date:         Tue, 27 Feb 1996 01:43:43 GMT
Reply-To:     Chien-Hua Wu <cwu@MATH.UNL.EDU>
Sender:       "SAS(r) Discussion" <SAS-L@UGA.CC.UGA.EDU>
From:         Chien-Hua Wu <cwu@MATH.UNL.EDU>
Organization: University of Nebraska--Lincoln
Subject:      Help for Fortran Programming

I generated 10 different values based on Poisson distribution by Fortran IMSL.
The number of the uniform's values is based on the value we generated by
Poisson distribution.  The following is my program.  Please let me know, if
you find anything wrong with this program.  Thanks.
 
 
C
C       GENERATE N FROM POISSON THEN GENERATE N OBSERVATIONS FROM UNIFORM
C
        INTEGER NS, NR
        PARAMETER (NS=10)
C
        INTEGER IR(NS), ISEED, SOUT
        REAL THETA
        EXTERNAL RNPOI, RNSET, UMACH
C
        CALL UMACH (2, SOUT)
        THETA = 169
        ISEED = 123457
        CALL RNSET (ISEED)
        CALL RNPOI (NS, THETA, IR)
        PRINT *, IR
C
        DO 10, J = 1, NS
                NR = IR(J)
        INTEGER ISEED, NOUT
        EXTERNAL RNSET, RNUN, UMACH
C
        CALL UMACH (2, NOUT)
        ISEED = 123457
        REAL R
        CALL RNSET (ISEED)
        CALL RNUN (NR, R)
10      CONTINUE
C
        PRINT *, R
        STOP
        END