Solveur EDO du Premier Ordre

Solveur EDO du Premier Ordre

Le Solveur EDO du Premier Ordre prend une équation différentielle ordinaire sous la forme dy/dx = f(x, y), classifie automatiquement sa structure (séparable, linéaire, autonome, exacte ou générale), et produit à la fois une solution symbolique sous forme fermée lorsque c'est possible et une solution numérique de haute précision partout. Une visualisation en direct du champ de pentes avec la courbe de solution superposée rend le sens géométrique de l'équation immédiatement évident : les solutions sont exactement les courbes tangentes à chaque flèche.

Qu'est-ce qu'une EDO du premier ordre ?

Une équation différentielle ordinaire du premier ordre implique une fonction inconnue y(x) et sa dérivée première y'(x) uniquement. La forme explicite standard est :

Combinée à une condition initiale y(x₀) = y₀, cela définit un problème de valeur initiale (PVI). Le théorème de Picard-Lindelöf garantit une solution unique dans un certain voisinage de x₀ tant que f est continue de Lipschitz en y près de (x₀, y₀). Géométriquement, le PVI demande la courbe unique passant par (x₀, y₀) dont la pente en chaque point correspond à f en ce point — exactement la courbe tangente au champ de pentes.

Six classes reconnues par le solveur

Méthode de la forme fermée : Linéaire à coefficients constants

Lorsque le membre de droite se simplifie en dy/dx = a·y + b avec des constantes a et b, le facteur intégrant μ(x) = e^(-a·x) donne une solution exacte. La solution générale est :

L'application de la condition initiale y(x₀) = y₀ fixe la constante C et donne l'unique solution particulière. Cette classe unique couvre un nombre énorme de problèmes de manuels scolaires :

• Croissance exponentielle — dy/dx = k·y, solution particulière y(t) = y₀·e^(k·t).
• Décroissance exponentielle — dy/dx = -k·y, demi-vie ln 2 / k.
• Loi de refroidissement de Newton — dy/dx = -k·(y - T_amb), la température du corps se détend de manière exponentielle vers l'ambiante.
• Charge d'un circuit RC — dV/dt = (1/RC)·(V_in - V), la tension du condensateur s'approche de la source.
• Clairance d'un médicament — pharmacocinétique de premier ordre avec taux d'élimination k.

Lecture d'un champ de pentes

À chaque point de la grille (x, y), l'outil trace un court segment de droite dont la pente est égale à f(x, y). Trois observations utiles :

• Les équilibres sont des points où f(x, y) = 0 — le champ de pentes est horizontal. Pour les équations autonomes, ce sont des points fixes y* satisfaisant f(y*) = 0 ; les trajectoires proches s'en approchent (stable) ou s'en éloignent (instable) y*.
• Les isoclines sont des courbes où f(x, y) est égal à une constante c, de sorte que toutes les flèches le long de la courbe ont la même pente c.
• Les courbes de solution ne se croisent jamais (lorsque f est Lipschitz) — visuellement évident car deux courbes qui se croisent auraient besoin de pentes différentes à l'intersection.

Méthode numérique : Runge-Kutta classique (RK4)

Étant donné (x_n, y_n), la valeur suivante est calculée en faisant la moyenne de quatre estimations de pente :

RK4 a une erreur de troncature locale O(h⁵) et une erreur globale O(h⁴), donnant une précision d'environ six chiffres avec le nombre d'étapes par défaut pour les équations non raides. Le solveur intègre vers l'extérieur à partir du point initial dans les deux directions x et s'arrête proprement si la magnitude de y dépasse 10¹⁵ — typique des solutions qui explosent en un temps fini, comme dy/dx = y².

Comment utiliser ce calculateur

1. Entrez le membre de droite dans le champ dy/dx = .... Utilisez x et y comme variables, * pour la multiplication, ^ ou ** pour les puissances, et les fonctions standard comme sin, cos, exp, log, sqrt. Les constantes pi et e sont reconnues.
2. Spécifiez la condition initiale (x₀, y₀) — la courbe de solution unique passera par ce point.
3. Choisissez la plage de x sur laquelle tracer le champ de pentes et la courbe de solution. La plage y est automatiquement ajustée à partir de la solution intégrée.
4. Cliquez sur Résoudre & Visualiser. Le classificateur s'exécute d'abord ; si votre équation correspond à un modèle de forme fermée (linéaire à coefficients constants), vous obtenez la réponse symbolique. Le champ de pentes et la courbe de solution sont toujours rendus.
5. Basculez le champ de pentes pour l'activer ou le désactiver afin de vous concentrer sur la courbe de solution, ou rejouez l'animation du tracé de la courbe pour voir comment l'intégration progresse à partir du point initial.

Exemple concret : Loi de refroidissement de Newton

Une tasse de café à 80 °C refroidit dans une pièce à 20 °C. Le taux de transfert de chaleur est proportionnel à la différence de température :

C'est une équation linéaire à coefficients constants (a = -0.1, b = 2). La forme fermée est :

Après 30 minutes : T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. La vue du champ de pentes rend le comportement limite évident — chaque courbe de solution, quelle que soit la température de départ, tend vers la ligne horizontale T = 20.

Applications courantes

• Dynamique des populations — modèles exponentiels, logistiques, effet Allee.
• Pharmacocinétique — absorption et élimination des médicaments, calculs de demi-vie.
• Transfert de chaleur — loi de refroidissement de Newton, modèles à capacité thermique localisée.
• Circuits RC et RL — transitoires électriques linéaires du premier ordre.
• Désintégration radioactive — chaînes de désintégration d'un seul isotope.
• Réservoirs de mélange — concentration d'un soluté sous flux entrant / sortant.
• Objet en chute avec traînée — analyse de la vitesse terminale dv/dt = g - kv.

Foire Aux Questions

Qu'est-ce qu'une équation différentielle ordinaire du premier ordre ?

Une équation différentielle ordinaire (EDO) du premier ordre est une équation de la forme dy/dx = f(x, y) qui implique la fonction inconnue y(x) et sa dérivée première. Résoudre l'EDO signifie trouver la fonction y(x) dont la dérivée correspond au membre de droite. Avec une condition initiale y(x₀) = y₀, la solution est unique sous des hypothèses de régularité modérées (théorème de Picard-Lindelöf).

Qu'est-ce qu'un champ de pentes ?

Un champ de pentes (ou champ de directions) trace un petit segment de droite à chaque point de grille (x, y) dont la pente est égale à f(x, y). Les courbes de solution de l'EDO sont exactement les courbes qui sont tangentes à ces segments en chaque point. Le champ de pentes donne une intuition visuelle instantanée du comportement global des solutions sans résoudre l'équation symboliquement.

Quelles classes d'EDO du premier ordre cet outil résout-il ?

L'outil classifie automatiquement l'équation parmi les suivantes : intégrable (dépend uniquement de x, résolue par intégration directe), linéaire à coefficients constants y' = a·y + b (forme fermée complète fournie), autonome (dépend uniquement de y), séparable (se factorise en g(x)·h(y)), linéaire à coefficients variables (P(x)·y + Q(x)), ou générale. Pour chaque classe, une solution numérique Runge-Kutta de haute précision et une visualisation du champ de pentes sont produites.

Quelle méthode numérique est utilisée ?

La méthode classique de Runge-Kutta d'ordre 4 (RK4) est appliquée avec 300 sous-étapes dans chaque direction à partir du point initial. RK4 a une erreur de troncature locale O(h⁵) et constitue l'outil standard pour les EDO non raides à cette échelle. Le solveur détecte la divergence (dépassement de capacité ou NaN) et arrête l'intégration proprement afin que le tracé reste valide.

Qu'est-ce que la méthode du facteur intégrant pour les EDO linéaires ?

Pour une EDO linéaire du premier ordre y' + P(x)·y = Q(x), multipliez les deux côtés par le facteur intégrant μ(x) = e^∫P(x) dx. Le côté gauche devient la dérivée exacte d/dx[μ·y], donc y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C). Lorsque P et Q sont des constantes, cela se réduit à la forme fermée y = -b/a + C·e^(a·x), que l'outil renvoie automatiquement.

Cet outil peut-il gérer des équations raides ou des systèmes d'EDO ?

Ce solveur est destiné aux EDO scalaires du premier ordre non raides. Les problèmes très raides (où la solution a plusieurs échelles de temps différant de plusieurs ordres de grandeur) peuvent nécessiter une méthode implicite comme Euler arrière ou Rosenbrock ; les systèmes couplés nécessitent un solveur à valeur vectorielle. Dans ces cas, utilisez un package dédié tel que solve_ivp de SciPy ou des solveurs spécialisés pour les EDO raides.

Lectures complémentaires

• Équation différentielle ordinaire — Wikipédia
• Champ de directions — Wikipédia
• Méthodes de Runge-Kutta — Wikipédia
• Facteur intégrant — Wikipédia
• Théorème de Picard-Lindelöf — Wikipédia