Mercurial > octave-nkf
view libcruft/ranlib/genchi.f @ 14193:72aebe619641 stable rc-3-6-0-0
3.6.0-rc0 release candidate
* configure.ac (AC_INIT): Version is now 3.6.0-rc0.
(OCTAVE_RELEASE_DATE): Now 2012-01-10.
(OCTAVE_API_VERSION_NUMBER): Now 47.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 10 Jan 2012 16:43:41 -0500 |
parents | 0e71ead7359d |
children |
line wrap: on
line source
REAL FUNCTION genchi(df) C********************************************************************** C C REAL FUNCTION GENCHI( DF ) C Generate random value of CHIsquare variable C C C Function C C C Generates random deviate from the distribution of a chisquare C with DF degrees of freedom random variable. C C C Arguments C C C DF --> Degrees of freedom of the chisquare C (Must be positive) C REAL DF C C C Method C C C Uses relation between chisquare and gamma. C C********************************************************************** C .. Scalar Arguments .. REAL df C .. C .. External Functions .. C REAL gengam C EXTERNAL gengam REAL sgamma EXTERNAL sgamma C .. C .. Executable Statements .. IF (.NOT. (df.LE.0.0)) GO TO 10 WRITE (*,*) 'DF <= 0 in GENCHI - ABORT' WRITE (*,*) 'Value of DF: ',df CALL XSTOPX ('DF <= 0 in GENCHI - ABORT') C JJV changed this to call sgamma directly C 10 genchi = 2.0*gengam(1.0,df/2.0) 10 genchi = 2.0*sgamma(df/2.0) RETURN END