RNUND

Generates pseudorandom numbers from a discrete uniform distribution.

Required Arguments

K — Parameter of the discrete uniform distribution. (Input)
The integers 1, 2, …, K occur with equal probability. K must be positive.

IR — Vector of length NR containing the random discrete uniform deviates. (Output)

Optional Arguments

NR — Number of random numbers to generate. (Input)
Default: NR = size (IR,1).

FORTRAN 90 Interface

Generic: CALL RNUND (K, IR [, …])

Specific: The specific interface name is S_RNUND.

FORTRAN 77 Interface

Single: CALL RNUND (NR, K, IR)

Description

Routine RNUND generates pseudorandom numbers from a discrete uniform distribution over the integers 1, 2, …, K. A random integer is generated by multiplying K by a uniform (0, 1) random number, adding 1.0, and truncating the result to an integer. This, of course, is equivalent to sampling with replacement from a finite population of size K. To do the equivalent of sampling without replacement, the routine RNSRI can be used.

Comments

The routine RNSET can be used to initialize the seed of the random number generator. The routine RNOPT can be used to select the form of the generator.

Example

In this example, RNUND is used to generate five pseudorandom deviates from a discrete uniform distribution over the integers from 1 to 6.

 

USE RNUND_INT

USE UMACH_INT

USE RNSET_INT

 

IMPLICIT NONE

INTEGER IR(5), ISEED, K, NOUT, NR

!

CALL UMACH (2, NOUT)

K = 6

NR = 5

ISEED = 123457

CALL RNSET (ISEED)

CALL RNUND (K, IR)

WRITE (NOUT,99999) IR

99999 FORMAT (' Discrete uniform (1,6) random deviates: ', 5I7)

END

Output

 

Discrete uniform (1,6) random deviates: 6 2 5 4 6