FUNCTION PHI(Z)
C VERTEILUNGSFUNKTION DER NORMIERTEN NORMALVERTEILUNG
C STANDARDABWEICHUNG Z IST ARGUMENT
C FORTRAN II
C NACH DECUS #8-705
	TOL=0.1E-6
	TWO=2.
	ZSQ=Z*Z/TWO
	PHI=0.
	IF(ZSQ-8.0)1004,1004,1002
1002	IF(Z)1020,1003,1003
1003	PHI=1.
	GOTO 1020
1004	CON=1.1283792
	ZPR=ZSQ*ZSQ
	ZPO=0.70710636*Z
	FOU=4.
	FAC=1.
	A1=FAC
	A2=FAC
	A3=FAC
1005	H1=A1*(A2+TWO)-A2*ZSQ
	H2=A2*(A2+A3)*(A2+TWO)
	DEL=(ZPO*H1)/(FAC*H2)
	PHI=PHI+DEL
	IF(DEL)1011,1012,1012
1011	DEL=-DEL
1012	IF(DEL-TOL)1015,1015,1010
1010	ZPO=ZPO*ZPR
	FAC=FAC*A1*(A1+A3)
	A1=A1+TWO
	A2=A2+FOU
	GOTO 1005
1015	PHI=0.5+PHI*CON
1020	RETURN
	END