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