% on efface tout !
clear all
% definition de la fonction f(t,y)
f=@(t,y)[y];
% definition de la condition initiale ya
ya=1;
% definition de l'intervalle sur lequel on resoud
a=0;
b=2;
% on definit la solution exacte du probleme de Cauchy
sol=@(t)[exp(t)];
% on resoud pour les valeurs n=2^k, avec k=2,...,10,
% et on stocke la difference |u(2)-sol(2)| dans le vecteur erreur
erreur=[];
for k=2:13
  n=2**k;
  [x u]=euler_explicite1(f,ya,a,b,n);
  erreur = [ erreur abs(sol(2)-u(length(u)))];
endfor
% on trace l'erreur en fonction de n sur une échelle logarithmique
k=2:13;
err=log10(erreur);
plot(k,erreur,'b-');