      SUBROUTINE GAUSS(IX,X)
C
C       THE ROUTINE GAUSS GENERATES A RANDOM NUMBER
C       IN A GAUSSIAN DISTRIBUTION, WITH A MEAN OF ZERO
C       AND A STANDARD DEVIATION OF ONE
C
C       NOTE: THIS ROUTINE DOES NOT USE THE LOG FEATURE
C
C       VAX/VMS VERSION
C
C       VARIABLES:
C
C       IX - INITIAL SEED INTEGER; GAUSS UPDATES IT WITH A NEW SEED VALUE
C
C       X - THE OUTPUT GAUSSIAN VARIABLE
C
      DOUBLE PRECISION   X
      DOUBLE PRECISION   Y
      DOUBLE PRECISION   SUM
      SUM=0.0
      DO 10 I=1,5
         CALL RAND(IX,IY,Y)
         IX = IY
         SUM=SUM+Y
10       CONTINUE
      X=6.0D0*(SUM/5.0D0)-3.0D0
      RETURN
      END
