MatLab实现改进的欧拉算法

分别取步长为0.1，0.05，0.01求解以上定解问题

代码如下:

function example_Euler
T = 1;
h = 0.05;
t = 0:h:T;
N = length(t)-1;
solu = exp(-5.0*t);
u0 = 1;
f = @f1;
u_euler = euler(f,u0,t,h,N);
u_implicit = implicit_euler(f,u0,t,h,N);
u_modified = update_euler(f,u0,t,h,N);
 
figure (1)
plot(t,u_euler,'*r',t,u_implicit,'o',t,u_modified,'d',t,solu,'r')
legend('显示欧拉法','隐式欧拉法','改进欧拉法','精确解')
end
 
function u = euler(f,u0,t,h,N)
u    = zeros(N+1,1);
u(1) = u0;
for n = 1:N
    fn     = f(t(n),u(n));
    u(n+1) = u(n)+h*fn;
end end
function u = implicit_euler(f,u0,t,h,N)
u    = zeros(N+1,1);
u(1) = u0;
 
eps_in = 1e-6;
  K_in = 6;
  
for n = 1:N
    s1 = u(n);
    du = 1;
    k  = 1;
    while abs(du)>eps_in && k<K_in
        s2 = u(n)+h*f(t(n+1),s1);
        du = s2-s1;
        s1 = s2;
         k = k+1;
    end
    u(n+1) = s2;end end
function u = update_euler(f,u0,t,h,N)
u    = zeros(N+1,1);
u(1) = u0;
eps_in = 1e-6;
  K_in = 6;
for n = 1:N
    s1 = u(n);
    du = 1;
    k  = 1;
    while abs(du)>eps_in && k<K_in
        s2 = u(n)+h/2*(f(t(n+1),s1)+f(t(n),u(n)));
        du = s2-s1;
        s1 = s2;
         k = k+1;end
    u(n+1) = s2;
end end

结果:分别对应三种步长

在h变小的过程中，改进欧拉方法的精度和误差逐渐降低，且均高于显示欧拉方法和隐式欧拉方法。