这是我的程序,即求变步长积分

题目

这是我的程序,即求变步长积分
function T = Qvs(f,a,b,eps)
t=zeors(50);
T0=0.5*(b-a)*(feval(f,a)+feval(f,b));
t(1)=0.5*T0+(b-a)*(feval(f,a+(b-a)/2))/2;
for K=2:100
h=(b-a)/(2^K);
y=zeros(n);
for i=1:2^(K-1)
y(i)=feval(f,a+(2*i-1)*h);
end
t(K)=0.5*t(K-1)+h*sum(y);
if abs(t(K)-t(K-1))

答案