euler.m

来自「关于Euler、jacobi、rungkutta4数值算法的matlab实现」· M 代码 · 共 9 行

function[x,y]=euler(dyfun,xspan,y0,h)

x=xspan(1):h:xspan(2)
y(1)=y0;
for n=1:length(x)-1
    y(n+1)=y(n)+h*feval(dyfun,x(n),y(n));
end
x=x';y=y';