my_ode_rough.m

来自「解偏微分方程」· M 代码 · 共 17 行

function my_ode_rough
% 一步Euler法求解

%   $[e-mail] $  
%   $Date: 2005/07/20    $

h   = 0.1;    
t   = 0:h:20
x(1)= -4;
for k=1:(length(t)-1)
    x(k+1)=x(k)+h*(-0.2*x(k)+2) ;
end

plot(t,x)