Ode45

Coucou tout le monde!

Est-ce que quelqu'un aurait compris comment on fait pour utiliser la fonction ode45 pour l'exercice 1 de la séance 5.
Ce que j'essaye n'a pas l'air de fonctionner...

Apparemment les autres années, ils avaient une équation avec des x et des y, est-ce que nous, c'est le t et le u?
Mais alors on doit intégrer l'équation (4)?

Bref, c'est un peu confus...

Merci à ceux qui me répondront!

apparemment dans les corrections qu'il y donné il y a que le u de l'équa qui varie

bon j'estime que ma torture de cerveau a assez duré :p
si qqun pouvait m'aider ou me corriger --> jcrois que c'est la super fonction qui pause des bugs (on a qu'une variable au lieu de 2 comme il faudrait...)

j'ai 2 exemples pour ODE45

Code:

clear;clc;clf;
a=0;
b=200;
N=50;
h=(b-a)/N;
y(1)=-1;
f=@(x,y) -5/y.^2+sin(x/5);

[t u]=ode45(f,[a b],y(1),[]);
plot(t,u,'r')
grid on

celui la il marche parfaitement mais la fonction c'est pas celle de l'énoncé --> c'est une fonction habituelle à 2 variables (comme dans la théorie...)

et l'autre le graphe bug mais ya plus d'erreur (ouf)

Code:

clear;clc;clf;
a=0;
b=200;
N=50;
h=(b-a)/N;
y(1)=-1;
rho=3.3;
q=10;
t=a:h:b;
f=@(t,u) rho*u.*(1-u/q)-((u.^2)/(1+u.^2));

[t u]=ode45(f,[a b],y(1),[]);
plot(t,u,'r')
grid on

qqun a une idée pour qu'on ai une belle courbe comme d'hab au 2e??????????????
HELP : s

Question:
ode45 sert à vérifier la méthode qu'on a utilisé(comme quad et fzero)?
Donc il faut moduler nos conditions initiales pour que les courbes se superposent? ou rien à voir?