Next: , Previous: DRANDSTUDENTST, Up: Continuous Univariate Distributions


DRANDTRIANGULAR / SRANDTRIANGULAR

Generates a vector of random variates from a Triangular distribution with probability density function, f(X), where: f(X) = [2 * (X - XMIN)] / [(XMAX - XMIN) * (XMED - XMIN)], if XMIN < X <= XMED, else f(X) = [2 * (XMAX - X)] / [(XMAX - XMIN) * (XMAX - XMED)], if XMED < X <= XMAX, otherwise f(X) = 0.

(Note that SRANDTRIANGULAR is the single precision version of DRANDTRIANGULAR. The argument lists of both routines are identical except that any double precision arguments of DRANDTRIANGULAR are replaced in SRANDTRIANGULAR by single precision arguments - type REAL in FORTRAN or type float in C).

— SUBROUTINE: DRANDTRIANGULAR (N,XMIN,XMED,XMAX,STATE,X,INFO)
— Input: INTEGER N

On input: number of variates required.
Constraint: N>=0.

— Input: DOUBLE PRECISION XMIN

On input: minimum value for the distribution.

— Input: DOUBLE PRECISION XMED

On input: median value for the distribution.
Constraint: XMIN<=XMED <=XMAX.

— Input: DOUBLE PRECISION XMAX

On input: maximum value for the distribution.
Constraint: XMAX>=XMIN.

— Input/Output: INTEGER STATE(*)

The STATE vector holds information on the state of the base generator being used and as such its minimum length varies. Prior to calling DRANDTRIANGULAR STATE must have been initialized. See Initialization of the Base Generators for information on initialization of the STATE variable.
On input: the current state of the base generator.
On output: the updated state of the base generator.

— Output: DOUBLE PRECISION X(N)

On output: vector of variates from the specified distribution.

— Output: INTEGER INFO

On output: INFO is an error indicator. On successful exit, INFO contains 0. If INFO = -i on exit, the i-th argument had an illegal value.

     Example:
     

     Generate 100 values from the Triangular distribution
                 INTEGER LSTATE,N
                 PARAMETER (LSTATE=16,N=100)
                 INTEGER I,INFO,SEED(1),STATE(LSTATE)
                 DOUBLE PRECISION XMIN,XMAX,XMED
                 DOUBLE PRECISION X(N)
          C      Set the seed
                 SEED(1) = 1234
          
          C      Read in the distributional parameters
                 READ(5,*) XMIN,XMAX,XMED
          
          C      Initialize the STATE vector
                 CALL DRANDINITIALIZE(1,1,SEED,1,STATE,LSTATE,INFO)
          C      Generate N variates from the Triangular distribution
                 CALL DRANDTRIANGULAR(N,XMIN,XMAX,XMED,STATE,X,INFO)
          
          C      Print the results
                 WRITE(6,*) (X(I),I=1,N)