前向欧拉法、后向欧拉法、梯形方法、改进欧拉方法MATLAB

采用下述方法，求解常微分方程初值问题 y’=y-2x/y，y(0)=1，计算区间为[0, 1]， 步 长为 0.1。 (1)前向欧拉法。 (2)后向欧拉法。 (3)梯形方法。 (4)改进欧拉方法。

function  [x y] = BackwordEuler（f y0 a b h）
y（1） = y0;
n = （b - a）/h;
x = a:h:b;
for i=1:n
    yt = y（i） + h * feval（f x（i） y（i））;
    done = 0;
    while  ~done
        y （ i + 1 ） = y （ i ） + h * feval （ f x （ i+1 ） yt ）;
        done = （ abs （ y （ i + 1 ） - yt ） < 1e-6 ）;
        yt = y （ i + 1 ）;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-03-18 00:45  6\
     文件         334  2018-07-05 19:36  6\BackwordEuler.m
     文件         187  2018-07-05 03:01  6\ForwardEuler.m
     文件         382  2018-07-05 03:23  6\Trapezium.m
     文件         256  2018-07-05 03:28  6\ImprovedEuler.m