Ode fichier omer

Code:

clc;clear;
tspan=[0 100];
u0=[20000;100;0];
ode45(@jojo,tspan,u0)


dans une autre fenetre de l'éditeur enregistrer la fonction:

Code:

function yp=jojo(t,u)
yp=[-1/20000*u(1)*u(2);1/20000*u(1)*u(2)-1/10*u(2);1/10*u(2)];
end

Cool, mais on fait quoi si ça ne marche pas?

Est-ce que tu aurais l'exemple pour l'exo 3 de la séance 7?

Moi je le fais en une ligne je comprends rien a tout votre bordel xD
voila pour la séance 7 ex 3 Corentin. (il marche aussi pour l'exemple d'examen en ligne sur son site, suffit de remplacer la fonction. les graphiques se superposent à chaque fois avec Euler et Heun, a peu de chose près)

[x,y]=ode45(f,[0 6],y(1));
plot(x,y,'m')

et f c'est juste la fonction déclarée plus haut
f=@(x,y) 2*y*sin(x)+cos(x)*exp(x/(1+x^2));