function s=tixing(a,b,n)
% a，b为积分区间
% n为区间划分细度
h=(b-a)/n;
s=0;
for k=1:(n-1)
   x= a + h*k;
   s=s + fun(x);
end
s=h*(fun(a)+fun(b))/2+h*s;