ode23 à programmer

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

Je pense qu'il est possible de trouver le fichier source de cette fonction sur le site d'octave ou sur le site de scilab...
J'ai jeté un coup d'oeil sur le site d'octave et je n'ai pas trouvée...
J'ai surement fait une mauvaise recherche car je pense que ceci doit y etre

La référence à l'algorithme utilisé par MATLAB est donnée dans la documentation (cliquer sur Algorithms)

merci DUT je crois que j'ai tous les éléments à présent