求助：把这个Matlab的Euler公式改进为Euler两步法

hynjtgw

function E2=Euler_2(fun,x0,y0,xN,N)
% Euler 向前公式，其中
% fun 为一阶微分方程的函数
% x0,y0 为初始条件
% xN 为取值范围的一个端点
% h 为区间步长
% N 为区间个数
% x 为 Xn 构成的向量
% y 为 Yn 构成的向量
x=zeros(1,N+1);y=zeros(1,N+1);
x(1)=x0;y(1)=y0;
h=(xN-x0)/N;
for n=1:N
    x(n+1)=x(n)+h;
    y(n+1)=y(n)+h*feval(fun,x(n),y(n));
end
T=[x',y']
新手上路 希望各位高手帮我解决问题。十分感谢！
另外也可以写出各位对Euler两步法编写的高见.再次感谢各位的帮助。



[ 本帖最后由 hynjtgw 于 2008-7-7 04:35 编辑 ]