5 * Revision 1.1.1.1 1996/04/01 15:02:42 mclareni
10 #if defined(CERNLIB_DOUBLE)
13 C Computes a "Normal Deviate"
14 C Based on G.W. Hill & A.W. Davis, Algorithm 442 Normal Deviate
15 C Collected Algorithms from CACM
19 PARAMETER (NAME = 'GAUSIN')
20 PARAMETER (C = 2.50662 827)
22 IF(P .LE. 0 .OR. P .GE. 1) THEN
25 CALL MTLPRT(NAME,'G105.1',ERRTXT)
26 ELSEIF(P .EQ. 0.5) THEN
32 X=X-((7.45551*X+450.636)*X+1271.059)/
33 1 (((X+110.4212)*X+750.365)*X+500.756)
35 H=C*(P-FREQ(X))*EXP(0.5*X**2)+X
39 101 FORMAT('ARGUMENT P =',1P,E15.5,' NOT IN RANGE')