function T=rctrap(fun,a,b,m)
n=1;h=b-a; T=zeros(1,m+1); x=a; T(1)=h*(feval(fun,a)+feval(fun,b))/2;
for i=1:m
           h=h/2; n=2*n; s=0;
          for k=1:n/2
           x=a+h*(2*k-1); s=s+feval(fun,x);
end
T(i+1)=T(i)/2+h*s;
end                    
T=T(1+m)


T=rctrap(fun, 0,pi/2,14), syms t
fi=int(exp((-t^2)/2)/(sqrt(2*pi)),t,0, pi/2);
Fs= double(fi), wT= double(abs(fi-T))
这两段程序怎样才能运行出来啊？ 我是matlab的初学者，迫切等指导，谢谢啊