Pour une racine : vérifiez le changement de signe pour la dichotomie, ou calculez f et f′ pour Newton, et gardez six chiffres ou plus dans un tableau d’itérés.
Pour une valeur de tableau : interpolez linéairement entre les deux points voisins.
Pour une intégrale : trouvez h = (b − a)/n, listez f en chaque nœud, appliquez les poids (1, 2, …, 2, 1 pour les trapèzes ; 1, 4, 2, …, 4, 1 pour Simpson avec n pair) et multipliez par h/2 ou h/3.
Pour une équation différentielle : avance depuis t₀ avec yₙ₊₁ = yₙ + h·f(tₙ, yₙ) (Euler) ou la pente moyenne (Heun), une ligne par pas, jusqu’à atteindre le t demandé.
2Exercices
1.Équations différentiellesMoyen
Une tasse de café à 94 °C est laissée dans une pièce à 21 °C. Au bout de 14 min, sa température est de 63 °C. Elle suit la loi de refroidissement de Newton, dT/dt = −k·(T − Tₐ). Combien de temps après le début, t, lui faut-il pour atteindre 85 °C ?
Soit f(x) = x·sin(x), avec x en radians. Calcule f′(1,5).
f(x) = x·sin(x)x = 1,5 rad
f′(1.5) =
6.DérivéesDifficile
La tangente à f(x) = 2x2 + 2x + 2 en x = 2 coupe l'axe des y en y = b. Trouve b.
f(x) = 2x2 + 2x + 2x = 2
b =
7.Méthodes numériquesMoyen
Approchez ∫02 1/√(1 + x⁴) dx par la méthode de Simpson avec n = 6 bandes. Donnez S.
f(x) = 1/√(1 + x⁴)[a, b] = [0, 2]n = 6
S =
8.Méthodes numériquesDifficile
Faites un pas de la méthode de Runge–Kutta classique (RK4) avec h = 0,2 pour y′ = cos t − y, y(0) = 2. Trouvez y₁ ≈ y(0,2).
y′ = cos t − yy(0) = 2h = 0,2
y1 =
3Solutions
1.Équations différentielles
Réponse: t = 3,332 min
Équations différentielles du premier ordre : dy/dt = k·y a pour solution y = y₀·ekt (croissance pour k > 0, décroissance pour k < 0) ; une demi-vie T donne k = ln 2/T. La loi de refroidissement de Newton, dT/dt = −k·(T − Tₐ), donne T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Une équation linéaire y′ + p(x)·y = q(x) se résout avec le facteur intégrant μ = e∫p dx : (μy)′ = μq, donc μy = ∫μq dx + C.
Constante de refroidissement
k = ln(73/42)/14 = 0,03948 min⁻¹
Temps pour atteindre la valeur visée
t = ln(73/64)/0,03948 = 3,332 min
Réponse
t = 3,332 min
2.Équations différentielles
Réponse: y(2) = 1,472
Équations différentielles du premier ordre : dy/dt = k·y a pour solution y = y₀·ekt (croissance pour k > 0, décroissance pour k < 0) ; une demi-vie T donne k = ln 2/T. La loi de refroidissement de Newton, dT/dt = −k·(T − Tₐ), donne T(t) = Tₐ + (T₀ − Tₐ)·e−kt. Une équation linéaire y′ + p(x)·y = q(x) se résout avec le facteur intégrant μ = e∫p dx : (μy)′ = μq, donc μy = ∫μq dx + C.
Facteur intégrant
μ = e3t
(e3ty)′ = 10e2,5t
Solution générale
y = 10/2,5·e−0,5t + C·e−3t = 4·e−0,5t + C·e−3t
Condition initiale
4 = 4 + C → C = 0
Évaluer
y(2) = 4·0,3679 + 0·0,002479 = 1,472
Réponse
y(2) = 1,472
3.Intégrales
Réponse: I = 3,219
Intégrales définies : ∫ₐᵇ f(x) dx est l'aire algébrique entre le graphe et l'axe des x de a à b. Trouve une primitive F (avec F′ = f) et calcule F(b) − F(a). Règle des puissances : ∫xⁿ dx = xⁿ⁺¹/(n + 1). L'aire entre deux courbes est ∫(haute − basse) dx.
Primitive
F(x) = 2·ln x
Évaluer aux bornes
I = 2·(ln 10 − ln 2) = 2·ln 5 = 3,219
Réponse
I = 3,219
4.Intégrales
Réponse: I = 3,602
Intégration par parties : ∫u·dv = u·v − ∫v·du, issue de la règle du produit. Choisissez pour u le facteur qui se simplifie quand on le dérive (une puissance de x, ou ln x), et pour dv le facteur que vous savez intégrer (eˣ, sin x, cos x, ou x lorsque u = ln x). Pour une intégrale définie, évaluez u·v aux deux bornes : ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. Un facteur x² demande deux étapes.
Première étape
∫x²eˣ dx = x²eˣ − ∫2x eˣ dx
Deuxième étape
∫2x eˣ dx = 2x eˣ − 2eˣ
∫x²eˣ dx = eˣ(x² − 2x + 2)
Limites
4,482·(1,25) − 1·2 = 3,602
Réponse
I = 3,602
5.Dérivées
Réponse: f′(1.5) = 1,104
Dérivation : la dérivée f′(x) est la pente du graphe en x. Règle des puissances : (xⁿ)′ = n·xⁿ⁻¹. Règle du produit : (uv)′ = u′v + uv′. Règle de dérivation en chaîne : (f(g(x)))′ = f′(g(x))·g′(x). Là où f′(x) = 0, le graphe a une tangente horizontale : un maximum, un minimum ou un point selle.
Règle du produit et en chaîne
f′(x) = sin(x) + x·cos(x)
Remplacer
f′(1,5) = 0,9975 + 1,5·0,07074 = 1,104
Réponse
f′(1,5) = 1,104
6.Dérivées
Réponse: b = −6
Dérivation : la dérivée f′(x) est la pente du graphe en x. Règle des puissances : (xⁿ)′ = n·xⁿ⁻¹. Règle du produit : (uv)′ = u′v + uv′. Règle de dérivation en chaîne : (f(g(x)))′ = f′(g(x))·g′(x). Là où f′(x) = 0, le graphe a une tangente horizontale : un maximum, un minimum ou un point selle.
Point et pente
f(2) = 14
f′(x) = 4x + 2 → f′(2) = 10
Tangente en x = 0
b = f(x₀) − f′(x₀)·x₀ = 14 − 10·2 = −6
Réponse
b = −6
7.Méthodes numériques
Réponse: S = 1,357
Intégration numérique : découpe [a, b] en n bandes de largeur h = (b − a)/n avec les nœuds x₀ … xₙ. Méthode des trapèzes : T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. Méthode de Simpson (n pair) : S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. L’erreur des trapèzes décroît comme h² (majoration (b − a)h²·M/12 avec M = max|f″|), celle de Simpson comme h⁴, donc Simpson est en général bien plus précise pour les mêmes nœuds.
Cette intégrale n’a pas de primitive élémentaire, ce qui est exactement le cas où une méthode numérique est nécessaire.
8.Méthodes numériques
Réponse: y1 = 1,817
Résolution numérique de y′ = f(t, y) : la méthode d’Euler avance selon la pente au début de chaque pas : yn+1 = yₙ + h·f(tₙ, yₙ). La méthode de Heun (Euler amélioré) moyenne la pente au début et à une prédiction d’Euler : k₁ = f(tₙ, yₙ), k₂ = f(tₙ + h, yₙ + h·k₁), yn+1 = yₙ + (h/2)(k₁ + k₂). La méthode de Runge–Kutta classique (RK4) combine quatre pentes : yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). L’erreur globale d’Euler est proportionnelle à h, celle de Heun à h², celle de RK4 à h⁴.