Algorithme de la Méthode d'Euler pour résoudre les équations différentielles de premier ordre en Matlab

Algorithme de la Méthode d'Euler pour résoudre les équations différentielles de premier ordre en Matlab :

%nm610: Euler method to solve a 1st-order differential equation clear, clf a = 1; r = 1; y0 = 0; tf = 2; t = [0:0.01:tf]; yt=1- exp(-a*t); %Eq.(6.1.5): true analytical solution plot(t,yt,’k’), hold on klasts = [8 4 2]; hs = tf./klasts; y(1) = y0; for itr = 1:3 %with various step size h = 1/8,1/4,1/2 klast = klasts(itr); h = hs(itr); y(1)=y0; for k = 1:klast y(k + 1) = (1 - a*h)*y(k) +h*r; %Eq.(6.1.3): plot([k - 1 k]*h,[y(k) y(k+1)],’b’, k*h,y(k+1),’ro’) if k < 4, pause; end end end