% on efface tout !
clear all
% definition de la fonction "1"
un=@(t,y)[1];
% definition de la fonction f(t,y)
f=@(t,y)[t.*y];
%definition du champ de vitesses
norm=@(t,y)[sqrt(1+f(t,y).**2)]
xx=@(t,y)[1./norm(t,y)];
ff=@(t,y)[f(t,y)./norm(t,y)];
% on definit la grille pour le champ de vitesses
[t,y]=meshgrid(-3:0.2:3,-2:0.2:2);
% on trace le champ de vitesses
quiver(t,y,xx(t,y),ff(t,y),0.8); %le 0.4 fixe la taille des fleches
% on tracera un autre graphique donc on conserve le plot
hold on;
% definition de la condition initiale ya
ya=3;
% definition de l'intervalle sur lequel on resoud
a=-2;
b=2;
% on cree les vecteurs des abscisses et ordonnees pour la solution exacte
g=@(y,t)[f(t,y)];
%sol=@(t)[exp(t)];
xxx=linspace(a,b,100);
sol=lsode(g,ya,xxx);
% on trace cette solution sur le meme graphique
plot(xxx,sol,'r-',a,ya,'ko');
hold off;