matlab高手请进````是Runge-Kutta算法的请问怎么做,下面的代码总提示有错````

Question

h=0.05;
a=0;
b=0.5;
x=[a:h:b];
y(1)=y1;***************************************提示这里有错
n=(b-a)/h+1;
for j=2:n
fk1=f(x(i-1),y(i-1));
fk2=f(x(i-1)+h/2,y(i-1)+fk1*h/2);
fk3=f(x(i-1)+h/2,y(i-1)+fk2*h/2);
fk4=f(x(i-1)+h,y(i-1)+fk3*h);
y(i)=y(i-1)+h*(fk1+2*fk2+2*fk3+fk4)/6;
end

另一个文件:文件名是f.m
function z=f(x,y)
z=(-y+x^2+4*x-1)/2;
end

Answer 1

你那行
y(1)=y1直接改成初值。
即x=a的时候，y的值

然后将你的for j=2:n
改成for i=2:n
应当就可以了。