File UGFI.FT (FORTRAN source file)

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

	SUBROUTINE UGFI(P,U,IE)
C BERECHNET U=F(P), DIE UMKEHRFUNKTION DES GAUSS'SCHEN
C FEHLER-INTEGRALS MIT EINEM NAEHERUNGS-POLYNOM
C	IE=0 : KEIN FEHLER
C	IE=1 : P NICHT IM INTERVALL <0,1>
C MAXIMALER FEHLER 0.00045
C FORTRAN II
	IE=0
	IF(P)1001,1004,1002
1001	IE=1
	GOTO 1012
1002	IF(P-1.0)1007,1006,1001
1004	U=-.999999E+38
	GOTO 1012
1006	U=+.999999E+38
	GOTO 1012
1007	TEM=P
	IF(TEM-0.5)1009,1009,1008
1008	TEM=1.0-TEM
1009	TEM=ALOG(1.0/(TEM*TEM))
	T=SQRT(TEM)
	U=2.515517+0.802853*T+0.010328*TEM
	U=T-U/(1.0+1.432788*T+0.189269*TEM+0.001308*T*TEM)
	IF(P-0.5)1010,1010,1012
1010	U=-U
1012	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