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é.
Suites et séries
Suite arithmétique
aₙ = a₁ + (n − 1)d, Sₙ = n(a₁ + aₙ)/2
Suite géométrique
aₙ = a₁rⁿ⁻¹, Sₙ = a₁(1 − rⁿ)/(1 − r)
Série géométrique infinie
S = a₁/(1 − r), |r| < 1
Série de Riemann
Σ1/kᵖ converges only for p > 1
Test du rapport
L = lim |ak+1/ak|: L < 1 converges, L > 1 diverges
Rayon de convergence
R = lim |ck/ck+1|
Polynôme de Taylor
Pₙ(x) = Σ f⁽ᵏ⁾(a)(x − a)ᵏ/k!
Série de Maclaurin
eˣ = Σxᵏ/k!, sin x = x − x³/3! + …, cos x = 1 − x²/2! + …, ln(1 + x) = x − x²/2 + …
Majoration de l’erreur de Lagrange
|Rₙ| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!
Méthode
Identifiez le type : arithmétique (différence constante), géométrique (rapport constant), une série usuelle, une série entière ou un polynôme de Taylor.
Pour les sommes, utilisez Sₙ = n(a₁ + aₙ)/2, Sₙ = a₁(1 − rⁿ)/(1 − r), ou S = a₁/(1 − r) quand |r| < 1, avec a₁ le premier terme de la somme.
Pour la convergence : d’abord le test du terme général, puis comparez avec Σ1/kᵖ ou une série géométrique, et utilisez le test du rapport ou de la racine pour les factorielles et les puissances k-ièmes.
Pour les polynômes de Taylor, listez f(a), f′(a), f″(a), …, divisez la dérivée k-ième par k!, et majore l’erreur par M·|x − a|ⁿ⁺¹/(n + 1)!.
Matrices et déterminants
déterminant 2×2
det = ad − bc
Inverse d’une matrice 2×2
A⁻¹ = (1/det)·[d −b; −c a]
Produit
(AB)ij = Σk AikBkj
Valeurs propres
det(A − λI) = 0
Règle de Cramer
xi = det(Ai)/det(A)
Équations différentielles
Croissance/décroissance exponentielle
dy/dt = ky → y = y0ekt
Demi-vie
k = ln 2 / t½
Loi de refroidissement de Newton
T(t) = Ta + (T0 − Ta)e−kt
Facteur intégrant
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Cuve de mélange
dm/dt = r·c − (r/V)·m
2Exercices
1.Méthodes numériquesMoyen
Approchez ∫01,5 e−x² dx par la méthode de Simpson avec n = 4 bandes. Donnez S.
f(x) = e−x²[a, b] = [0, 1,5]n = 4
S =
2.Méthodes numériquesDifficile
Utilisez la méthode d’Euler avec h = 0,2 pour y′ = e−ty, y(0) = 1 afin d’estimer y4 ≈ y(0,8).
y′ = e−tyy(0) = 1h = 0,2
y4 =
3.Suites et sériesMoyen
Trouvez la somme S = Σk=1∞ 3/(k(k + 4)).
ak = 3/(k(k + 4))
S =
4.Suites et sériesDifficile
f(x) = cos x est approchée par son polynôme de Maclaurin de degré 3 en x = 1. Utilisez la majoration de l’erreur de Lagrange M·|x|ⁿ⁺¹/(n + 1)! avec M = 1 pour trouver la plus grande erreur possible |R₃|.
f(x) = cos xn = 3x = 1M = 1
|R₃| =
5.Matrices et déterminantsMoyen
Calcule le déterminant det A de la matrice 3×3 A représentée.
A = [−4 5 4; 2 1 −4; −3 4 −1]
det A =
6.Matrices et déterminantsDifficile
Trouve la plus grande valeur propre λ de la matrice A représentée.
A = [1 −3; 5 −7]
λ =
7.Équations différentiellesMoyen
Résolvez y″ + 4y′ + 4y = 0 avec y(0) = 2 et y′(0) = -2. Trouvez y(2).
y″ + 4y′ + 4y = 0y(0) = 2y′(0) = -2
y(2) =
8.Équations différentiellesDifficile
Un système obéit à y″ + y′ + 25,25y = 126,3 pour t ≥ 0, à partir du repos : y(0) = 0, y′(0) = 0. Trouvez y(0,25) ; angles en radians.
y″ + y′ + 25,25y = 126,3y(0) = 0y′(0) = 0
y(0.25) =
3Solutions
1.Méthodes numériques
Réponse: S = 0,8561
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.
2.Méthodes numériques
Réponse: y4 = 1,556
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⁴.
Pas d’Euler
y₁ = 1 + 0,2·(1) = 1,2
y₂ = 1,2 + 0,2·(0,786628) = 1,357326
y₃ = 1,357326 + 0,2·(0,581043) = 1,473534
y₄ = 1,473534 + 0,2·(0,413077) = 1,55615
Graphique
y4 = 1,556
Réponse
y4 = 1,556
3.Suites et séries
Réponse: S = 1,562
Convergence : si les termes ne tendent pas vers 0, la série diverge (test du terme général). Σ1/kᵖ ne converge que pour p > 1, et Σrᵏ que pour |r| < 1. Comparez une série à l’une d’elles quand ses termes se comportent comme celles-ci pour k grand. Test du rapport : L = lim |ak+1/ak| ; test de la racine : L = lim |ak|1/k. L < 1 converge, L > 1 diverge, L = 1 ne décide rien. Une série entière Σck(x − x₀)ᵏ converge pour |x − x₀| < R avec R = lim |ck/ck+1| ; une série alternée dont les termes décroissent vers 0 converge.
Décomposition en éléments simples
3/(k(k + 4)) = (3/4)·(1/k − 1/(k + 4))
Annuler
S = (3/4)·(1 + 1/2 + 1/3 + 1/4) = 1,562
Réponse
S = 1,562
4.Suites et séries
Réponse: |R₃| = 0,04167
Polynômes de Taylor : près de x = a, une fonction régulière est approchée par Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n! ; pour a = 0, c’est le polynôme de Maclaurin. Séries usuelles : eˣ = Σxᵏ/k!, sin x = x − x³/3! + x⁵/5! − …, cos x = 1 − x²/2! + x⁴/4! − …, ln(1 + x) = x − x²/2 + x³/3 − …, 1/(1 − x) = Σxᵏ. Majoration de l’erreur de Lagrange : |f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!, où M est la plus grande valeur de |f⁽ⁿ⁺¹⁾| entre a et x.
Majoration de la dérivée
M = 1
Majoration de Lagrange
|R₃| ≤ 1·1⁴/4! = 0,04167
Réponse
|R₃| = 0,04167
L’erreur réelle est plus petite : |f(1) − P₃(1)| = 0,0403.
5.Matrices et déterminants
Réponse: det A = 54
Matrices : le déterminant d'une matrice 2×2 est ad − bc. Pour une 3×3, développe selon une ligne : chaque coefficient multiplié par le déterminant de la matrice 2×2 qui reste après avoir rayé sa ligne et sa colonne, avec des signes alternés + − +. Les valeurs propres λ vérifient det(A − λI) = 0 ; en 2×2 c'est λ² − (trace)·λ + det = 0.
Mineurs 2×2 de la première ligne
M₁₁ = 15, M₁₂ = −14, M₁₃ = 11
Développement par cofacteurs
det A = (−4)·15 − 5·(−14) + 4·11 = 54
Réponse
det A = 54
6.Matrices et déterminants
Réponse: λ = −2
Matrices : le déterminant d'une matrice 2×2 est ad − bc. Pour une 3×3, développe selon une ligne : chaque coefficient multiplié par le déterminant de la matrice 2×2 qui reste après avoir rayé sa ligne et sa colonne, avec des signes alternés + − +. Les valeurs propres λ vérifient det(A − λI) = 0 ; en 2×2 c'est λ² − (trace)·λ + det = 0.
Équation caractéristique
λ² − (−6)λ + (8) = 0
Résoudre
λ = (−6 ± √(4))/2 → λ = −4, −2
λ = −2
Réponse
λ = −2
7.Équations différentielles
Réponse: y(2) = 0,1099
Équations linéaires du second ordre : pour y″ + b·y′ + c·y = 0, essayez y = ert : on obtient l’équation caractéristique r² + b·r + c = 0. Deux racines réelles r₁ ≠ r₂ donnent y = A·er₁t + B·er₂t (suramorti) ; une racine double r donne y = (A + Bt)·ert (amortissement critique) ; des racines complexes −α ± jω donnent y = e−αt(A cos ωt + B sin ωt) (sous-amorti, oscillant). Les valeurs initiales y(0) et y′(0) fixent A et B. Écrite sous la forme y″ + 2ζωny′ + ωn²y = 0, ωn est la pulsation propre et ζ le coefficient d’amortissement : ζ < 1 donne des oscillations. Avec une entrée constante, y″ + by′ + cy = c·K, la solution est la valeur stationnaire K plus la réponse libre.
Racine double
(r + 2)² = 0 → r = −2
Constantes
A = y(0) = 2
y′(0) = B − 2A → B = -2 + 4 = 2
Évaluer
y(2) = (2 + 2·2)·e−4 = 0,1099
Réponse
y(2) = 0,1099
8.Équations différentielles
Réponse: y(0.25) = 3,19
Équations linéaires du second ordre : pour y″ + b·y′ + c·y = 0, essayez y = ert : on obtient l’équation caractéristique r² + b·r + c = 0. Deux racines réelles r₁ ≠ r₂ donnent y = A·er₁t + B·er₂t (suramorti) ; une racine double r donne y = (A + Bt)·ert (amortissement critique) ; des racines complexes −α ± jω donnent y = e−αt(A cos ωt + B sin ωt) (sous-amorti, oscillant). Les valeurs initiales y(0) et y′(0) fixent A et B. Écrite sous la forme y″ + 2ζωny′ + ωn²y = 0, ωn est la pulsation propre et ζ le coefficient d’amortissement : ζ < 1 donne des oscillations. Avec une entrée constante, y″ + by′ + cy = c·K, la solution est la valeur stationnaire K plus la réponse libre.