FUNCTION ASGAUSS(X0,SIG,k,j)
	DATA NTRY/2/
	NTRY=3-NTRY
	GO TO (1,2),NTRY
1	CONTINUE
	RR=RAN(K,J)
	R=SQRT(-2.*ALOG(RR))*SIG
	RR=RAN(K,J)
	ANG=6.283185308*RR
	ASGAUSS=X0+R*COS(ANG)
	RETURN
2	CONTINUE
	ASGAUSS=X0+R*SIN(ANG)
	RETURN
	END
	FUNCTION GAUSS(IX)
	IMPLICIT REAL*8 G
	k=ix
	j=ix+3
	x=ASGAUSS(0.0,1.0,k,j)
	ix=k
	gauss=x
	RETURN
	END
	function urand(ix)
	implicit real*8 u
	urand=ran(ix)
	return
	end