% 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 resoud pour la valeur n=2^k, faire les differentes valeurs avec k=2,...,10
k=10;
n=2**k;
[x u]=euler_explicite1(f,ya,a,b,n);
% on cree les vecteurs des abscisses et ordonnees pour la solution exacte
sol=@(t)[exp(t)];
xx=linspace(a,b,100);
v=sol(xx);
% on trace les deux fonctions sur un meme dessin
plot(x,u,'b-',xx,v,'r-');
legend(["u";"sol"]);