gaijingeuler.m

来自「一些数值计算用的MATLAB代码,适合初学MATLAB和数字计算的人使用,根据M」· M 代码 · 共 13 行

a=input('x1=');
b=input('x2=');
h=input('h=');
y(1)=input('y(x1)=');
b=b-h;
i=1;
for x=a:h:b
    i=i+1;
k1=h*f(x,y(i-1));
k2=h*f(x+h,y(i-1)+k1);
y(i)=y(i-1)+(k1+k2)/2;
end
y