期末作业,水博客
欧拉法
function result = Euler(f,y0,x0,x1,h)
F = matlabFunction(f);
l = (x1-x0)/h;
result = zeros(1,l);
for i=1:l
K1 = F(x0,y0);
y0 = y0 + h * K1;
x0 = x0 + h;
result(i) = y0;
end
改进欧拉法
function result = Euler_(f,y0,x0,x1,h)
F = matlabFunction(f);
l = (x1-x0)/h;
result = zeros(1,l);
for i=1:l
K1 = F(x0,y0);
K2 = F(x0+h, y0+K1*h);
y0 = y0 + h / 2 * (K1+K2);
x0 = x0 + h;
result(i) = y0;
end