FUNCTION TWERT(NF,OGR)
C INTEGRIERT DIE STUDENT-VERTEILUNG H(T) FUER DIE OBERE
C GRENZE OGR UND FREIHEITSGRADE NF>=4 MIT HILFE DER SIMPSON-REGEL
C FORTRAN II
	PI=3.1415927
C TEST OB NF>4
	IF(NF-4)1030,1002,1002
C TEST OB NF GERADE ODER UNGERADE
1002	A=FLOAT(NF)/2.
	IA=A
	B=A-FLOAT(IA)
	IF(B)1006,1006,1007
1006	A=1.
C A FUER NF GERADE
	DO 1008 I=4,NF,2
	FI=FLOAT(I)
1008	A=A*(FI-1.)/(FI-2.)
	A=.5*A
	GOTO 1010
1007	A=1.
C A FUER NF UNGERADE
	DO 1009 I=3,NF,2
	FI=FLOAT(I)
1009	A=A*(FI-1.)/(FI-2.)
	A=A/PI
1010	F=FLOAT(NF)
	C1=A/SQRT(F)
	C2=-.5*(F+1.)
	TWERT=0.
	XFT=0.
C INTEGRATIONS-SCHLEIFE-SCHRITTWEITE 0.1
1011	CONTINUE
	XA=-100.+XFT*1.E-1
	XM=XA+5.E-2
	XE=XA+1.E-1
	HA=C1*(1.+XA*(XA/F))**C2
	HM=C1*(1.+XM*(XM/F))**C2
	HE=C1*(1.+XE*(XE/F))**C2
	TWERT=TWERT+(HA+4.*HM+HE)
	XFT=XFT+1.
	IF(XE-OGR-5.E-6)1011,1012,1012
C INTEGRATION BEENDET
1012	RETURN
1030	WRITE(1,1031)
1031	FORMAT('BEDINGUNG NF>=4 NICHT ERFUELLT')
	GOTO 1035
1035	RETURN
	END