Exemple : Résolution d'un système d'EDO de premier ordre
Utilisez un bloc de résolution et la fonction odesolve pour résoudre des équations différentielles ordinaires de premier ordre.
1. Définissez le point d'extrémité de l'intervalle de la solution.
2. Définissez le problème à l'aide des dérivées et d'un ensemble de conditions initiales :
3. Tracez les solutions sur un seul intervalle :
Utilisation de la fonction Rkadapt
Utilisez la fonction Rkadapt pour résoudre le même système d'équations différentielles :
1. Définissez une fonction qui détermine un vecteur de valeurs dérivées à un point de la solution (t,Y) :
2. Définissez des arguments supplémentaires pour le solveur d'EDO :
a. Valeur initiale d'une variable indépendante
b. Vecteur des valeurs de fonction initiales
c. Nombre de valeurs de la solution sur [t0, t1]
3. Utilisez la fonction Rkadapt pour trouver la matrice Solution :
| Vous pouvez également utiliser les fonctions rkfixed, Bulstoer ou Radau. |
4. Extrayez les valeurs de la variable indépendante :
5. Extrayez les valeurs de la première fonction de la solution :
6. Extrayez les valeurs de la deuxième fonction de la solution :
7. Extrayez les valeurs de la troisième fonction de la solution :
8. Tracez les trois solutions :
| Le tracé de la solution odesolve est quasiment identique à celui de la solution Rkadapt. |
