euler_f.m

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function [x,y]=Euler_f(ydot_fun,x0,y0,h,N)
x=zeros(1,N+1);y=zeros(1,N+1);
x(1)=x0;y(1)=y0;
for n=1:N
    x(n+1)=x(n)+h;
    y(n+1)=y(n)+h*feval(ydot_fun,x(n),y(n));
end