
TRAPRULE(F,A,B,N):=
   (B-A)/N * ((F(B)+F(A))/2 + SUM(F(A+(B-A)/N*I),I,1,N-1))$

SIMPSON(F,A,B,N):=
   (B-A)/N/3 * (F(A)+F(B) + 4*SUM(F(A+(B-A)/N*(2*I-1)),I,1,N/2)
			  + 2*SUM(F(A+(B-A)/N*(2*I)),I,1,(N-2)/2))$
