euler.m

来自「把微分方程的一维经典初值问题的有限插分法：Euler法、改进的Euler法、Ad」· M 代码 · 共 10 行

%euler法解一维微分方程的通用子函数
%函数的输入变量为：T为x的取值上限，n为把[0，T]区间等份的数目，u0为方程的初始值
function u=euler(T,n,u0)
h=T/n;
u=[];
u(1)=u0;
for i=1:n-1
    u(i+1)=u(i)+h*f(h*i,u(i));
end