Procedures

Sample program:

Sample program:

ISEED An integer iseed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

ISEED An integer iseed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

CDF The cumulative distribution function value for the binomial distribution.

CDF The cumulative distribution function value for the binomial distribution.

PPF The percent point function value.

PPF The percent point function value.

N The desired integer number of random numbers to be generated.

N The desired integer number of random numbers to be generated.

CDF The cumulative distribution function value.

CDF The cumulative distribution function value.

PDF The probability density function value.

PDF The probability density function value.

Sample program:

Sample program:

Sample program:

Sample program:

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

–X A vector (of dimension at least N) into which the generated random sample of size N from the chi-squared distribution will be placed.

–X A vector (of dimension at least N) into which the generated random sample of size N from the chi-squared distribution will be placed.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

(On the UNIVAC 1108, EXEC 8 system at NBS)

(On the UNIVAC 1108, EXEC 8 system at NBS)

Sample program:

Sample program:

Sample program:

Sample program:

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

ISEED An integer iseed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

ISEED An integer iseed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

ISEED An integer seed value. Should be set to a non-negative value to start a new sequence of values. Will be set to -1 on return to indicate the next call should continue the current random sequence walk.

Sample program:

Sample program:

Sample program:

Sample program:

X The vector of (unsorted or sorted) observations.

X The vector of (unsorted or sorted) observations.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

X The vector of (unsorted or sorted) observations.

X The vector of (unsorted or sorted) observations.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

! INPUT ARGUMENTS–Y = THE VECTOR OF ! (UNSORTED OR SORTED) OBSERVATIONS ! TO BE PLOTTED VERTICALLY. ! –X = THE VECTOR OF ! (UNSORTED OR SORTED) OBSERVATIONS ! TO BE PLOTTED HORIZONTALLY. ! –N = THE INTEGER NUMBER OF OBSERVATIONS ! IN THE VECTOR Y. ! –D = THE VECTOR ! WHICH ‘DEFINES’ THE VARIOUS ! POSSIBLE SUBSETS. ! –DMIN = THE VALUE ! WHICH DEFINES THE LOWER BOUND ! (INCLUSIVELY) OF THE PARTICULAR ! SUBSET OF INTEREST TO BE PLOTTED. ! –DMAX = THE VALUE ! WHICH DEFINES THE UPPER BOUND ! (INCLUSIVELY) OF THE PARTICULAR ! SUBSET OF INTEREST TO BE PLOTTED. ! OUTPUT–A NARROW-WIDTH (71-CHARACTER) TERMINAL PLOT ! OF Y(I) VERSUS X(I), ! FOR ONLY OF A SPECIFIED SUBSET OF THE DATA. ! THE BODY OF THE PLOT (NOT COUNTING AXIS VALUES ! AND MARGINS) IS 25 ROWS (LINES) AND 49 COLUMNS.

! INPUT ARGUMENTS–Y = THE VECTOR OF ! (UNSORTED OR SORTED) OBSERVATIONS ! TO BE PLOTTED VERTICALLY. ! –X = THE VECTOR OF ! (UNSORTED OR SORTED) OBSERVATIONS ! TO BE PLOTTED HORIZONTALLY. ! –N = THE INTEGER NUMBER OF OBSERVATIONS ! IN THE VECTOR Y. ! –D = THE VECTOR ! WHICH ‘DEFINES’ THE VARIOUS ! POSSIBLE SUBSETS. ! –DMIN = THE VALUE ! WHICH DEFINES THE LOWER BOUND ! (INCLUSIVELY) OF THE PARTICULAR ! SUBSET OF INTEREST TO BE PLOTTED. ! –DMAX = THE VALUE ! WHICH DEFINES THE UPPER BOUND ! (INCLUSIVELY) OF THE PARTICULAR ! SUBSET OF INTEREST TO BE PLOTTED. ! OUTPUT–A NARROW-WIDTH (71-CHARACTER) TERMINAL PLOT ! OF Y(I) VERSUS X(I), ! FOR ONLY OF A SPECIFIED SUBSET OF THE DATA. ! THE BODY OF THE PLOT (NOT COUNTING AXIS VALUES ! AND MARGINS) IS 25 ROWS (LINES) AND 49 COLUMNS.

Y The vector of (unsorted or sorted) observations to be plotted vertically.

Y The vector of (unsorted or sorted) observations to be plotted vertically.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

The rank of the first element of the vector X will be placed in the first position of the vector XR, the rank of the second element of the vector X will be placed in the second position of the vector XR, etc.

The rank of the first element of the vector X will be placed in the first position of the vector XR, the rank of the second element of the vector X will be placed in the second position of the vector XR, etc.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

SORTP(3f) sorts (in ascending order) the N elements of the precision precision vector X, puts the resulting N sorted values into the precision precision vector Y; and puts the position (in the original vector X) of each of the sorted values into the REAL vector XPOS.

SORTP(3f) sorts (in ascending order) the N elements of the precision precision vector X, puts the resulting N sorted values into the precision precision vector Y; and puts the position (in the original vector X) of each of the sorted values into the REAL vector XPOS.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

call unit_check(‘autoco’, 0.eq.0, ‘checking’,100)

call unit_check(‘betran’, 0.eq.0, ‘checking’,100)

call unit_check(‘bincdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘binppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘binran’, 0.eq.0, ‘checking’,100)

call unit_check(‘caucdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘caupdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘cauplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘cauppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘cauran’, 0.eq.0, ‘checking’,100)

call unit_check(‘causf’, 0.eq.0, ‘checking’,100)

call unit_check(‘chscdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘chsplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘chsppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘chsran’, 0.eq.0, ‘checking’,100)

call unit_check(‘code’, 0.eq.0, ‘checking’,100)

call unit_check(‘corr’, 0.eq.0, ‘checking’,100)

call unit_check(‘count’, 0.eq.0, ‘checking’,100)

call unit_check(‘decomp’, 0.eq.0, ‘checking’,100)

call unit_check(‘define’, 0.eq.0, ‘checking’,100)

call unit_check(‘delete’, 0.eq.0, ‘checking’,100)

call unit_check(‘demod’, 0.eq.0, ‘checking’,100)

call unit_check(‘dexcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘dexpdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘dexplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘dexppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘dexran’, 0.eq.0, ‘checking’,100)

call unit_check(‘dexsf’, 0.eq.0, ‘checking’,100)

call unit_check(‘discr2’, 0.eq.0, ‘checking’,100)

call unit_check(‘discr3’, 0.eq.0, ‘checking’,100)

call unit_check(‘discre’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev1cdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev1plt’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev1ppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev1ran’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev2cdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev2plt’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev2ppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘ev2ran’, 0.eq.0, ‘checking’,100)

call unit_check(‘expcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘exppdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘expplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘expppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘expran’, 0.eq.0, ‘checking’,100)

call unit_check(‘expsf’, 0.eq.0, ‘checking’,100)

call unit_check(‘extrem’, 0.eq.0, ‘checking’,100)

call unit_check(‘fcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘fourie’, 0.eq.0, ‘checking’,100)

call unit_check(‘fran’, 0.eq.0, ‘checking’,100)

call unit_check(‘freq’, 0.eq.0, ‘checking’,100)

call unit_check(‘gamcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘gamplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘gamppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘gamran’, 0.eq.0, ‘checking’,100)

call unit_check(‘geocdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘geoplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘geoppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘georan’, 0.eq.0, ‘checking’,100)

call unit_check(‘hfncdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘hfnplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘hfnppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘hfnran’, 0.eq.0, ‘checking’,100)

call unit_check(‘hist’, 0.eq.0, ‘checking’,100)

call unit_check(‘invxwx’, 0.eq.0, ‘checking’,100)

call unit_check(‘lamcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘lampdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘lamplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘lamppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘lamran’, 0.eq.0, ‘checking’,100)

call unit_check(‘lamsf’, 0.eq.0, ‘checking’,100)

call unit_check(‘lgncdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘lgnplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘lgnppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘lgnran’, 0.eq.0, ‘checking’,100)

call unit_check(‘loc’, 0.eq.0, ‘checking’,100)

call unit_check(‘logcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘logpdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘logplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘logppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘logran’, 0.eq.0, ‘checking’,100)

call unit_check(‘logsf’, 0.eq.0, ‘checking’,100)

call unit_check(‘median’, 0.eq.0, ‘checking’,100)

call unit_check(‘move’, 0.eq.0, ‘checking’,100)

call unit_check(‘nbcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘nbppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘nbran’, 0.eq.0, ‘checking’,100)

call unit_check(‘norcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘norout’, 0.eq.0, ‘checking’,100)

call unit_check(‘norpdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘norplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘norppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘norran’, 0.eq.0, ‘checking’,100)

call unit_check(‘norsf’, 0.eq.0, ‘checking’,100)

call unit_check(‘parcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘parplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘parppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘parran’, 0.eq.0, ‘checking’,100)

call unit_check(‘plot’, 0.eq.0, ‘checking’,100)

call unit_check(‘plot10’, 0.eq.0, ‘checking’,100)

call unit_check(‘plot6’, 0.eq.0, ‘checking’,100)

call unit_check(‘plot7’, 0.eq.0, ‘checking’,100)

call unit_check(‘plot8’, 0.eq.0, ‘checking’,100)

call unit_check(‘plot9’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotc’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotco’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotct’, 0.eq.0, ‘checking’,100)

call unit_check(‘plots’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotsc’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotsp’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotst’, 0.eq.0, ‘checking’,100)

call unit_check(‘plott’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotu’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotx’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotxt’, 0.eq.0, ‘checking’,100)

call unit_check(‘plotxx’, 0.eq.0, ‘checking’,100)

call unit_check(‘pltsct’, 0.eq.0, ‘checking’,100)

call unit_check(‘pltxxt’, 0.eq.0, ‘checking’,100)

call unit_check(‘poicdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘poiplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘poippf’, 0.eq.0, ‘checking’,100)

call unit_check(‘poiran’, 0.eq.0, ‘checking’,100)

call unit_check(‘poly’, 0.eq.0, ‘checking’,100)

call unit_check(‘propor’, 0.eq.0, ‘checking’,100)

call unit_check(‘range’, 0.eq.0, ‘checking’,100)

call unit_check(‘rank’, 0.eq.0, ‘checking’,100)

call unit_check(‘ranper’, 0.eq.0, ‘checking’,100)

call unit_check(‘read’, 0.eq.0, ‘checking’,100)

call unit_check(‘readg’, 0.eq.0, ‘checking’,100)

call unit_check(‘relsd’, 0.eq.0, ‘checking’,100)

call unit_check(‘replac’, 0.eq.0, ‘checking’,100)

call unit_check(‘retain’, 0.eq.0, ‘checking’,100)

call unit_check(‘runs’, 0.eq.0, ‘checking’,100)

call unit_check(‘sampp’, 0.eq.0, ‘checking’,100)

call unit_check(‘scale’, 0.eq.0, ‘checking’,100)

call unit_check(‘sd’, 0.eq.0, ‘checking’,100)

call unit_check(‘skipr’, 0.eq.0, ‘checking’,100)

call unit_check(‘sortp’, 0.eq.0, ‘checking’,100)

call unit_check(‘spcorr’, 0.eq.0, ‘checking’,100)

call unit_check(‘stmom3’, 0.eq.0, ‘checking’,100)

call unit_check(‘stmom4’, 0.eq.0, ‘checking’,100)

call unit_check(‘subse1’, 0.eq.0, ‘checking’,100)

call unit_check(‘subse2’, 0.eq.0, ‘checking’,100)

call unit_check(‘subset’, 0.eq.0, ‘checking’,100)

call unit_check(‘tail’, 0.eq.0, ‘checking’,100)

call unit_check(‘tcdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘time’, 0.eq.0, ‘checking’,100)

call unit_check(‘tol’, 0.eq.0, ‘checking’,100)

call unit_check(‘tplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘tppf’, 0.eq.0, ‘checking’,100)

call unit_check(‘tran’, 0.eq.0, ‘checking’,100)

call unit_check(‘trim’, 0.eq.0, ‘checking’,100)

call unit_check(‘unimed’, 0.eq.0, ‘checking’,100)

call unit_check(‘unipdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘uniplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘unippf’, 0.eq.0, ‘checking’,100)

call unit_check(‘uniran’, 0.eq.0, ‘checking’,100)

call unit_check(‘unisf’, 0.eq.0, ‘checking’,100)

call unit_check(‘weib’, 0.eq.0, ‘checking’,100)

call unit_check(‘weicdf’, 0.eq.0, ‘checking’,100)

call unit_check(‘weiplt’, 0.eq.0, ‘checking’,100)

call unit_check(‘weippf’, 0.eq.0, ‘checking’,100)

call unit_check(‘weiran’, 0.eq.0, ‘checking’,100)

call unit_check(‘wind’, 0.eq.0, ‘checking’,100)

call unit_check(‘write’, 0.eq.0, ‘checking’,100)

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

CDF the REAL cumulative distribution function value.

CDF the REAL cumulative distribution function value.

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program:

Sample program: