function s=simpr1(f1432,a,b,n)
%  f是被積函數(shù)。
%  a, b分別為積分的上下限。
%  n是子區(qū)間的個(gè)數(shù)。
%  s是梯形總面積。
h=(b-a)/(2*n);
s1=0;
s2=0;
for k=1:n
    x=a+h*(2*k-1);
    s1=s1+feval('f1432',x);
end
for k=1:(n-1);
    x=a+h*2*k;
    s2=s2+feval('f1432',x);
end
s=h*(feval('f1432',a)+feval('f1432',b)+4*s1+2*s2)/3;