FUNCTION SIMPSON(N,H,F)
C	N=(UNGERADE) ANZAHL DER AEQUIDISTANTEN
C	STUETZSTELLEN (N GROESSER 3)
C	H=SCHRITTWEITE
C	F=FELD DER FUNKTIONSWERTE AN DEN STUETZSTELLEN
	DIMENSION F(N)
2	SIMPSON=F(1)+F(N)+4*F(N-1)
	M=N-3
	DO 8 I=2,M,2
8	SIMPSON=SIMPSON+2*F(I+1)+4(F(I)
	SIMPSON=SIMPSON*H/3.
	RETURN
	END