UNIVERSITE MAPON

KINDU

FACULTE POLYTECHNIQUE

COURS DE MATHEMATIQUES APPLIQUEES

RESOLUTIONS DES EXERCICES DES TRAVAUX DIRIGES

TRAVAIL FAIT PAR :

1. BAHATI KADUKU Cédrick
2. JONATHAN BYANDI Emmanuel
3. LUMUMBA MUSHIKO Martyr
4. MALIANI KAMBULU Adolphe
5. MATALA MAKIESE
6. MIRINDI KWABWINE Ezéchias
7. MOBI KITAMBALA Joseph
8. MURHULA NGUVUMALI Christian
9. OMARI BIN SEFU Prudent
10. ONEMA PENGE Josué

PROMOTIONS : Master 1 Electro-Energétique et Master 1 Génie Informatique

ENCADREURS :
Professeur Joseph BUKWELI
Assistant MUSA RACHIDI Thomas

MARS 2023
 1

QUESTION I. Equations différentielles à coefficients constants.

1. Résoudre les équations ci-après en utilisant la méthode appropriée :

a) 𝟐𝒚′′ − 𝒚′ + 𝒚 = (𝒕𝟐 + 𝟏)𝒆𝒕 (1)
Solution
Equation homogène : 2𝑦 ′′ − 𝑦 ′ + 𝑦 = 0
Équation caractéristique : 2𝑟 2 − 𝑟 + 1 = 0

∆ = 1 − 8 = −7 √∆= 𝑖 √7
1+𝑖√7
𝑟1 =
4
{
1−𝑖√7
𝑟2 =
4
𝑡 𝑡
√7 √7
𝑆𝐹𝑆 = {𝑦1 = 𝑒 4 cos ( 𝑡) ; 𝑦2 = 𝑒 4 sin ( 𝑡)}
4 4

Solution homogène :
𝑡 √7 𝑡 √7
𝑦ℎ = 𝐶1 𝑦1 + 𝐶2 𝑦2 = 𝐶1 𝑒 4 cos ( 𝑡) + 𝐶2 𝑒 4 sin ( 𝑡)
4 4
Cherchons la solution particulière par la méthode d’identification ; elle est de la forme :
𝑦𝑝 = (𝑎𝑡 2 + 𝑏𝑡 + 𝑐 )𝑒 𝑡

𝑦𝑝′ = (2𝑎𝑡 + 𝑏)𝑒 𝑡 + (𝑎𝑡 2 + 𝑏𝑡 + 𝑐 )𝑒 𝑡 = (𝑎𝑡 2 + (2𝑎 + 𝑏)𝑡 + 𝑏 + 𝑐 )𝑒 𝑡

𝑦𝑝′′ = (𝑎𝑡 2 + (4𝑎 + 𝑏)𝑡 + 2𝑎 + 2𝑏 + 𝑐 )𝑒 𝑡

En remplaçant dans (1) et en divisant membre à membre par 𝑒 𝑡 on a :

2(𝑎𝑡 2 + (4𝑎 + 𝑏)𝑡 + 2𝑎 + 2𝑏 + 𝑐 ) − (𝑎𝑡 2 + (2𝑎 + 𝑏)𝑡 + 𝑏 + 𝑐 ) + (𝑎𝑡 2 + 𝑏𝑡 + 𝑐 ) = (𝑡 2 + 1)
⇒ (2𝑎 − 𝑎 + 𝑎)𝑡 2 + (8𝑎 + 2𝑏 − 2𝑎 − 𝑏 + 𝑏)𝑡 + (4𝑎 + 4𝑏 + 2𝑐 − 𝑏 − 𝑐 + 𝑐 ) = 𝑡 2 + 1
⇒ 2𝑎𝑡 2 + (6𝑎 + 2𝑏)𝑡 + (4𝑎 + 3𝑏 + 2𝑐 ) = 𝑡 2 + 1

Par identification on a :
2𝑎 = 1
{ 6𝑎 + 2𝑏 = 0
4𝑎 + 3𝑏 + 2𝑐 = 1
1 3 7
En résolvant ce système on trouve : 𝑎 = , 𝑏= − et 𝑐= . D’où la solution particulière :
2 2 4

1 3 7
𝑦𝑝 = ( 𝑡 2 − 𝑡 + ) 𝑒 𝑡
2 2 4
Ainsi la solution générale à l’EDO (1) est : 𝑦 = 𝑦ℎ + 𝑦𝑝 :

𝑡 √7 𝑡 √7 1 3 7
𝑦 = 𝐶1 𝑒 4 cos ( 𝑡) + 𝐶2 𝑒 4 sin ( 𝑡) + ( 𝑡 2 − 𝑡 + ) 𝑒 𝑡
4 4 2 2 4
 2

b) 𝟑𝒚′′ + 𝟒𝒚′ + 𝒚 = (𝐬𝐢𝐧 𝒕)𝒆−𝒕 (1)

Solution
Equation homogène : 3𝑦 ′′ + 4𝑦 ′ + 𝑦 = 0
Équation caractéristique : 3𝑟 2 + 4𝑟 + 1 = 0

∆ = 16 − 12 = 4 √∆= 2
−4+2
𝑟1 = = −1
6
{ −4−2 1
𝑟2 = =−
6 3
1
𝑆𝐹𝑆 = {𝑦1 = 𝑒 −𝑡 , 𝑦2 = 𝑒 −3𝑡 }

Solution homogène :
1
𝑦ℎ = 𝐶1 𝑦1 + 𝐶2 𝑦2 = 𝐶1 𝑒 −𝑡 + 𝐶2 𝑒 −3𝑡
Cherchons la solution particulière par la méthode d’identification. Elle est de la forme :
𝑦𝑝 = (𝑎 cos 𝑡 + 𝑏 sin 𝑡)𝑒 −𝑡

𝑦𝑝′ = (𝑏 cos 𝑡 − 𝑎 sin 𝑡)𝑒 −𝑡 − (𝑎 cos 𝑡 + 𝑏 sin 𝑡)𝑒 −𝑡 = ((𝑏 − 𝑎) cos 𝑡 − (𝑎 + 𝑏) sin 𝑡)𝑒 −𝑡

𝑦𝑝′′ = −2(𝑏 cos 𝑡 − 𝑎 sin 𝑡)𝑒 −𝑡

En remplaçant dans (1) et en divisant membre à membre par 𝑒 −𝑡 on a :

3(−2(𝑏 cos 𝑡 − 𝑎 sin 𝑡)) + 4((𝑏 − 𝑎) cos 𝑡 − (𝑎 + 𝑏) sin 𝑡) + (𝑎 cos 𝑡 + 𝑏 sin 𝑡) = sin 𝑡
⇒ (6𝑎 − 4𝑎 − 4𝑏 + 𝑏) sin 𝑡 + (−6𝑏 + 4𝑏 − 4𝑎 + 𝑎) cos 𝑡 = sin 𝑡
Par identification on a :
2𝑎 − 3𝑏 = 1
{
−2𝑏 − 3𝑎 = 0
2 3
En résolvant ce système on trouve : 𝑎 = , 𝑏= − . D’où la solution particulière :
13 13
1
𝑦𝑝 = (2 cos 𝑡 − 3 sin 𝑡)𝑒 −𝑡
13

Ainsi la solution générale est : 𝑦𝑔 = 𝑦ℎ + 𝑦𝑝

1 1
𝑦 = 𝐶1 𝑒 −𝑡 + 𝐶2 𝑒 −3𝑡 + (2 cos 𝑡 − 3 sin 𝑡)𝑒 −𝑡
13
 3

2. Trouver la solution générale de l’équation différentielle

𝑡 𝑑𝑦 𝑡
𝑒 𝑦 (𝑦 − 𝑡) ( ) − 𝑦 (1 + 𝑒 𝑦 ) = 0
𝑑𝑡
Solution
𝑡 𝑑𝑡 𝑑𝑢 𝑑𝑦 𝑑𝑦
Posons = 𝑢 ⇒ 𝑡 = 𝑢𝑦 ⇒ 𝑑𝑡 = 𝑢𝑑𝑦 + 𝑦𝑑𝑢 ⇒ =𝑢+𝑦 ⇒ =
𝑦 𝑑𝑦 𝑑𝑦 𝑑𝑡 𝑢𝑑𝑦+𝑦𝑑𝑢
L’équation devient donc :
𝑑𝑦
𝑒 𝑢 (𝑦 − 𝑢𝑦) − 𝑦(1 + 𝑒 𝑢 ) = 0 ⇒ 𝑦𝑒 𝑢 (1 − 𝑢)𝑑𝑦 − 𝑦(1 + 𝑒 𝑢 )(𝑢𝑑𝑦 + 𝑦𝑑𝑢 ) = 0
𝑢𝑑𝑦 + 𝑦𝑑𝑢
⇒ 𝑒 𝑢 (1 − 𝑢)𝑑𝑦 − (𝑢𝑑𝑦 + 𝑦𝑑𝑢 + 𝑢𝑒 𝑢 𝑑𝑦 + 𝑦𝑒 𝑢 𝑑𝑢) = 0
⇒ 𝑒 𝑢 𝑑𝑦 − 𝑢𝑒 𝑢 𝑑𝑦 − 𝑢𝑑𝑦 − 𝑦𝑑𝑢 − 𝑢𝑒 𝑢 𝑑𝑦 − 𝑦𝑒 𝑢 𝑑𝑢 = 0
⇒ (𝑒 𝑢 − 2𝑢𝑒 𝑢 − 𝑢)𝑑𝑦 − 𝑦(1 + 𝑒 𝑢 )𝑑𝑢 = 0

Equation différentielle à variables séparables :

1 + 𝑒𝑢 𝑑𝑦 1 + 𝑒𝑢
⇒ 𝑦(1 + 𝑒 𝑢 )𝑑𝑢 = (𝑒 𝑢 − 2𝑢𝑒 𝑢 − 𝑢)𝑑𝑦 ⇒ 𝑑𝑢 = ⇒ 𝑑𝑢
𝑒 𝑢 − 2𝑢𝑒 𝑢 − 𝑢 𝑦 (2𝑢 − 1)𝑒 𝑢 + 𝑢
𝑑𝑦
=−
𝑦

En intégrant membre à membre on a :

1 + 𝑒𝑢 𝑑𝑦
∫ 𝑑𝑢 = ∫ −
(2𝑢 − 1)𝑒 𝑢 + 𝑢 𝑦
𝑢
1+𝑒
𝐼1 = ∫ 𝑑𝑢
(2𝑢 − 1)𝑒 𝑢 + 𝑢
𝑑𝑦 𝑐
𝐼2 = ∫ − = − ln 𝑐𝑦 ⇒ 𝐼2 = ln ⇒ 𝑦 = 𝐶1 𝑒 −𝐼1
𝑦 𝑦

1+𝑒 𝑢 𝑡
−∫ 𝑑𝑢
⇒ 𝑦 = 𝐶1 𝑒 (2𝑢−1)𝑒 𝑢 +𝑢 , 𝑢=
𝑦

3. Résoudre le problème aux valeurs initiales

𝑑𝑡
𝑡 ( ) = 𝑦 + √(𝑡 2 + 𝑦 2 ), 𝑦(1) = 0
𝑑𝑦
Solution
Ce problème peut se mettre sous la forme :

𝑑𝑦 𝑡2 𝑑𝑦 𝑡2 𝑡 𝑑𝑦 𝑡2
𝑡 = 𝑦 + 𝑦 √1 + 2 ⇒ 𝑡 = 𝑦 (1 + √1 + 2 ) ⇒ = (1 + √1 + 2 )
𝑑𝑡 𝑦 𝑑𝑡 𝑦 𝑦 𝑑𝑡 𝑦

𝑡 𝑑𝑦 𝑑𝑦
Posons = 𝑢 ⇒ 𝑡 = 𝑢𝑦 ⇒ 𝑑𝑡 = 𝑢𝑑𝑦 + 𝑦𝑑𝑡 ⇒ = . L’équation devient :
𝑦 𝑑𝑡 𝑢𝑑𝑦+𝑦𝑑𝑢
 4

𝑑𝑦
𝑢( ) = 1 + √1 + 𝑢2 ⇒ 𝑢𝑑𝑦 = 𝑢𝑑𝑦 + 𝑦𝑑𝑢 + (𝑢𝑑𝑦 + 𝑦𝑑𝑢)√1 + 𝑢2
𝑢𝑑𝑦 + 𝑦𝑑𝑢
⇒ 𝑦𝑑𝑢 + (𝑢𝑑𝑦 + 𝑦𝑑𝑢)√1 + 𝑢2 = 0 ⇒ 𝑦 (1 + √1 + 𝑢2 ) 𝑑𝑢 + 𝑢√1 + 𝑢2 𝑑𝑦 = 0
Qui est une EDO à variable séparables :
1 + √1 + 𝑢2 𝑑𝑦
𝑦 (1 + √1 + 𝑢2 ) 𝑑𝑢 = −𝑢√1 + 𝑢2 𝑑𝑦 ⇒ 𝑑𝑢 = −
𝑢√1 + 𝑢2 𝑦

1 + √1 + 𝑢2 𝑑𝑦
∫ 𝑑𝑢 = − ∫
𝑢√1 + 𝑢2 𝑦
Premier membre :

1 + √1 + 𝑢2 𝑑𝑢 √1 + 𝑢2 𝑑𝑢 𝑑𝑢
∫ 𝑑𝑢 = ∫ +∫ 𝑑𝑢 = ∫ +∫ = 𝐼1 + 𝐼2
𝑢√1 + 𝑢2 𝑢√1 + 𝑢2 𝑢√1 + 𝑢2 𝑢√1 + 𝑢2 𝑢

𝑑𝑢
𝐼1 = ∫
𝑢√1 + 𝑢2
Posons 𝑢 = tan 𝑣 ⇒ 𝑑𝑢 = sec 2 𝑣 𝑑𝑣
𝑑𝑢 sec 2 𝑣 sec 2 𝑣 sec 𝑣 𝑑𝑣
𝐼1 = ∫ =∫ 𝑑𝑣 = ∫ 𝑑𝑣 = ∫ 𝑑𝑣 = ∫
𝑢√1 + 𝑢2 tan 𝑣 √1 + tan2 𝑣 tan 𝑣 sec 𝑣 tan 𝑣 sin 𝑣
= ln|cosec 𝑣 − cot 𝑣|

1 sec 𝑣 √1+𝑢2 1 1
Ayant posé 𝑢 = tan 𝑣 ⇒ √1 + 𝑢2 = sec 𝑣, on a : cosec 𝑣 = = = et cot 𝑣 = = ;
sin 𝑣 tan 𝑣 𝑢 tan 𝑣 𝑢
alors :
√1 + 𝑢2 − 1
𝐼1 = ln|cosec 𝑣 − cot 𝑣| = ln | |
𝑢
𝑑𝑢
𝐼2 = ∫ = ln|𝑢|
𝑢
D’où
1 + √1 + 𝑢2 √1 + 𝑢2 − 1 √1 + 𝑢2 − 1
∫ 𝑑𝑢 = 𝐼1 + 𝐼2 = ln | | + ln|𝑢| = ln | 𝑢| = ln |√1 + 𝑢2 − 1|
𝑢√1 + 𝑢2 𝑢 𝑢

Deuxième membre :
𝑑𝑦 𝐶
−∫ = − ln|𝑦| + 𝐶1 = ln | |
𝑦 𝑦

D’où la solution générale de l’EDO est :

𝐶 𝐶 𝑡
ln |√1 + 𝑢2 − 1| = ln | | ⇒ √1 + 𝑢2 − 1 = ; 𝑢 =
𝑦 𝑦 𝑦

𝑡2 𝐶 1 2 2−1=
𝐶
√1 + − 1 = ⇒ √𝑡 + 𝑦 ⇒ √𝑡 2 + 𝑦 2 − 𝑦 = 𝐶; 𝑦(1) = 0
𝑦2 𝑦 𝑦 𝑦
√12 + 02 − 0 = 𝐶 ⇒ 𝐶 = 1

D’où la solution finale :

 5

√𝑡 2 + 𝑦 2 − 𝑦 = 1

𝑡2 − 1
𝑡 2 + 𝑦 2 = (1 + 𝑦)2 = 1 + 2𝑦 + 𝑦 2 ⇒ 1 + 2𝑦 = 𝑡 2 ⇒ 𝑦 =
2

4. Déterminer toutes les fonctions f(x) pour lesquelles u(t)=t soit un facteur intégrant de l’équation
différentielle
𝑑𝑦
𝑓 (𝑡 ) ( ) + 𝑡2 + 𝑦 = 0
𝑑𝑡
Solution
𝑑𝑦
𝑢(𝑡) = 𝑡 est un facteur intégrant de l’EDO 𝑓(𝑡) + 𝑡 2 + 𝑦 = 0 si,
𝑑𝑡
𝑑𝑦 𝑑𝑦
𝑢(𝑡) [𝑓(𝑡) + 𝑡 2 + 𝑦] = 𝑡𝑓(𝑡) + 𝑡 3 + 𝑡𝑦 = 0 est une différentielle exacte. Se référant à la forme
𝑑𝑡 𝑑𝑡
générale 𝑀(𝑡, 𝑦)𝑑𝑦 + 𝑁(𝑡, 𝑦)𝑑𝑡 = 0, l’EDO peut se mettre sous la forme :

𝑡𝑓 (𝑡)𝑑𝑦 + 𝑡(𝑡 2 + 𝑦)𝑑𝑡 = 0

Avec 𝑀(𝑡, 𝑦) = 𝑡𝑓 (𝑡) et 𝑁(𝑡, 𝑦) = 𝑡(𝑡 2 + 𝑦) = 𝑡 3 + 𝑡𝑦.

𝜕𝑀 𝜕𝑁
La différentielle étant totale alors la condition = est vérifiée.
𝜕𝑡 𝜕𝑦
𝜕𝑀 𝜕(𝑡𝑓 (𝑡))
= = 𝑓(𝑡) + 𝑡𝑓 ′ (𝑡)
𝜕𝑡 𝜕𝑡
𝜕𝑁 𝜕(𝑡 3 + 𝑡𝑦)
= =𝑡
𝜕𝑦 𝜕𝑦

𝜕𝑀 𝜕𝑁
= ⇒ 𝑓(𝑡) + 𝑡𝑓 ′ (𝑡) = 𝑡
𝜕𝑡 𝜕𝑦
Pour trouver 𝑓 (𝑡), il faut résoudre l’EDO

𝑡𝑓 ′ (𝑡) + 𝑓(𝑡) = 𝑡

𝑓 ′ (𝑡) 1 𝑑𝑓(𝑡) 𝑑𝑡
Partie homogène : 𝑡𝑓 ′ (𝑡) + 𝑓(𝑡) = 0 ⇒ 𝑡𝑓 ′ (𝑡) = −𝑓(𝑡) ⇒ =− ⇒ =−
𝑓(𝑡) 𝑡 𝑓(𝑡) 𝑡
𝐶 𝐶
⇒ ln 𝑓(𝑡) = − ln 𝐶1 𝑡 = ln ⇒ 𝑓(𝑡) =
𝑡 𝑡
Par identification la solution particulière est de la forme 𝑓(𝑡) = 𝑎𝑡 + 𝑏 ⇒ 𝑓 ′ (𝑡) = 𝑎. En substituant
dans l’EDO, on a :
𝑎𝑡 + 𝑎𝑡 + 𝑏 = 𝑡 ⇒ 2𝑎𝑡 + 𝑏 = 𝑡
1
⇒ 𝑎 = ;𝑏 = 0
2
1
Il vient que la solution particulière est 𝑓 (𝑡) = 𝑡
2
D’où les fonctions f(t) sont :

𝐶 1
𝑓 (𝑡 ) = + 𝑡
𝑡 2
 6

QUESTION II. Equations différentielles linéaires à coefficients non constants.

1. Trouver la solution générale de

a) 𝒙𝟐 𝒚′′ − 𝒙𝒚′ + 𝒚 = 𝟎
Solution
Équation indicielle : 𝑟(𝑟 − 1) − 𝑟 + 1 = 0 → (𝑟 − 1)2 = 0 → 𝑟 = 1 (𝑅𝑎𝑐𝑖𝑛𝑒 𝑑𝑜𝑢𝑏𝑙𝑒)
𝑆𝐹𝑆 = {𝑦1 = 𝑥, 𝑦2 = 𝑥 ln 𝑥 }

𝑦 = 𝐶1 𝑥 + 𝐶2 𝑥 ln 𝑥

b) 𝟐𝒙𝟐 𝒚′′ + 𝟏𝟎𝒙𝒚′ + 𝟗𝒚 = 𝟎

Solution

Équation indicielle : 2𝑟(𝑟 − 1) + 10𝑟 + 9 = 0 → 2𝑟 2 + 8𝑟 + 9 = 0 ∆ = −8 √∆= 𝑖2√2

−8+𝑖√8 𝑖√2
𝑟1 = = −2 +
4 2
{
−8−𝑖√8 𝑖√2
𝑟1 = = −2 −
4 2

ln 𝑥 ln 𝑥
𝐶1 cos ( ) + 𝐶2 sin ( )
𝑦= √2 √2
𝑥2

c) 𝒙𝟐 𝒚′′ − 𝟕𝒙𝒚′ + 𝟖𝒚 = 𝟎
Solution

Équation indicielle : 𝑟(𝑟 − 1) − 7𝑟 + 8 = 0 → 𝑟 2 − 8𝑟 + 8 ∆ = 32 √∆= 4√2

8+4√2
𝑟1 = = 2(2 + √2)
2
{
8−4√2
𝑟1 = = 2(2 − √2)
2

𝑦 = 𝐶1 𝑥 2(2+√2) + 𝐶2 𝑥 2(2−√2)

d) 𝟏𝟐𝒙𝟐 𝒚” – 𝟓𝒙𝒚′ + 𝟔𝒚 = 𝟎
Solution

Équation indicielle : 12𝑟 (𝑟 − 1) − 5𝑟 + 6 = 0 → 12𝑟 2 − 17𝑟 + 6 ∆ = 1 √∆= 1

 7

17+1 3
𝑟1 = =
24 4
{ 17−1 2
𝑟1 = =
24 3

D’où la solution à l’EDO :

3 2
𝑦 = 𝐶1 𝑥 4 + 𝐶2 𝑥 3

2. Trouver le SFS de Frobenius pour chacune des équations ci-après :

a) 𝟒𝒙𝟐 (𝟏 + 𝒙)𝒚′′ + 𝟒𝒙(𝟏 + 𝟐𝒙)𝒚′ − (𝟏 + 𝟑𝒙)𝒚 = 𝟎
Solution
Cette équation peut s’écrire aussi sous la forme de :
𝑥 2 (4 + 4𝑥 )𝑦 ′′ + 𝑥 (4 + 8𝑥 )𝑦 ′ − (1 + 3𝑥 )𝑦 = 0
𝛼0 = 𝛼1 = 4, 𝛼2 = 0, 𝛽0 = 4, 𝛽1 = 8, 𝛽2 = 0, 𝛾0 = −1, 𝛾1 = −3, 𝛾2 = 0
𝑃0 (𝑟) = 𝛼0 𝑟(𝑟 − 1) + 𝛽0 𝑟 + 𝛾0 = 4𝑟 2 − 1
𝑃1 (𝑟) = 𝛼1 𝑟(𝑟 − 1) + 𝛽1 𝑟 + 𝛾1 = −4𝑟 2 − 4𝑟 − 3
𝑃2 (𝑟) = 𝛼2 𝑟(𝑟 − 1) + 𝛽2 𝑟 + 𝛾2 = 0
1 1
𝑃0 (𝑟) = 0 ⇒ 4𝑟 2 − 1 = 0 → 𝑟1 = 𝑒𝑡 𝑟2 = − 𝑟1 − 𝑟2 = 1 𝑃0 (𝑟1 ) = 0
2 2

1 1 2
𝑃0 (𝑛 + 𝑟) = 𝑃0 (𝑛 + ) = 4 (𝑛 + ) − 1 ≠ 0
2 2
∞ ∞
𝑟1 𝑛
1 1
𝑦1 = 𝑥 ∑ 𝑎𝑛 (𝑟1 )𝑥 = 𝑥2 ∑ 𝑎𝑛 ( ) 𝑥 𝑛
2
𝑛=0 𝑛=0

𝑎0 (𝑟 ) = 1
𝑃1 (𝑟+𝑛−1)
𝑎𝑛 (𝑟 ) = − 𝑎𝑛−1 (𝑟)
𝑃0 (𝑟+𝑛)

𝑃1 (𝑟 + 𝑛 − 1) = 4(𝑟 + 𝑛 − 1)(𝑟 + 𝑛) − 3
𝑃0 (𝑟 + 𝑛) = 4(𝑟 + 𝑛)2 − 1
4(𝑟 + 𝑛 − 1)(𝑟 + 𝑛) − 3
𝑎𝑛 (𝑟 ) = − 𝑎𝑛−1 (𝑟)
4 (𝑟 + 𝑛 )2 − 1
En développant on trouve :
𝑛
4(𝑟 + 𝑘 − 1)(𝑟 + 𝑘) − 3
𝑎𝑛 (𝑟) = (−1 )𝑛 ∏
4 (𝑟 + 𝑘 )2 − 1
𝑘=1
∞ 𝑛
1 4(𝑟 + 𝑘 − 1)(𝑟 + 𝑘) − 3 𝑛
𝑦1 (𝑥 ) = 𝑥2 ∑ (−1 )𝑛 ∏ 𝑥
4 (𝑟 + 𝑘 ) 2 − 1
𝑛=0 𝑘=1

𝑃1 (𝑟1 − 1)
𝐶= − 𝑎𝑘−1 (𝑟2 ) 𝑃1 (𝑟1 − 1) = −4 𝑘 = 1 𝛼0 = 4 𝑎𝑘−1 (𝑟2 ) = 𝑎0 (𝑟2 ) = 1
𝑘𝛼0
𝐶=1
 8

Trouvons 𝑎𝑛′ (𝑟)

𝑛
4(𝑟 + 𝑘) + 4(𝑟 + 𝑘 − 1) 8(𝑟 + 𝑘)
𝑎𝑛′ (𝑟) = 𝑎𝑛 (𝑟 ) ∑ −
4(𝑟 + 𝑘 − 1)(𝑟 + 𝑘) − 3 4(𝑟 + 𝑘)2 − 1
𝑘=1

𝑛 1 1 1
4 (𝑘 + ) + 4 (𝑘 − ) 8 (𝑘 + )
𝑎𝑛′ (𝑟1 ) = 𝑎𝑛 (𝑟1 ) ∑ 2 2 − 2
1 1 1 2
𝑘=1 4 (𝑘 − ) (𝑘 + ) − 3 4 (𝑘 + ) − 1
2 2 2
𝑛 3 1
4 (𝑘 − ) (𝑘 − ) − 3
𝑎𝑛 (𝑟2 ) = (−1)𝑛 ∏ 2 2
2
1
𝑘=1 4 (𝑘 − ) − 1
2
𝑘−1 ∞
𝑟2 𝑛 𝑟1
𝑦2 = 𝑥 ∑ 𝑎𝑛 (𝑟2 )𝑥 + 𝐶(𝑦1 (𝑥 ) ln 𝑥 + 𝑥 ∑ 𝑎𝑛′ (𝑟1 )𝑥 𝑛
𝑛=0 𝑛=1

Alors

𝑦(𝑥 ) = 𝐶1 𝑦1 (𝑥 ) + 𝐶2 𝑦2 (𝑥)

b) 𝒙𝒚′′ + 𝒙𝒚′ + 𝒚 = 𝟎
Solution
Cette équation peut encore s’écrire, en multipliant membre à membre par 𝑥 :
𝑥 2 𝑦 ′′ + 𝑥 2 𝑦 ′ + 𝑥𝑦 = 0
Il vient que, par identification à la forme générale de l’équation de Frobenius
(𝑥 − 𝑥0 )2 𝐴(𝑥 )𝑦 ′′ + (𝑥 − 𝑥0 )𝐵 (𝑥 )𝑦 ′ + 𝐶 (𝑥 )𝑦 = 0
Avec
𝑥0 = 0, 𝐴(𝑥 ) = 1, 𝐵 (𝑥 ) = 𝑥 = 𝐶 (𝑥 )
⇒ 𝛼0 = 1, 𝛼1 = 𝛼2 = 0, 𝛽0 = 𝛽2 = 0, 𝛽1 = 1, 𝛾0 = 𝛾2 = 0, 𝛾1 = 1
𝑝0 (𝑟) = 𝛼0 𝑟(𝑟 − 1) + 𝛽0 𝑟 + 𝛾0 = 𝑟(𝑟 − 1)
𝑝1 (𝑟) = 𝛼1 𝑟(𝑟 − 1) + 𝛽1 𝑟 + 𝛾1 = 𝑟 + 1
𝑝2 (𝑟) = 𝛼2 𝑟(𝑟 − 1) + 𝛽2 𝑟 + 𝛾2 = 0
𝑟1 = 1
𝑝0 (𝑟) = 0 ⇒ 𝑟(𝑟 − 1) = 0 ⇒ {
𝑟2 = 0
𝑟1 = 𝑟2 + 𝑘 ⇒ 𝑘 = 1
𝑝0 (𝑟 + 𝑛) = (𝑟 + 𝑛)(𝑟 + 𝑛 − 1) = (𝑟 + 𝑛)2 − (𝑟 + 𝑛) ≠ 0
𝑎0 (𝑟 ) = 1
𝑝1 (𝑟 + 𝑛 − 1) 1
𝑎𝑛 (𝑟 ) = − 𝑎𝑛−1 (𝑟) = − 𝑎 (𝑟 )
𝑝0 (𝑟 + 𝑛) 𝑟 + 𝑛 − 1 𝑛−1
 9

1
𝑎1 = − 𝑎 (𝑟 )
𝑟+1−1 0
1 1
𝑎2 = − (− 𝑎 (𝑟 ) )
𝑟+2−1 𝑟+1−1 0
⋮
𝑛
1 1
𝑎𝑛 (𝑟) = (−1 )𝑛 ∏ = (−1)𝑛
𝑟+𝑘−1 (𝑟 + 𝑛 − 1)!
𝑘=1

Ainsi ;
∞ ∞ ∞
𝑥𝑛
𝑦1 (𝑥 ) = 𝑥 𝑟1 ∑ 𝑎𝑛 (𝑟1 )𝑥 𝑛 = 𝑥 0 ∑ 𝑎𝑛 (0)𝑥 𝑛 = ∑ (−1)𝑛
(𝑛 − 1)!
𝑛=0 𝑛=0 𝑛=0

∞ ∞
𝑟2 𝑛 𝑟1
𝑝1 (𝑟1 − 1)
𝑦2 (𝑥 ) = 𝑥 ∑ 𝑎𝑛 (𝑟2 )𝑥 + 𝐶 (𝑦1 (𝑥 ) ln 𝑥 + 𝑥 ∑ 𝑎𝑛′ (𝑟1 )𝑥 𝑛 ) , 𝐶=− 𝑎𝑘−1 (𝑟2 )
𝑘𝛼𝑜
𝑛=0 𝑛=1

𝑝1 (𝑟2 + 𝑛 − 1)
𝑎0 (𝑟2 ) = 1, 𝑎𝑛 (𝑟2 ) = − 𝑎𝑛−1 (𝑟2 ), 1≤𝑛 ≤𝑘−1
𝑝0 (𝑟2 + 𝑛)
Recherche de 𝑎𝑛′ (𝑟)
𝑛 𝑛 𝑛
1 1 1
ln|𝑎𝑛 (𝑟)| = ln |(−1 )𝑛 ∏ | = ln |∏ | = ∑ ln | |
𝑟+𝑘−1 𝑟+𝑘−1 𝑟+𝑘−1
𝑘=1 𝑘=1 𝑘=1
𝑛

= − ∑ ln|𝑟 + 𝑘 − 1|
𝑘=1

En dérivant membre à membre on a :

𝑛 𝑛 𝑛
𝑎𝑛′ (𝑟) 1 1 1
= −∑ ⇒ 𝑎𝑛′ (𝑟) = −𝑎𝑛 (𝑟) ∑ ⇒ 𝑎𝑛′ (𝑟) = 𝑎𝑛 (𝑟) ∑
𝑎𝑛 (𝑟 ) 𝑟+𝑘−1 𝑟+𝑘−1 1 − (𝑟 + 𝑘 )
𝑘=1 𝑘=1 𝑘=1

Ainsi,
𝑝1 (𝑟1 − 1) 𝑟1 − 1 + 1 1−1+1
𝐶=− 𝑎𝑘−1 (𝑟2 ) = − 𝑎0 (𝑟2 ) = − 1=1
𝑘𝛼𝑜 1.1 1
Comme 𝑘 = 1, alors 𝑎𝑛 (𝑟2 ) = 0 étant donné que l’on doit avoir 1 ≤ 𝑛 ≤ 𝑘 − 1.
∞ ∞
𝑟2 𝑛 𝑟1
𝑦2 (𝑥 ) = 𝑥 ∑ 𝑎𝑛 (𝑟2 )𝑥 + 𝐶 (𝑦1 (𝑥 ) ln 𝑥 + 𝑥 ∑ 𝑎𝑛′ (𝑟1 )𝑥 𝑛 )
𝑛=0 𝑛=1
∞ 𝑛
1
= 0 + (𝑦1 (𝑥 ) ln 𝑥 + 𝑥 ∑ 𝑎𝑛 (𝑟1 ) ∑ 𝑥𝑛)
1 − (𝑟1 + 𝑘)
𝑛=1 𝑘=1
𝑛 𝑛
1 1
= 𝑦1 (𝑥 ) ln 𝑥 − 𝑥 ∑ (−1)𝑛 ∑ 𝑥𝑛
𝑛! 𝑘
𝑘=1 𝑘=1
∞ 𝑛 𝑛
𝑛
𝑥 1 1
= ∑ (−1 )𝑛 ln 𝑥 − 𝑥 ∑ (−1)𝑛 ∑ 𝑥 𝑛
(𝑛 − 1)! 𝑛! 𝑘
𝑛=0 𝑘=1 𝑘=1
 10

D’où le système fondamental de solution :

∞ ∞ 𝑛 𝑛
𝑛
𝑥𝑛 𝑛
𝑥𝑛 1 1
(
𝑆𝐹𝑆 = {𝑦1 = ∑ −1) , (
𝑦2 = ∑ −1) ln 𝑥 − 𝑥 ∑ (−1)𝑛 ∑ 𝑥 𝑛 }
(𝑛 − 1)! (𝑛 − 1)! 𝑛! 𝑘
𝑛=0 𝑛=0 𝑘=1 𝑘=1

3. Déterminer la fonction de green ainsi que la solution pour le PVB suivant :

a) 𝒚′′ − 𝒚′ = 𝒙 𝒚(𝟎) + 𝒚′ (𝟎) = 𝟑 𝒚(𝟏) + 𝒚′ (𝟏) = 𝟎
Solution
𝑃0 (𝑥 )𝑦 ′′ + 𝑃1 (𝑥 )𝑦 ′ + 𝑃2 (𝑥 )𝑦 = 𝐹(𝑥)
Forme { 𝛼𝑦(𝑎) + 𝛽𝑦 ′ (𝑎) = 𝑘1
𝛾𝑦(𝑏) + 𝛿𝑦 ′ (𝑏) = 𝑘2
Pour notre cas : 𝛼 = 1 𝛽= 1 𝛾 = 1 𝛿 = −1 𝑘1 = 3 𝑘2 = 0
𝐵1 (𝑦) = 3 𝐵2 (𝑦) = 0
𝑦 ′′ − 𝑦 = 0
Le PVB devient {𝐵1 (𝑦) = 3
𝐵2 (𝑦) = 0
𝑍1 (𝑥 ) = 𝑒 𝑥 𝑍2 (𝑥 ) = 𝑒 −𝑥 𝑒𝑠𝑡 𝑢𝑛 𝑆𝐹𝑆 𝑑𝑒 𝑦 ′′ − 𝑦 = 0
2
𝐵1 (𝑍1 ) = 2 𝐵1 (𝑍2 ) = 0 𝐵2 (𝑍1 ) = 0 𝐵2 (𝑍2 ) =
𝑒

𝐵 (𝑍 ) 𝐵2 (𝑍1 ) 4
| 1 1 |= ≠ 0
𝐵1 (𝑍2 ) 𝐵2 (𝑍2 ) 𝑒
Alors l’EDO admet une solution unique
𝑦1 (𝑥 ) = 𝐵1 (𝑍2 )𝑍1 (𝑥 ) − 𝐵1 (𝑍1 )𝑍2 (𝑥 ) = −2𝑒 −𝑥
𝑦2 (𝑥 ) = 𝐵2 (𝑍2 )𝑍1 (𝑥 ) − 𝐵2 (𝑍1 )𝑍2 (𝑥 ) = 2𝑒 𝑥−1

𝑦1 𝑦2 8
𝑊 = |𝑦 ′ 𝑦2′ | = − 2 ≠ 0
1

{𝑦1 , 𝑦2 } 𝑒𝑠𝑡 𝑢𝑛 𝑆𝐹𝑆

𝑏
𝑘2 𝑘1
𝑦(𝑥 ) = ∫ 𝐺 (𝑥, 𝑡)𝐹 (𝑡) 𝑑𝑡 + 𝑦1 (𝑥 ) + 𝑦 (𝑥 )
𝐵2 (𝑦1 ) 𝐵1 (𝑦2 ) 2
𝑎

𝑦1 (𝑡)𝑦2 (𝑥 )
0≤𝑡 ≤𝑥
𝑃0 (𝑡)𝑊 (𝑡)
𝐺 (𝑥, 𝑡) =
𝑦1 (𝑥 )𝑦2 (𝑡)
𝑥≤𝑡 ≤1
{ 𝑃0 (𝑡)𝑊 (𝑡)
 11

∫ 𝐺 (𝑥, 𝑡)𝐹 (𝑡) 𝑑𝑡

𝑎
𝑥 1 𝑥
𝑦1 (𝑡)𝑦2 (𝑥 )
= ∫ 𝐺 (𝑥, 𝑡)𝐹 (𝑡)𝑑𝑡 + ∫ 𝐺 (𝑥, 𝑡)𝐹 (𝑡)𝑑𝑡 = ∫ 𝐹 (𝑡)𝑑𝑡
𝑃0 (𝑡)𝑊 (𝑡)
0 𝑥 0
1 𝑥 𝑥
𝑥
𝑦1 (𝑥 )𝑦2 (𝑡) 𝑒 𝑒 −𝑥 1
+ ∫ 𝐹 (𝑡)𝑑𝑡 = ∫ 𝑡𝑒 −𝑡 𝑑𝑡 + ∫ 𝑡𝑒 𝑡 𝑑𝑡 = − 𝑒 𝑥 − 𝑥
𝑃0 (𝑡)𝑊 (𝑡) 2 2 2
𝑥 0 0

𝑘2
𝑦 (𝑥 ) = 0
𝐵2 (𝑦1 ) 1
𝑘1 3
𝑦2 (𝑥 ) = 𝑒 𝑥
𝐵1 (𝑦2 ) 2
1 3
𝑦 (𝑥 ) = − 𝑒 𝑥 − 𝑥 + 𝑒 𝑥
2 2
𝑦 (𝑥 ) = 𝑒 𝑥 − 𝑥

4. Résoudre le problème de Sturm-Liouville ci-après :

a) 𝒚′′ + 𝟐𝒚′ + 𝒚 + 𝝀𝒚 = 𝟎 𝒚′ (𝟎) = 𝟎, 𝒚′ (𝟏) = 𝟎 (𝟏)
Solution
Posons u(x) = 2, p′ (x) = p(x)u(x) p(x) = e2x
En multipliant (1) par p(x) le problème devient :
(𝑒 2𝑥 𝑦 ′ )′ + 𝑒 2𝑥 𝑦 + 𝜆𝑦 = 0 𝑦 ′ (0) = 0 𝑦 ′ (1) = 0 𝑆−𝐿
Équation caractéristique : 𝑟 2 + 2𝑟 + 1 + 𝜆 = 0
𝑃𝑜𝑢𝑟 𝜆 < 0, 𝜆 = −𝑘 2 , ∆= 4𝑘 2 > 0
𝑟1 = 𝑘 − 1
{
𝑟2 = −𝑘 − 1
𝑦 = 𝐶1 𝑒 (𝑘−1)𝑥 + 𝐶2 𝑒 (−𝑘−1)𝑥

𝑦′ = 𝐶1 (𝑘 − 1)𝑒 (𝑘−1)𝑥 + 𝐶2 (−𝑘 − 1)𝑒 (−𝑘−1)𝑥

C.B :
𝐶1 (𝑘 − 1) + 𝐶2 (−𝑘 − 1) = 0

𝐶1 (𝑘 − 1)𝑒 (𝑘−1) + 𝐶2 (−𝑘 − 1)𝑒 (−𝑘−1) = 0

𝑘−1 −𝑘 − 1
| | = (𝑘 2 − 1)(𝑒 (𝑘−1) − 𝑒 (−𝑘−1) ) ≠ 0
(𝑘 − 1)𝑒 (𝑘−1) (−𝑘 − 1)𝑒 (−𝑘−1)
D’où C1 = 0 et C2 = 0

𝑃𝑜𝑢𝑟 𝜆 = 0, ∆= 0
 12

𝑟 = −1 (𝑟𝑎𝑐𝑖𝑛𝑒 𝑑𝑜𝑢𝑏𝑙𝑒)
𝑦 = 𝐶1 𝑒 −𝑥 + 𝐶2 𝑥𝑒 −𝑥
𝑦 ′ = −𝐶1 𝑒 −𝑥 + 𝐶2 (𝑒 −𝑥 − 𝑥𝑒 −𝑥 )
C.B :
−𝐶1 + 𝐶2 = 0
−𝐶1 𝑒 −1 + 𝐶2 (𝑒 −1 − 𝑒 −1 ) = 0
D’où C1 = 0 et C2 = 0
𝑃𝑜𝑢𝑟 𝜆 > 0, 𝜆 = 𝑘2, ∆= −4𝑘 2 < 0

𝑟1 = −1 + 𝑖𝑘
{
𝑟2 = −1 − 𝑖𝑘
𝑦(𝑥 ) = 𝐶1 𝑒 −𝑥 cos 𝑘𝑥 + 𝐶2 𝑒 −𝑥 sin 𝑘𝑥

𝑦 ′(𝑥) = 𝐶1 (−𝑒 −𝑥 cos 𝑘𝑥 − 𝑘𝑒 −𝑥 sin 𝑘𝑥 ) + 𝐶2 (−𝑒 −𝑥 sin 𝑘𝑥 + 𝑒 −𝑥 cos 𝑘𝑥)

C.B :
−𝐶1 + 𝐶2 = 0 → 𝐶1 = 𝐶2
{
𝐶1 (cos 𝑘 + sin 𝑘 ) = 0
On prend le système
𝐶1 = 𝐶2
{
(cos 𝑘 + sin 𝑘) = 0
(cos 𝑘 + sin 𝑘 ) = 0 → cos 𝑘 = sin 𝑘
𝜋 𝜋 2
𝑘= + 𝜋𝑛 , 𝜆 = 𝑘 2 , 𝜆 = ( + 𝜋𝑛 )
4 4

𝜋
𝑦(𝑥 ) = 𝐶2 𝑒 −𝑥 sin (( + 𝜋𝑛) 𝑥)
4

QUESTION III. Equations aux dérivées partielles

1. Résoudre le PVB ci-après par la méthode des caractéristiques.

𝒖𝒕 − 𝒖𝒖𝒙 = 𝟎, 𝒕 > 𝟎, 𝒙 ∈ 𝑹
a) { 𝝅
𝒖(𝟎, 𝒙) = − 𝐚𝐭𝐚𝐧 𝒙 , 𝒙 ∈ 𝑹
𝟐

Solution
Le système fondamental d’EDO associé à ce problème est :
𝑑𝑡 𝑑𝑥 𝑑𝑢
= =
1 −𝑢 0
𝑑𝑡 𝑢𝑑𝑥 𝑡𝑢𝑑𝑢 𝑑𝑡 − 𝑢𝑑𝑥 − 𝑡𝑢𝑑𝑢
⇒ = = =
1 −𝑢2 0 1 + 𝑢2
𝑑𝑡 𝑑𝑡−𝑢𝑑𝑥−𝑡𝑢𝑑𝑢
L’EDO = donne :
1 1+𝑢2
 13

𝑑𝑡 𝑑𝑡 − 𝑢𝑑𝑥 − 𝑡𝑢𝑑𝑢 𝑑𝑡 𝑑𝑡 𝑢 𝑡𝑢𝑑𝑢 𝑢2 𝑑𝑡

= ⇒ − = − 𝑑𝑥 − ⇒
1 1 + 𝑢2 1 1 + 𝑢2 1 + 𝑢2 1 + 𝑢2 1 + 𝑢2
𝑢 𝑡𝑢𝑑𝑢
=− 2
𝑑𝑥 −
1+𝑢 1 + 𝑢2

1+𝑢2
En multipliant membre à membre par on a :
𝑢

𝑢𝑑𝑡 = −𝑑𝑥 − 𝑡𝑑𝑢 ⇒ 𝑢𝑑𝑡 + 𝑡𝑑𝑢 = −𝑑𝑥 ⇒ 𝑑 (𝑢𝑡) = −𝑑𝑥 ⇒ 𝑢𝑡 = −𝑥 + 𝑐1 ⇒ 𝑐1 = 𝑢𝑡 + 𝑥

D’où ; 𝑐1 = 𝑢𝑡 + 𝑥 = Φ(𝑥, 𝑦, 𝑢)

𝑑𝑡 𝑑𝑢
L’EDO = ⇒ 0𝑑𝑡 = 𝑑𝑢 a pour solution 𝑐2 = 𝑢 = Ψ(𝑥, 𝑦, 𝑢)
1 0

Ainsi, la solution au problème est :

𝑓 (Ψ, Φ) = 0 = Ψ − 𝑔(Φ) ⇒ 𝑢(𝑡, 𝑥 ) = 𝑔(𝑢𝑡 + 𝑥 ) ⇒ 𝑢(𝑡, 𝑥 ) = 𝑔(𝑢𝑡 + 𝑥 )
𝜋 𝜋
Tenant compte de la condition initiale : 𝑢(0, 𝑥 ) = − atan 𝑥 ⇒ 𝑔(𝑥 ) = − atan 𝑥
2 2

D’où :
𝜋
𝑢(𝑡, 𝑥 ) = 𝑔 ( − atan(𝑢𝑡 + 𝑥 ))
2

2. Résoudre le problème aux valeurs initiales et des bords ci-après par la méthode de séparation des
variables :
𝒖𝒕 = 𝟗𝒖𝒙𝒙 , 𝒕 > 𝟎, 𝟎 < 𝒙 < 𝟒
a) {𝒖(𝒕, 𝟎) = 𝟎, 𝒖(𝒕, 𝟒) = 𝟎, 𝒕 > 𝟎
𝒖𝒙 (𝟎, 𝒙) = 𝟏, 𝟎≤𝒙≤𝟒
Solution
Supposons que la solution non triviale 𝑢 de ce problème peut se mettre sous la forme :

𝑢(𝑡, 𝑥 ) = 𝑇(𝑡)𝑋(𝑥 )
⇒ 𝑢𝑡 = 𝑇 ′ (𝑡)𝑋(𝑥 ), 𝑢𝑥𝑥 = 𝑇(𝑡)𝑋 ′′ (𝑥 )
En substituant dans l’équation on trouve :
𝑇 ′ (𝑡)𝑋(𝑥 ) = 9𝑇(𝑡)𝑋 ′′ (𝑥 )
En divisant l’équation par 9𝑇(𝑡)𝑋(𝑡) on trouve :
𝑇 ′ (𝑡) 𝑋 ′′ (𝑥 )
= =𝜇
9𝑇(𝑡) 𝑋 (𝑥 )
D’où le système
𝑇 ′ (𝑡) − 9𝜇𝑇(𝑡) = 0
{
𝑋 ′′ (𝑥 ) − 𝜇𝑋(𝑥 ) = 0

Conditions initiales : puisque la solution 𝑢(𝑡, 𝑥 ) = 𝑇(𝑡)𝑋(𝑥 ) ne doit pas être triviale alors on a :
 14

𝑢(𝑡, 0) = 0 ⇒ 𝑇(𝑡)𝑋(0) = 0 ⇒ 𝑋(0) = 0 et 𝑢(𝑡, 4) = 0 ⇒ 𝑇(𝑡)𝑋(4) = 0 ⇒ 𝑋(4) = 0

D’où le système d’EDO :

𝑇 ′ (𝑡) − 9𝜇𝑇(𝑡) = 0
{ 𝑋 ′′ (𝑥 ) − 𝜇𝑋(𝑥 ) = 0
𝑋(0) = 0, 𝑋(4) = 0, 𝑡 > 0

Le deuxième problème 𝑋 ′′ (𝑥 ) − 𝜇𝑋 (𝑥 ) = 0 avec les conditions 𝑋(0) = 0, 𝑋(4) = 0 est un

problème de S-L :
En posant −𝜇 = 𝜆, l’équation caractéristique associée est 𝑟 2 + 𝜆 = 0 dont le discriminant est Δ =
−4𝜆.
• Si Δ = 0 ⇒ λ = 0 = −μ ⇒ X ′′ (x) = 0 ⇒ 𝑋(𝑥 ) = 𝑐1 𝑥 + 𝑐2
Conditions de bord :
𝑋(0) = 0 ⇒ 𝑐2 = 0
{
𝑋(4) = 0 ⇒ 𝑐1 = 0
La solution est triviale.

• Si Δ > 0 ⇒ −4𝜆 > 0 ⇒ 𝜆 < 0, nous avons deux racines réelles distinctes :
𝑟1,2 = ±√−𝜆, 𝜆<0
La solution est :

𝑋(𝑥 ) = 𝑐1 𝑒 √−𝜆𝑥 + 𝑐2 𝑒 −√−𝜆𝑥

Conditions de bord :
𝑋(0) = 0 ⇒ 𝑐1 + 𝑐2 = 0
{
𝑋(4) = 0 ⇒ 𝑐1 𝑒 4√−𝜆 + 𝑐2 𝑒 −4√−𝜆 = 0
1 1
Comme | 4√−𝜆 −4√−𝜆
| = 𝑒 −4√−𝜆 − 𝑒 4√−𝜆 = 𝑒 −4√−𝜆 (1 − 𝑒 8√−𝜆 ) ≠ 0 alors la solution est
𝑒 𝑒
triviale.

• Si Δ < 0 ⇒ −4𝜆 < 0 ⇒ 𝜆 > 0, nous avons deux racines complexes conjuguées :
𝑟1,2 = ±𝑖√𝜆, 𝜆>0
La solution est :
𝑋(𝑥 ) = 𝑐1 cos √𝜆𝑥 + 𝑐2 sin √𝜆𝑥

Conditions de bord :
𝑋(0) = 0 ⇒ 𝑐1 = 0
{
𝑋(4) = 0 ⇒ 𝑐2 sin 4√𝜆 = 0
 15

Pour la deuxième équation du système, la solution 𝑐2 = 0 sera rejetée car elle conduit à une
solution triviale pour 𝑋(𝑥 ). Par contre, sin 4√𝜆 = 0 permet d’écrire :
𝑘𝜋 2
4√𝜆 = 𝑘𝜋 ⇒ 𝜆𝑘 = ( )
4
Ainsi la solution est :

𝑘𝜋 𝑘𝜋 2
{𝑋𝑘 (𝑥 ) = 𝑐2 sin √𝜆𝑘 𝑥 = 𝑐2 sin ( 𝑥) , 𝜆𝑘 = −𝜇𝑘 = ( ) , 𝑘, 𝑥, 𝑐2 ≠ 0}
4 4

Le premier problème 𝑇 ′ (𝑡) − 9𝜇𝑇(𝑡) = 0 a pour solution 𝑇(𝑡) = 𝑐𝑒 9𝜇𝑡

𝑘𝜋 2
𝑘𝜋 2 −9( ) 𝑡
Comme 𝜇𝑘 = −𝜆𝑘 = − ( ) alors on a : 𝑇𝑘 (𝑡) = 𝑐𝑒 4
4

Ainsi, on peut déterminer 𝑢(𝑡, 𝑥 ) ;

2
𝑘𝜋𝑥 −9(𝑘𝜋) 𝑡
𝑢𝑘 (𝑡, 𝑥 ) = 𝐶 sin ( )𝑒 4
4
Par le principe de la superposition des solutions ;
∞
𝑘𝜋 2 𝑘𝜋𝑥
−9( ) 𝑡
𝑢(𝑡, 𝑥 ) = ∑ 𝑏𝑘 𝑒 4 sin ( ), 𝑡 > 0, 0<𝑥<4
4
𝑘=1

Où 𝑏𝑘 sont les coefficients de Fourier en sinus donnés par :

2 4 𝑘𝜋𝑥
𝑏𝑘 = ∫ 𝑢(0, 𝑥 ) sin ( ) 𝑑𝑥
4 0 4
Avec 𝑢𝑥 (0, 𝑥 ) = 1 ⇒ 𝑢(0, 𝑥 ) = 𝑥, on a :
1 4 𝑘𝜋𝑥
𝑏𝑘 = ∫ 𝑥 sin ( ) 𝑑𝑥
2 0 4
𝑎 = 𝑥 ⇒ 𝑑𝑎 = 𝑑𝑥
𝑘𝜋𝑥 4 𝑘𝜋𝑥
𝑣 = ∫ sin ( ) 𝑑𝑥 = − cos ( )
4 𝑘𝜋 4
1 4𝑥 𝑘𝜋𝑥 4 4 4 𝑘𝜋𝑥 1 4𝑥 𝑘𝜋𝑥 4 8
𝑏𝑘 = {[− cos ( )] + ∫ cos ( ) 𝑑𝑥 } = [− cos ( )] = − cos 𝑘𝜋
2 𝑘𝜋 4 0 𝑘𝜋 0 4 2 𝑘𝜋 4 0 𝑘𝜋
8
=− (−1)𝑘
𝑘𝜋
8
⇒ 𝑏𝑘 = − (−1)𝑘 , 𝑘 = 1,2,3, …
𝑘𝜋
D’où
∞ 2
8 𝑘𝜋
−9( ) 𝑡 𝑘𝜋𝑥
𝑢(𝑡, 𝑥 ) = − ∑ (−1)𝑘 𝑒 4 sin ( ),
𝑘𝜋 4
𝑘=1
 16

𝒖𝒕 = 𝒖𝒙𝒙 + 𝝅𝟐 𝐬𝐢𝐧(𝝅𝒙) , 𝒕 > 𝟎, 𝟎 < 𝒙 < 𝟏

b) { 𝒖(𝒕, 𝟎) = 𝟑, 𝒖𝒙 (𝒕, 𝟏) = −𝝅, 𝒕 > 𝟎
𝒖(𝟎, 𝒙) = 𝟐 𝐬𝐢𝐧(𝝅𝒙) , 𝟎≤𝒙≤𝟏
Solution
Il s’agit d’un problème non homogène qu’on doit essayer de transformer en problème homogène.
Pour ce faire posons 𝑢(𝑡, 𝑥 ) = 𝑣 (𝑡, 𝑥 ) + 𝑞 (𝑥 ) où 𝑞(𝑥 ) est une fonction à déterminer. On a :
𝑢𝑡 = 𝑣𝑡 , 𝑢𝑥𝑥 = 𝑣𝑥𝑥 + 𝑞 ′′ (𝑥 ), 𝑢𝑥 = 𝑣𝑥 + 𝑞 ′ (𝑥 )
Ainsi 𝑢(𝑡, 𝑥 ) est une solution au problème donné si 𝑣(𝑡, 𝑥 ) satisfait :
𝑣𝑡 = 𝑣𝑥𝑥 + 𝑞 ′′ (𝑥 ) + 𝜋 2 sin(𝜋𝑥 ) , 𝑡 > 0, 0 < 𝑥 < 1
{𝑣(𝑡, 0) + 𝑞 (0) = 3, 𝑣𝑥 (𝑡, 1) + 𝑞 ′ (1) = −𝜋, 𝑡 > 0
𝑣(0, 𝑥 ) + 𝑞(𝑥 ) = 2 sin(𝜋𝑥 ) , 0≤𝑥≤1
Ceci se réduit à :
𝑣𝑡 = 𝑣𝑥𝑥 , 𝑡 > 0, 0 < 𝑥 < 1
( )
{𝑣 𝑡, 0 = 0, 𝑣𝑥 (𝑡, 1) = 0, 𝑡 > 0
𝑣(0, 𝑥 ) + 𝑞 (𝑥 ) = 2 sin(𝜋𝑥 ) , 0≤𝑥≤1

A condition que 𝑞(𝑥 ) soit solution du PVB :

𝑞 ′′ (𝑥 ) + 𝜋 2 sin(𝜋𝑥 ) = 0, 𝑞(0) = 3, 𝑞 ′ (1) = −𝜋
On peut donc trouver 𝑞(𝑥 ) en intégrant deux fois :

𝑞 ′′ (𝑥 ) = −𝜋 2 sin(𝜋𝑥 ) ⇒ 𝑞 ′ (𝑥 ) = 𝜋 cos(𝜋𝑥 ) + 𝑐1 ⇒ 𝑞(𝑥 ) = sin(𝜋𝑥 ) + 𝑐1 𝑥 + 𝑐2

Tenant compte des conditions initiales on a :

𝑞(0) = 3 ⇒ 𝑐2 = 3, 𝑞 ′ (1) = −𝜋 ⇒ 𝜋 cos 𝜋 + 𝑐1 = −𝜋 ⇒ 𝑐1 = 0

Ainsi ; on a la fonction 𝑞(𝑥 ) = 3 + sin(𝜋𝑥 )

On doit donc chercher à résoudre le PVB :

𝑣𝑡 = 𝑣𝑥𝑥 , 𝑡 > 0, 0 < 𝑥 < 1
{ 𝑣(𝑡, 0) = 0, 𝑣𝑥 (𝑡, 1) = 0, 𝑡 > 0
𝑣(0, 𝑥 ) = −3 + sin(𝜋𝑥 ) , 0≤𝑥≤1
Supposons que la solution non triviale 𝑢 de ce problème peut se mettre sous la forme :
𝑣(𝑡, 𝑥 ) = 𝑇(𝑡)𝑋(𝑥 ) ⇒ 𝑣𝑡 = 𝑇 ′ (𝑡)𝑋(𝑥 ), 𝑣𝑥𝑥 = 𝑇(𝑡)𝑋 ′′ (𝑥 ), 𝑣𝑥 = 𝑇(𝑡)𝑋 ′ (𝑥 )
En substituant dans l’équation on trouve :
𝑇 ′ (𝑡)𝑋(𝑥 ) = 𝑇(𝑡)𝑋 ′′ (𝑥 )
En divisant l’équation par 𝑇(𝑡)𝑋(𝑡) on trouve :
𝑇 ′ (𝑡) 𝑋 ′′ (𝑥 )
= =𝜇
𝑇 (𝑡 ) 𝑋 (𝑥 )
D’où le système
𝑇 ′ (𝑡) − 𝜇𝑇(𝑡) = 0
{
𝑋 ′′ (𝑥 ) − 𝜇𝑋(𝑥 ) = 0
 17

Conditions initiales : puisque la solution 𝑣(𝑡, 𝑥 ) = 𝑇(𝑡)𝑋(𝑥 ) ne doit pas être triviale alors on a :
𝑣(𝑡, 0) = 0 ⇒ 𝑇(𝑡)𝑋(0) = 0 ⇒ 𝑋(0) = 0 et 𝑣𝑥 (𝑡, 1) = 0 ⇒ 𝑇(𝑡)𝑋 ′ (1) = 0 ⇒ 𝑋 ′ (1) = 0
D’où le système d’EDO :
𝑇 ′ (𝑡) − 𝜇𝑇(𝑡) = 0
{ 𝑋 ′′ (𝑥 ) − 𝜇𝑋(𝑥 ) = 0
𝑋(0) = 0, 𝑋 ′ (1) = 0, 𝑡 > 0

Le deuxième problème 𝑋 ′′ (𝑥 ) − 𝜇𝑋 (𝑥 ) = 0 avec les conditions 𝑋(0) = 0, 𝑋 ′ (1) = 0 est un

problème de S-L :
En posant −𝜇 = 𝜆, l’équation caractéristique associée est 𝑟 2 + 𝜆 = 0 dont le discriminant est Δ =
−4𝜆.
• Si Δ = 0 ⇒ λ = 0 = −μ ⇒ X ′′ (x) = 0 ⇒ 𝑋(𝑥 ) = 𝑐1 𝑥 + 𝑐2
𝑋(0) = 0 ⇒ 𝑐2 = 0
Conditions de bord : {
𝑋 ′ (1) = 0 ⇒ 𝑐1 = 0
La solution est triviale.

• Si Δ > 0 ⇒ −4𝜆 > 0 ⇒ 𝜆 < 0, nous avons deux racines réelles distinctes :
𝑟1,2 = ±√−𝜆, 𝜆<0
La solution est :

𝑋(𝑥 ) = 𝑐1 𝑒 √−𝜆𝑥 + 𝑐2 𝑒 −√−𝜆𝑥

Conditions de bord :
𝑋(0) = 0 ⇒ 𝑐1 + 𝑐2 = 0
{
𝑋 ′ (1) = 0 ⇒ 𝑐1 √−𝜆𝑒 √−𝜆 − 𝑐2 √−𝜆𝑒 −√−𝜆 = 0
Le déterminant principal :
1 1 1 1
| −√−𝜆
| = −√−𝜆 | √−𝜆 | = −√−𝜆 (𝑒 −√−𝜆 + 𝑒 √−𝜆 ) < 0 ∀𝜆 < 0
√−𝜆𝑒 √−𝜆 −√−𝜆𝑒 −𝑒 𝑒 −√−𝜆

La solution associée à cette condition est donc triviale.

• Si Δ < 0 ⇒ −4𝜆 < 0 ⇒ 𝜆 > 0, nous avons deux racines complexes conjuguées :
𝑟1,2 = ±𝑖√𝜆, 𝜆>0
La solution est :
𝑋(𝑥 ) = 𝑐1 cos √𝜆𝑥 + 𝑐2 sin √𝜆𝑥
Conditions de bord :
𝑋(0) = 0 ⇒ 𝑐1 = 0
{ ′
𝑋 (1) = 0 ⇒ 𝑐2 √𝜆 cos √𝜆 = 0
Pour la deuxième équation du système, la solution 𝑐2 = 0 sera rejetée car elle conduit à une
solution triviale pour 𝑋(𝑥 ). Par contre, √𝜆 cos √𝜆 = 0 permet d’écrire :
 18

𝜋 2𝑘 + 1 2
√𝜆 = 0 ⇒ 𝜆 = 0 à 𝑟𝑒𝑗 𝑐𝑎𝑟 𝜆 > 0 , 𝑜𝑢 √𝜆 = + 𝑘𝜋 ⇒ 𝜆𝑘 = ( 𝜋)
2 2
Ainsi la solution est :

2𝑘 + 1 2𝑘 + 1 2
{𝑋𝑘 (𝑥 ) = 𝑐2 sin √𝜆𝑘 𝑥 = 𝑐2 sin ( 𝜋𝑥) , 𝜆𝑘 = −𝜇𝑘 = ( 𝜋) ; 𝑥, 𝑐2 ≠ 0}
2 2

Le premier problème 𝑇 ′ (𝑡) − 𝜇𝑇(𝑡) = 0 a pour solution 𝑇(𝑡) = 𝑐𝑒 𝜇𝑡

2𝑘+1 2
Comme 𝜇𝑘 = − ( 𝜋) ;
2
2𝑘+1 2
−( 𝜋) 𝑡
𝑇𝑘 (𝑡) = 𝑐𝑒 2

Ainsi, on peut déterminer 𝑢(𝑡, 𝑥 ) ;

2
2𝑘 + 1 −(
2𝑘+1
𝜋) 𝑡
𝑣𝑘 (𝑡, 𝑥 ) = 𝐶 sin ( 𝜋𝑥) 𝑒 2
2
Par le principe de la superposition des solutions ;
∞
2𝑘+1 2 2𝑘 + 1
−( 𝜋) 𝑡
𝑣(𝑡, 𝑥 ) = ∑ 𝑏𝑘 𝑒 2 sin ( 𝜋𝑥) , 𝑡 > 0, 0<𝑥<1
2
𝑘=1

Où 𝑏𝑘 sont les coefficients de Fourier en sinus donnés par :

2 1 2𝑘 + 1
𝑏𝑘 = ∫ 𝑣 (0, 𝑥 ) sin ( 𝜋𝑥) 𝑑𝑥
1 0 2
Avec 𝑣(0, 𝑥 ) = −3 + sin(𝜋𝑥 ), on a :
1
2𝑘 + 1
𝑏𝑘 = 2 ∫ (−3 + sin(𝜋𝑥 )) sin ( 𝜋𝑥) 𝑑𝑥
0 2
1 1
2𝑘 + 1 2𝑘 + 1
= −6 ∫ sin ( 𝜋𝑥) 𝑑𝑥 + 2 ∫ sin(𝜋𝑥 ) sin ( 𝜋𝑥) 𝑑𝑥
0 2 0 2
1
1 2𝑘 + 1 2𝑘 + 1
= 0 + 2∫ [cos (𝜋𝑥 − 𝜋𝑥) − cos (𝜋𝑥 + 𝜋𝑥)] 𝑑𝑥
0 2 2 2
1 1
1 − 2𝑘 3 + 2𝑘
= ∫ cos ( 𝜋𝑥) 𝑑𝑥 − ∫ cos ( 𝜋𝑥) 𝑑𝑥
0 2 0 2
1 1
2 1 − 2𝑘 2 3 + 2𝑘
= [sin ( 𝜋𝑥)] − [sin ( 𝜋𝑥)]
(1 − 2𝑘)𝜋 2 0 (3 + 2𝑘)𝜋 2 0
2 1 − 2𝑘 2 3 + 2𝑘
= [sin ( 𝜋) − 0] − [sin ( 𝜋) − 0]
(1 − 2𝑘)𝜋 2 (3 + 2𝑘)𝜋 2
2 2𝑘 − 1 2 2𝑘 + 3
= sin ( 𝜋) − sin ( 𝜋)
(2𝑘 − 1)𝜋 2 (2𝑘 + 3)𝜋 2
2 2 2 1 1
= (−1)𝑘 − (−1)𝑘+1 = ( + ) (−1)𝑘
(2𝑘 − 1 )𝜋 (2𝑘 + 3)𝜋 𝜋 2𝑘 − 1 2𝑘 + 3
2 2𝑘 + 3 + 2𝑘 − 1 2 4𝑘 + 2 4 2𝑘 + 1
= (−1)𝑘 = =
𝜋 (2𝑘 − 1)(2𝑘 + 3) 𝜋 (2𝑘 − 1)(2𝑘 + 3) 𝜋 (2𝑘 − 1)(2𝑘 + 3)
 19

4 2𝑘 + 1
𝑏𝑘 =
𝜋 (2𝑘 − 1)(2𝑘 + 3)

D’où
∞ ∞
2𝑘+1 2 2𝑘 + 1 4 2𝑘 + 1 2𝑘+1 2𝑘 + 1
2
−( 𝜋) 𝑡 −( 𝜋) 𝑡
𝑣 (𝑡, 𝑥 ) = ∑ 𝑏𝑘 𝑒 2 sin ( 𝜋𝑥) = ∑ 2 𝑒 2 sin ( 𝜋𝑥)
2 𝜋 4𝑘 + 4𝑘 − 3 2
𝑘=1 𝑘=1

Ainsi, la solution 𝑢(𝑡, 𝑥 ) = 𝑣(𝑡, 𝑥 ) + 𝑞 (𝑥 ) est :

∞ 2
4 2𝑘 + 1 −(
2𝑘+1
𝜋) 𝑡 2𝑘 + 1
𝑢(𝑡, 𝑥 ) = 3 + sin(𝜋𝑥 ) + ∑ 2 𝑒 2 sin ( 𝜋𝑥)
𝜋 4𝑘 + 4𝑘 − 3 2
𝑘=1

3. Résoudre le problème aux valeurs initiales et des bords ci-dessous :

𝒖𝒕𝒕 = 𝟗𝒖𝒙𝒙 , 𝒕 > 𝟎, 𝟎<𝒙<𝟏
𝒖(𝒕, 𝟎) = 𝟎, 𝒖𝒙 (𝒕, 𝟏) = 𝟎, 𝒕 > 𝟎
𝒖(𝟎, 𝒙) = 𝟎, 𝟎≤𝒙≤𝟏
a) 𝟏
𝒙, 𝟎≤𝒙≤
𝟐
𝒖𝒕 (𝟎, 𝒙) = { 𝟏
{ 𝟏 − 𝒙, ≤𝒙≤𝟏
𝟐

Solution
Méthode de résolution : séparation des variables. Posons 𝑢(𝑡, 𝑥 ) = 𝑇(𝑡)𝑋(𝑥 )
⇒ 𝑢𝑥 = 𝑇(𝑡)𝑋 ′ (𝑥 ), 𝑢𝑡 = 𝑇 ′ (𝑡)𝑋(𝑥 ), 𝑢𝑥𝑥 = 𝑇(𝑡)𝑋 ′′ (𝑥 ), 𝑢𝑡𝑡 = 𝑇 ′′ (𝑡)𝑋(𝑥 )

En substituant dans le problème on a :

𝑇 ′′ (𝑡) 𝑋 ′′ (𝑥 )
𝑇 ′′ (𝑡)𝑋(𝑥 ) = 9𝑇(𝑡)𝑋 ′′ (𝑥 ) ⇒ = =𝜇 (1)
9𝑇(𝑡) 𝑋 (𝑥 )
Le système d’EDO (1) donne :
𝑋 ′′ (𝑥 )
=𝜇
𝑋 (𝑥 ) 𝑋 ′′ (𝑥 ) − 𝜇𝑋 (𝑥 ) = 0
′′ ( ) ⇒ { ′′
𝑇 𝑡 𝑇 (𝑡) − 9𝜇𝑇 (𝑡) = 0
=𝜇
{ 9𝑇(𝑡)
Puisque 𝑢(𝑡, 𝑥 ) ne doit pas être une solution triviale alors, les conditions de bord et initiales
permettent d’écrire :
𝑢(𝑡, 0) = 0 ⇒ 𝑇(𝑡)𝑋(0) = 0 ⇒ 𝑋(0) = 0, 𝑡 > 0, 0<𝑥<1
′ ′
𝑢𝑥 (𝑡, 1) = 0 ⇒ 𝑇(𝑡)𝑋 (1) = 0 ⇒ 𝑋 (1) = 0, 𝑡>0
( ) ( ) ( )
𝑢 0, 𝑥 = 0 ⇒ 𝑇 0 𝑋 𝑥 = 0 ⇒ 𝑇 0 = 0, ( ) 0≤𝑥≤1

On a donc deux problèmes qui sont :

𝑋 ′′ (𝑥 ) − 𝜇𝑋(𝑥 ) = 0 𝑇 ′′ (𝑡) − 9𝜇𝑇(𝑡) = 0

{ 𝑒𝑡, {
𝑋(0) = 0, 𝑋 ′ (1) = 0 𝑇(0) = 0
 20

𝑋 ′′ (𝑥 ) − 𝜇𝑋(𝑥 ) = 0
Le premier problème { est un problème de S-L. En posant 𝜆 = −𝜇 on a :
𝑋(0) = 0, 𝑋 ′ (1) = 0
𝑋 ′′ (𝑥 ) + 𝜆𝑋(𝑥 ) = 0 dont l’équation caractéristique associée est 𝑟 2 + 𝜆 = 0 et dont le discriminant
est Δ = −4𝜆.
On se ramène au cas de la question précédente pour laquelle la solution non triviale pour 𝑋(𝑥 ) est
trouvée pour Δ < 0 ⇒ 𝜆 > 0. Les 𝜆𝑘 sont telles que :
𝜋 2𝑘 + 1 2
√𝜆 = + 𝑘𝜋 ⇒ 𝜆𝑘 = ( 𝜋)
2 2
2𝑘+1 2𝑘+1 2
Et la solution est {𝑋𝑘 (𝑥 ) = 𝑐2 sin √𝜆𝑘 𝑥 = 𝑐2 sin ( 𝜋𝑥) , 𝜆𝑘 = −𝜇𝑘 = ( 𝜋) ; 𝑥, 𝑐2 ≠ 0}
2 2

𝑇 ′′ (𝑡) − 9𝜇𝑇(𝑡) = 0
Le deuxième problème { étant lié au premier alors, 𝜆 = −𝜇 > 0 ⇒ 𝜇 < 0. De
𝑇(0) = 0
ce fait, le polynôme caractéristique 𝑝(𝑟) = 𝑟 2 − 9𝜇 n’a que des racines complexes :
𝑟1,2 = ±3𝑖 √−𝜇

⇒ 𝑇(𝑡) = 𝑐1 cos 3√−𝜇𝑡 + 𝑐2 sin 3√−𝜇𝑡

On peut encore écrire 𝑇(𝑡) = 𝐴 cosh 3√−𝜇𝑡 + 𝐵 sinh 3√−𝜇𝑡
Tenant compte de la condition initiale 𝑇(0) = 0 ⇒ 𝐴 = 0 on a donc :
2𝑘 + 1
𝑇(𝑡) = 𝐵 sinh 3√−𝜇𝑡 ⇒ 𝑇𝑘 (𝑡) = 𝐵 sinh 3√−𝜇𝑘 𝑡 = 𝐵 sinh (3 𝜋𝑡)
2
Les solutions non triviales sont donc de la forme :
2𝑘 + 1 2𝑘 + 1
𝑢𝑘 (𝑡, 𝑥 ) = 𝐶 sin ( 𝜋𝑥) sinh (3 𝜋𝑡)
2 2
Par le principe de superposition la solution générale est la série :
∞
2𝑘 + 1 2𝑘 + 1
𝑢(𝑡, 𝑥 ) = ∑ 𝐵𝑘 sin ( 𝜋𝑥) sinh (3 𝜋𝑡)
2 2
𝑘=1

1
𝑥, 0 ≤ 𝑥 ≤
2
Cette solution satisfera la condition de bord 𝑢𝑡 (0, 𝑥 ) = { 1 si on a :
1 − 𝑥, ≤𝑥≤1
2
∞
2𝑘 + 1 2𝑘 + 1
𝑢(0, 𝑥 ) = ∑ 𝐵𝑘 sin ( 𝜋𝑥) sinh (3 𝜋𝑡)
2 2
𝑘=1
∞
2𝑘 + 1 2𝑘 + 1 2𝑘 + 1
⇒ 𝑢𝑡 (0, 𝑥 ) = ∑ (3 𝜋) 𝐵𝑘 sin ( 𝜋𝑥) cosh (3 𝜋𝑡) = 𝑓(𝑥 )
2 2 2
𝑘=1

Autrement dit si 𝑓(𝑥 ) peut se développer en demi série sinus de Fourier avec :
2𝑘 + 1 2𝑘 + 1 2 1 2𝑘 + 1
(3 𝜋) 𝐵𝑘 sin ( 𝜋. 1) = ∫ 𝑢𝑡 (0, 𝑥 ) cosh (3 𝜋𝑥) 𝑑𝑥
2 2 1 0 2
1
4 1 𝑘∫
2𝑘 + 1
⇒ 𝐵𝑘 = ( −1) 𝑢𝑡 (0, 𝑥 ) cosh (3 𝜋𝑥) 𝑑𝑥
3𝜋 2𝑘 + 1 0 2

Avec
 21

1
𝑥, 0≤𝑥≤
𝑢𝑡 (0, 𝑥 ) = { 2
1
1 − 𝑥, ≤𝑥≤1
2
∞
2𝑘 + 1 2𝑘 + 1
𝑢(𝑡, 𝑥 ) = ∑ 𝐵𝑘 sin ( 𝜋𝑥) sinh (3 𝜋𝑡)
2 2
𝑘=1