Bonjour tous,

pour des fins pédagogiques j'aimerai reprogrammer la fonction ODE23 de MATLAB, es ce que quelqu'un ici l'a déjà effectué ?

J'aimerai avoir quelques indications, pourriez vous m'aidez s'il vous plait ?

Voici ce que j'ai compris pour le moment et mes interrogations :
  1. il faut choisir un pas initial
  2. il faut faire une résolution avec un Runge Kutta ordre 3
  3. la même résolution avec un Runge Kutta ordre 2
  4. la différence des deux résultats nous donne l'erreur
    ---> mais comment exactement ?
  5. si l'erreur est acceptable on enregistre le résultat et on augment le
    pas de temps suivant
    ---> de combien ?
  6. si l'erreur n'est pas acceptable on accepte pas le pas de temps
    et on refait avec un pas plus petit
    ---> de combien ?
  7. ...Etc

J'espère que vous pourrez me donner des indications car ceci m'intéresse grandement/

merci