Résolution numérique d'équa. diff. avec Matlab

[egan]

Salut tout le monde,

J'ai commencé à prendre en main Matlab et je suis tombé sur quelque chose que je ne parviens pas à expliquer.

J'ai résolu numériquement l'équation différentielle suivante:

y' = 3y - 3t
y(0) = 1/3

avec les fonctions ode23 et ode 45.

On sait que la solution de ce problème est y(t) = t + 1/3.

J'obtiens les résultats suivants donc ça marche plutôt bien.

http://paf.im/Zlkjb

Par contre, ode45 plante avant ode23 alors que ode45 est censé être d'ordre 4 ou 5 et ode23 d'ordre 2 ou 3.
La seule explication rationnelle que j'ai pour justifier ce phénomène est le fait que ces fonctions utilisent des méthodes à pas adaptatifs et encore...

Si vous avez une idée, je suis prenneur. ^^

Merci d'avance.

@+ Boris.

[egan]

Personne ?

[Anonyme]

@egan

Pour info : je viens de télécharger ton fichier pdf
Il n'a été chargé qu'une seule fois , donc ton message ne passionne pas grand monde !

Ton message est peut être un peu trop technique
car visiblement ode23 et ode45 ne semblent pas trop connu par les internautes de Maths-Forum (et moi le premier)

ps)
ce message va servir comme "UP" comme disent certains d'entre nous !