+

来自「此文挡是用Matlab开发的几种计算数值微积分的方法」· 代码 · 共 25 行

part1
function T=euler(f,a,b,ya,N)
h=(b-a)/N;
T=zeros(1,N+1);
Y=zeros(1,N+1);
T=a:h:b;
Y(1)=ya;
for j=1:N
    Y(j+1)=Y(j)+h*feval('f',T(j),Y(j));
end
T=[T' Y'];
***************************
part2
function y1=f(x,y)
y1=-y-x.*y.^2;
***************************
test
>> E=euler('f',0,0.6,1,3)

E =

         0    1.0000
    0.2000    0.8000
    0.4000    0.6144
    0.6000    0.4613