这是我的程序,即求变步长积分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:100h=(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

来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/06 05:42:56
这是我的程序,即求变步长积分functionT=Qvs(f,a,b,eps)t=zeors(50);T0=0.5*(b-a)*(feval(f,a)+feval(f,b));t(1)=0.5*T0+(

这是我的程序,即求变步长积分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:100h=(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
这是我的程序,即求变步长积分
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))

这是我的程序,即求变步长积分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:100h=(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
大致的意思是:变量 zeors 的数据类型有问题,你定义 zeors 用的什么类型?