FUNCTION GAUSS(OGR)
C BERECHNET DAS GAUSS'SCHE FEHLERINTEGRAL FUER
C DIE OBERE GRENZE OGR MIT HILFE DER SIMPSON-REGEL
C FORTRAN II
	ID=0
	IF(OGR)1012,1010,1011
1010	GAUSS=.5
	RETURN
1012	ID=1
	OGR=ABS(OGR)
1011	IF(OGR-4.)1015,1015,1014
1014	GFI=.5
	GOTO 1018
1015	C=6.6490383E-4
	XFI=0.
	GFI=0.
C INTEGRATIONSSCHLEIFE - SCHRITTWEITE 0.01
1016	CONTINUE
	XA=XFI*1.E-2
	XM=XA+5.E-3
	XE=XA+1.E-2
	YA=EXP(-.5*XA*XA)
	YM=EXP(-.5*XM*XM)
	YE=EXP(-.5*XE*XE)
	GFI=GFI+(YA+4.*YM+YE)
	XFI=XFI+1.
	IF(XE-OGR+5.E-3)1016,1017,1017
C INTEGRATION BEENDET
1017	GFI=C*GFI
1018	IF(ID)1020,1019,1020
1019	GAUSS=.5+GFI
	RETURN
1020	GAUSS=.5-GFI
	RETURN
	END