GEANT321/ghrout/svgam7.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:12  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 19/05/94  14.28.55  by  S.Ravndal
  11 *-- Author :
  12       SUBROUTINE SVGAM7(A,X)
  13       DIMENSION U(2)
  14       REAL LA
  15       GA=A-1.
  16       LA=SQRT(2.*A-1.)
  17       EP=1.570796327+ATAN(GA/LA)
  18       RO=1.570796327-EP
  19     1 CALL GRNDM(U,2)
  20       X=GA+LA*TAN(EP*U(1)+RO)
  21       IF (X.EQ.0) GOTO 1
  22       Y=LOG(1.+((X-GA)/LA)**2) +GA*LOG(X/GA)-X+GA
  23       IF(LOG(U(2)).GT.Y) GOTO 1
  24       RETURN
  25       END