File PHI.FT (FORTRAN source file)

Directory of image this file is from
This file as a plain text file

	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

Feel free to contact me, David Gesswein djg@pdp8online.com with any questions, comments on the web site, or if you have related equipment, documentation, software etc. you are willing to part with. I am interested in anything PDP-8 related, computers, peripherals used with them, DEC or third party, or documentation.

PDP-8 Home Page PDP-8 Site Map PDP-8 Site Search