158 CORDINATEUR-BASEDNUMERICAL AND SSTATISTIQUET TECHNIQUES

En multipliant les racines de cette équation par 10, nous obtenons la nouvelle équation comme

f3(x) = 3x3+ 13080x2+ 28989600x– 18416000 = 0

Nous pouvons facilement voir que3(0) < 0 and f3(1) > 0, c'est-à-dire, la racine de3(x) = 0 ment
entre 0 et 1.
∴ La troisième chiffre après la virgule est zéro. Nous pouvons nous arrêter ici car le cas
requiert que la racine soit correcte à 2 décimales. Donc la racine est 1,32.

DEVOIR 3.8
1. Trouvez une racine des équations suivantes exactes à trois décimales en utilisant la méthode de Horner.
méthode.
(i)x3+ 3x2– 12x – 11 = 0 (ii)x4+x3- 4x2- 16 = 0
(iii)x3– 30 = 0.
2. Trouvez la racine positive de l'équation x3+x2x – 100 = 0, arrondi à quatre décimales
en utilisant la méthode de Horner.

3.29 Méthode de NEWTON-RAPHSON

Cette méthode est généralement utilisée pour améliorer le résultat obtenu par l'un des
previous methods. Letx0sois une racine approximative de f(x) = 0 et laisse x1=x0+hbe
la racine correcte afin que f(x1) = 0.
Développement de f(x)0+h) par la série de Taylor, nous obtenons

h2
f(x0) +hf′(x0) + f″ (x0) + ...... = 0
2!
Depuis qu'il est petit, il est négligé.2et des puissances supérieures de h, nous obtenons

f(x0)
f(x0) +hf′(x0) = 0 ouh= – (29)
f′ (x 0)
Une meilleure approximation merci0est donc donné parx1, où

f(x0)
x1=x0–
f′ (x 0)
Les approximations successives sont données par x2,x3, ......., xn +1, où
f(xn )
xn+ 1=xn– (30) (n= 0, 1, .......)
f′ (x n )
quelle est la formule de Newton-Raphson.
 UnALGÉBRIQUE ETTTRANSCENDENTALEQUATIONS 159

REMARQUE Cette méthode est utile dans les cas de grandes valeurs de f.′(x), c'est-à-dire, quand le
Le graphe de f(x) qui croise l'axe des x est presque vertical.
2.Si f′ (x) est zéro ou presque zéro, la méthode échoue.
3. La formule de Newton converge à condition que l'approximation initiale x0est
choisi suffisamment près de la racine.
Au début, nous devinons deux nombres b et c tels que f(b) et f(c)
sont de signes opposés. Puis la première racine approximative a se trouve entre b
et c.
Cette méthode est également utilisée pour obtenir des racines complexes.

3.30 CONVERGENCE

Comparer (30) avec xn+ 1 = φ (xn) de la méthode d'itération, nous obtenons

f(xn )
φ(xn) =xn +1=xn–
f′ (x n )

f(x)
En général, φ(x) = x -
f′ (x n )

f(x)f″ (x)
qui donne φ′(x) =
[ f ′ (x)] 2
Puisque la méthode d'itération converge si |φ′ (x) | < 1
∴ La méthode de Newton converge si
|f(x)f″ (x) | < [f′ (x)]2
dans l'intervalle considéré.
Supposonsf(x),f′(x), etf″(x) pour être continu, nous pouvons sélectionner un petit intervalle
dans les environs de la racineα dans lequel la condition ci-dessus est satisfaite.
Le taux auquel la méthode d'itération converge si l'approximation initiale
à la racine est suffisamment proche de la racine désirée est appelé le taux de
convergence.

3.31 ORDRE DE CONVERGENCE

Supposexndiffère de la racineα par une petite quantiténafin que

xn= α +enetxn+ 1= α +en+ 1
160 CORDINATEUR-BASEDNNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

f(α + en )
Alors (30) devient, en +1=en–
f′ (α + en )
en 2
f(α) + enf′ (α) + f″ (α) + ......
2!
=en– (Par l'expansion de Taylor)
f′ (α) + enf″ (α) + ......

en 2
enf′ (α) + f″ (α) + .......
=en– 2 | ∵ f(α) = 0
f′ (α) + enf″ (α) + ......
2
enf″ (α)
= Négligez les puissances élevées de en
2[f′ (α) + enf″ (α)]

en 2 f″ (α)
=
SR
2 f″ (α) Vous
f′ (α)
1+ e
T f′ (α) W V n

f″ (α) R f″ (α) U
2
e −1
=
n
.
2f′ (α) T
S 1+ e V
f′ (α) W
n

f″ (α) R U
2
e f″ (α)
2 f′ (α) S1 − e f′ (α) + .......V
n
=
T W
n

e 3Es-tu
2
en 2 f″ (α) f″ (α)
=
2f′ (α)
− n
2 f′ (α)
S
T V+ .......
W
en+1 1f″ (α) en Es-tu
f″ (α)
2

TSf′ (α) W
V+ .......
ou = −
en 2 2f′ (α) 2
f ″ (α)
≈ (Négligeant les termes contenant des puissances de e)n)
2f′ (α)
D'où, par définition, l'ordre de convergence de la méthode de Newton-Raphson
c'est 2, c'est-à-dire que la méthode de Newton-Raphson est de convergence quadratique.

Cela montre également que l'erreur subséquente à chaque étape est proportionnelle à la
square of the previous error and as such theconvergence is quadratic.
Ainsi, si à la première itération nous avons une réponse correcte à un décimal
lieu, alors il devrait être correct à deux décimales lors de la deuxième itération, et à
quatre places à la troisième itération.
Cela signifie que le nombre de décimales correctes à chaque itération est
presque doublé.
 UnALGÉBRIQUE ETTTRANSCENDENTALEQUATIONS 161

∴ La méthode converge très rapidement.

En raison de sa convergence quadratique, la formule (30) est également appelée asecond.
formule de commande.

3.32 GEOMETRICAL INTERPRETATION

La tex0soyez un point près de la racineα de l'équation f(x) = 0, puis la tangente en A{x0,f(x0)}
est
y–f(x0) =f′(x0) (x–x0)

Un{x0, f(x0 )}
f(x)
=
y

A1
A2

O x2x1 x0 X

a
f(x0)
Il coupe l'axe des x à x1=x0 –
f′ x 0)
quelle est une approximation de la racineα. Si A1correspond à tox1sur la courbe, puis
la tangente en A1coupera l'axe des x à x2, plus proche deα et est donc un autre
approximation à la racineα.
Répétant ce processus, nous approchons de la racineα assez rapidement. D'où le
la méthode consiste à remplacer la partie de la courbe entre A et l'axe des x par
la pente de la tangente à la courbe en A0.

3.33 ALGORITHME DE LA MÉTHODE DE NEWTON-RAPHSON

Étape 01. Début du programme

Étape 02. Entrez les variables x0, n pour la tâche
Étape 03. Epsilon et delta
Étape 04. pour i = 1 et répéter si i <= n
Étape 05. f0 = f(x0)
162 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

Étape 06. df0 = df(x1)

Étape 07. si |df0| <= delta
a. Imprimer la pente trop petite

b. Imprimez x0, f0, df0, i

c. Fin du programme
Étape 08. x1 = x0 - (f0 / df0)
Étape 09. si |(x1-x0)/x1| < epsilon
a. Imprimer convergent
b. Imprimer x1, f(x1), i
c. Fin du programme
Étape 10. x0 = x1
Étape 11. Fin de boucle

3.34 DIAGRAMME DE FLUX DE LA MÉTHODE DE NEWTON-RAPHSON

DÉMARRER

Définir la fonction f(x)

Définir la fonction df(x)

Obtenez les valeurs de x0,aerr,maxitr

Boucler pour itr= 1 à maxitr

h = f (x0) / df(x0)
x1= x0– h

Imprimer itr, x1

Isfabs(h) Non x0= x1

< aerr
Fin de boucle(itr)
Oui
Imprimer la solution
solution fait
ne converge pas
ARRÊTER
ARRÊTER
 AALGÉBRIQUE ETTTRANSCENDENTALEQUOTATIONS 163

3.35 FORMULES ITÉRATIVES DE NEWTON POUR TROUVER L'INVERSE,

RACINE CARRÉE
1. Inverse. Le réciproque ou l'inverse d'un nombre 'a' peut être considéré comme un
1
racine de l'équation -a= 0, ce qui peut être résolu par la méthode de Newton.
x

1 1
Depuis f(x) = –a,f′(x) = –
x x2
∴ La formule de Newton donne

F
G
1
− un I
xn + 1 =xn+ H xJ K
n

F1 I
G
Hx JKn
2

xn+ 1=xn(2 – axn)

2. Racine carrée. La racine carrée de 'a' peut être considérée comme une racine de l'équation
x2–a= 0, résoluble par la méthode de Newton.
Depuis f(x) = x2-a,f′(x) = 2x

xn−2 a
xn +1=xn–
2x n

1Fx + a I
2G JK
xn+ 1 =
H
n
xn

1
3. Racine carrée inverse. L'équation est –a= 0
x2
La formule itérative est

1
2 xn(3 –a xn)
xn+ 1 = 2

4. Formule générale pour pthracine.LaPthla racine d'un peut être considérée comme un
racine de l'équation xp– a= 0. Pour résoudre cela par la méthode de Newton, nous avons

f(x) = xp–et donc, f′(x) =pxp– 1

164 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

(x n−p a)
∴ La formule itérative est xn+ 1=xn– p− 1
pxn

( p − 1)xn+ ap
xn+ 1 = p− 1
pxn

Aussi, la formule générale pour l'inverse de pthracine de bois

Fp+ 1 − essieuJe. pn
x n+ 1 =x
n G
H p JK
3,36 TAUX DE CONVERGENCE DE LA RACINE CARRE DE NEWTON POUR-
MULA

Laisseza = α tel un a=α2. Si nous écrivons

F1+ e I
G
H 1− eJ K
n
x=α
n
n

F1 + e I n+ 1
ensuite, xn+ 1 =α
G
H1 − e JK n+ 1
(31)

= F
1 a I
Aussi, par formule, xn+ 1
G
2H
x+n
x JK n
, nous obtenons

1 LF 1 + e I a1− e F O I
= M
NH 1 − eJKα 1 +Ge H
G JP
n n
xn+ 1 α +
2 M Q K
Pn n

F1+ e I n
2
= α
G
H1 − e JK n
2 (32) ( ∵ un=α2)

En comparant (31) et (32), nous obtenons en+ 1=en 2

confirmation de la convergence quadratique de la méthode de Newton.

3.37 TAUX DE CONVERGENCE DE LA FORMULE INVERSÉE DE NEWTON

1 1
Laissez α= c'est-à-dire, a = . Si nous écrivons x n= α(1–en)
a α
alors, xn+ 1= α (1 –en+ 1 )
 UnALGÉBRIQUE ETTTRANSCENDENTALEQUATIONS 165

Par formule, xn+ 1=xn(2 –axn), nous obtenons

xn+ 1= α(1–en) [2 –aα (1–en)] = α(1–en) 2 | ∵ aα = 1

En comparant, nous obtenons en+1 =en, 2 la convergence est quadratique.

donc,

EXEMPLES
Exemple 1. En utilisant la méthode de Newton-Raphson, trouvez la racine réelle de la
Équation 3x = cos x + 1 corrigée à quatre décimales. Donnez un programme informatique.
utiliser 'C'.

Sol.Let f(x) = 3x - cos(x) - 1

Depuis f(0) = -2 = (-)ve;
f(1) = 1.4597 = (+)ve
∴ Une racine de f(x) = 0 se situe entre 0 et 1. Elle est plus proche de 1. Prenons
x0= 0,6.
Aussi, f′(x) = 3 + sinx
La formule d'itération de Newton donne,

f(xn )
xn+ 1=xn–
f′ (x n )

3x n− cosxn− 1 xnsinxn+ cosxn+ 1

=xn– =
3+ sinxn 3+ sinxn

Si n = 0, la première approximation x1est donné par,

x0sinx0+ cosx0+ 1
x1 =
3+ sinx0

0,6 sin 6+ cos 0.6+ 1

= = .6071
3+ sin 0,6
Si n = 1,
la deuxième approximation est
x1sinx1+ cosx1+ 1
x2=
3+ sinx1

.6071sin (.6071)+ cos(.6071)+ 1

= = 0.6071
3+ sin (0,6071)
Clair1=x2. Donc la racine désirée est 0.6071, corrigée à 4 décimales.
166 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

/* ********************************************************
Programme réalisé pour NEWTON RAPHSON pour résoudre l'équation
******************************************************* *\

//....y compris les fichiers d'en-tête source

# inclure <stdio.h>
# inclure <conio.h>
# inclure <math.h>
# inclure <process.h>
# inclure <string.h>

//....définition des formules

# définir f(x) 3*x - cos(x) - 1
# définir df(x)3+sin(x)
//...Déclaration de prototype de fonction
void NOUVEAU_RAP();

//... Fil d'exécution principal

vide principal()

{
effacer();
Solution par la méthode NEWTON RAPHSON
L'équation est :
3*X - COS X - 1=0
NOUVEAU_RAP();
getch();
}
//...Déclaration de fonction
vide NOUVEAU_RAP()
{
//...Champ de Déclaration Interne
long float x1,x0;
long float f0,f1;
float long df0;
int i=1;
int itr;
 AALGÉBRIQUE ETTRANSCENDENTALECITATIONS 167

flottant EPS;
erreur flottante;
/*Trouver une RACINE Approximative de l'Équation Donnée, Ayant
+ve Valeur*/
pour(x1=0;;x1 +=0.01)
{
f1=f(x1);
si (f1 > 0)
{
casser;
}
}
/*Trouver une racine approximative de l'équation donnée, ayant
-ve valeur*/
x0=x1-0,01;
f0=f(x0);
printf(" Entrez le nombre d'itérations : ");
scanf(" %d",&itr);
printf(" Entrez l'erreur maximale possible : ");
scanf("%f", &EPS);
si (fabs(f0) > f1)
{
La racine est proche de %.4f
}
Si (f1 > fabs(f(x0)))
{
\n\t\t La racine est proche de %.4f\n
}
x0=(x0+x1)/2;
pour(;i<=itr;i++)
{
f0=f(x0);
df0=df(x0);
x1=x0 - (f0/df0);
La %d approximation de la racine est :
%f",i,x1);
168 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

error=fabs(x1-x0);
si (erreur < EPS)
{
pause;
}
x0 = x1;
}
si(erreur > EPS)
{
REMARQUE:-
Le nombre d'itérations n'est pas suffisant.
}
------------------------------
The root is %.4f ",x1);
------------------------------
}

SORTIE
Solution par la méthode de NEWTON RAPHSON
L'équation est :
3*X - cos X - 1=0
Enter the number of iterations: 10
Entrez l'erreur maximale possible : 0,0000001
La racine est proche de 0,6100
L'approximation 1 de la racine est : 0,607102
L' approximation à la racine est : 0,607102
L'approximation à la racine est : 0,607102
--------------------------------
La racine est 0,6071
--------------------------------
Exemple 2. En utilisant la méthode itérative de Newton, trouvez la racine réelle de x log10x = 1.2
corriger à cinq décimales.
Sol. f(x) = xlog10x - 1,2

∵ f(1) = – 1.2 = (–)ve

f(3) = 3 log103 – 1,2 = (+)ve
 AALGÉBRIQUE ETTTRANSCENDANTALEQUATIONS 169

Donc, une racine de f(x) = 0 se situe entre 1 et 3.

Prenons x0 = 2.
et f ′(x) = log10x + log10e = log10x + 0,43429
La formule d'itération de Newton donne,

f(xn )
xn +1=xn–
f′ (x n )

x njournal 10xn− 1.2 .43429xn+ 1.2

=xn– = (33)
journal 10xn+ .43429 journal10xn+ .43429

Étant donné n = 0, la première approximation est

.43429x0+ 1.2
x1= = 2,81 (∵ x0= 2)
journal10 2+ .43429
De même, en prenant n = 1, 2, 3, 4 dans (33), nous obtenons

x2= 2.741,x3= 2.74064,x4= 2.74065, x5= 2.74065

Clairement, x4=x5
D'où la racine requise est 2.74065, arrondie à cinq décimales.

Exemple 3.Évaluer12 à quatre décimales par la méthode itérative de Newton.

Sol.Let x= 12 donc que x2- 12 = 0 (34)

Takef(x) = x2– 12, la formule d'itération de Newton donne,

f(xn ) xn−2 12 1Fx + 12 I

xn+ 1=xn–
f ′ (x n )
=xn –
2x n
= G
2H
n
xn JK (35)

Maintenant, puisque f(3) = - 3 (-)ve

f(4) = 4 (+)ve
∴ La racine de (34) se situe entre 3 et 4.
Donnéx0= 3,5, (35) donne,

1Fx + 12 I= 1 F 3,5+ Je= 3.4643

12
2 H x JK
G 2G JK
x1= 0
H 0 3.5

1F 12 I
x2= G
2H x
x+ J = 3.4641
1
K 1

x3= 3,4641
170 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

Depuis2=x3jusqu'à 4 décimales
nous avons 12 = 3.4641.
Exemple 4. En utilisant la méthode itérative de Newton, trouvez la racine réelle de

x sin x + cos x = 0 qui est près de x =π, corrigez à 3 décimales.

Sol. Nous avons
f(x) = xsinx + cosx et f′(x) = xcosx
La formule d'itération est
xnsinxn+ cosxn
xn +1=xn–
x ncosxn

x0sinx0+ cosx0 π péchéπ + cosπ

avecx0= π, x1=x0– =π– = 2.8233
x0cosx0 π cosπ
Les itératives successives sont

x2= 2.7986, x3= 2.7984, x4= 2,7984

Depuis3=x4, la racine requise est 2,798, arrondie à trois décimales.
Exemple 5. Trouvez une racine réelle de l'équation x = e–x en utilisant la méthode de Newton-Raphson
méthode.
Sol.Nous avons f(x) = xex-1
alors, f′(x) = (1+x)ex
Laissez x0= 1

alors, F
G
e− 1 SI
JG
1
1 + JeJ = 0.6839397
1
x=
1 1–
H K H eK
2e = 2
Maintenant, f(x1) = 0.3553424 et f′(x1) = 3.337012
0,3553424
afin que x2= 0.6839397 – = 0.5774545
3,337012
En procédant de cette manière, nous obtenons

x3= 0.5672297, x4= 0.5671433

Ainsi, la racine requise est 0,5671, arrondie à 4 décimales.
Exemple 6. Trouvez, jusqu'à quatre décimales, la plus petite racine de l'équation.
e–x=sin x.
Sol. L'équation donnée est
f(x)≡ e–x– sinx = 0
 AALGÉBRIQUE ETTTRANSCENDENTALEQUATIONS 171

e − x−
n
sinxn
afin que, xn +1=xn+
e − x+
n
cosxn
Takex0= .6 alors,
x1= .58848, x2= .588559
Ainsi, la valeur souhaitée de la racine est 0,5885.
Exemple 7.(i) Trouvez une valeur positive de (17) 1/3, corrigé à quatre décimales,
par la méthode de Newton-Raphson.
(ii) Trouvez la racine cubique de 10.
Sol.(i) La formule itérative est

1F 2x+ a Je
xn+ 1 =
G
3H xJ K
n
2 n
(36)

Ici a = 17
38= 2 et327= 3
Takex0= 2,5 ∵
En mettant n = 0 dans (36), nous obtenons

1F2x+ 17= Je 1 F
G 5+ Je J = 2.5733
17
x1=
G
3H J
x K H
0
3 6,25
2 0 K
En plaçant n = 1 dans (36), nous obtenons

1F 2x+ 17=Je 1 F
G
17 Je
G J 6,6220 J
x2=
x K 3H K= 2,5713
5.1466+
3H
1
2 1

En remettant n = 2 dans (36), nous obtenons

1F 2x+ 17=Je1 F
G
17 Je
G J 6.61158 J
x3=
x K3 H K= 2,57128
5,1426+
3H
2
2 2

En mettant n = 3 dans (36), nous obtenons

1F 2x+ 17= Je1 F

G
17 Je
3G 6.61148 J
H xJ K 3 H
x4= 3 5.14256+
2 3 K= 2,57128
Depuis x3etx4accord à quatre décimales, la racine requise est 2,5713
arrondir à quatre décimales.

2x n+3 un 1F 2x+ a Je
(ii) xn+ 1 =
3x n 2 =
G
3H xJ K
n
2 n
172 COrdinateur-BASEDNNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

38 327
Prendre x0= 2,5 (∵ = 2 et = 3)
∴ x1 = 2,2 (n =0)
x2= 2,155 (n =1)
x3= 2.15466 (n= 2)

∴ 310 ≈ 2.15466
Exemple 8. Montrer que les deux suites suivantes ont toutes deux une convergence de
le second ordre avec la même limite un.

1 F un I et, x = 1 x F3− . x I n
2
xn + 1 =
2 n G
x 1+
H xn 2 JK 2 GJ a
H n+1
K n

1 F1 + a I
Sol. Depuis, xn+ 1 = x
2 n G
Hx n
JK
nous avons
2

1 F1 + a I – a = 1 Fx+ a − 2un Je
xn+ 1 – a = x
2 n G
Hx n
2JK 2GH x JK n
n

1Fx − unJe 1 2
=
G
2H
n
x J
K n
=
2x
(x – un)
n
n
2

1
Ainsi, en+ 1 = e 2 (37)
2x n n
ce qui montre la convergence quadratique. De même pour le second,

F x I – un
1 n
2
xn +1 – a =
G
H a JK
x 3−
2 n

= x F I 2
1 x n
2 G
H J
1 −un
n + (x –
K n
a)

xn
L x O
N e x+ u n j P
M
n
= (a–xn) +2 (xn– un) = (x–n un)
2a 1−
2a Q n

xn− a
en+ 1 = [2a–xn2–xn a]
2a
xn− un
= [(a–xn) +2 (a–xn)] un
2a
 AALGÉBRIQUE ETTTRANSCENDENTALEQUATIONS 173

Fx − a
n I (x–
=– G
H 2a JKn a ) (x+
n 2a )

(x n− a) 2 (x n+ 2a)
en +1= – (xn+ 2a )=– .en 2 (38)
2a 2a
qui montre la convergence quadratique.

Exemple 9. Si xnest une approximation suffisamment proche de a , montrez que l'erreur

dans la formule
1 F Iparle de rd1 que dans la formule,
un
xn+ 1 =
2 G
xn 1+
H x n
JK 3
2

1 F x ,Iet déduire que la formule

n
2
xn+ 1 =
2 GH un JK
x 3−
n

x F 3a − x Je n2

8 G
H 6 x+ un JKdonne une séquence avec convergence d'ordre trois.
n
xn+ 1 = 2
n

Sol.Sincexnest très proche de a

Fx + 2x I e
n n
en+ 1 ~
− – G
H2x JK n
2 n
2 De (38)

1
=3. e 2 (39)
2x n n

Une simple observation montre qu'à partir de (37) (voir Ex. 8) et (39), l'erreur
1
dans la première formule avantn +1parle de rd de cela dans la deuxième formule.
3
Pour trouver le taux de convergence de la formule donnée, nous avons

xn F 6+3a xn 2 Je– 2
xn(6xnun+ 3a2− xn) 4
xn + 1 – a =
8 G
Hx n
2
−
unJ
K
a =
8axn 2 – un

2
6xnun+ 3a2− xn− 8x4na a − (x n+ 3a) (xn− a) 3
= =
8xn un 8x n un

Fx + 3a Jee
n
G
H 8x a JK
∴ en +1 3
=– n
n

Il montre que la formule ci-dessus a une convergence de troisième ordre.

174 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

Exemple 10. Appliquez la formule de Newton pour trouver les valeurs de (30)1/5 .
Sol. Pour trouver le pthracine ofa, nous avons

( p − 1)xn+ ap
xn+ 1 =
p x n p− 1
Ici, a= 30, p= 5, la première approximation est

4x0+5 30
x1= 4
5x0

Prendre x0= 1,9 nous obtenons1 = 1,98

Encore une fois, x2= 1,973
∴ Valeur = 1,973 (arrondi à 3 décimales).
Exemple 11. En utilisant la valeur de départ 2(1 + i), résoudre x4-5x3+ 20x2-40x + 60
= 0 par la méthode de Newton-Raphson, étant donné que toutes les racines de l'équation donnée
sont complexes.
Sol.Laissez f(x) = x4- 5x3+ 20x2- 40x + 60
afin que, f′(x) = 4x3- 15x2+ 40x –40
∴ La méthode de Newton-Raphson donne,

f(xn )
xn +1=xn–
f′ (x n )

xn−4 5x n + 20x
3 2
n− 40x n + 60
=xn –
4x n−3 15xn+ 40x
2
n− 40

3x n−4 10x n+ 20x

3
n− 60
2
=
4x n−3 15xn+ 40x
2
n− 40

Putn= 0, takex0= 2(1 + i) par essai, nous obtenons

x1= 1,92 (1 + i)
Encore, x2= 1,915 + 1,908i
Puisque les racines imaginaires se produisent en paires conjuguées, les racines sont 1,915 ± 1,908iup
à 3 décimales. En supposant que les autres paires de racines soientα ±iβ, alors

Fα + jeβ + α − jeβ Je
Sum = G
G + 1.915+ 1.908je J = 2α + 3,83 = 5
J
H + 1.915 − 1.908jeK
⇒ α 0,585
 UnALGÉBRIQUE ETTTRANSCENDANTALECITATIONS 175

De plus, le produit des racines = (α2+ β2) [(1.915)2+ (1.908)2] = 60

⇒ β = 2,805
Par conséquent, les deux autres racines sont 0,585 ± 2,805i.

Exemple 12. Obtenir la formule étendue de Newton-Raphson

f(x0) 1 {f(x0)}2. f″ (x0)

x1= x0– – .
f′ (x0) 2 {f ′ (x0)}3
pour la racine de l'équation f(x) = 0, également connue sous le nom de formule de Chebyshev de troisième ordre
commande.
Sol. Développer f(x) par la série de Taylor dans le voisinage de x0, nous obtenons
f(x) = 0⇒ f(x0) + (x–x0)f′(x0) = 0
f(x0)
⇒ x=x0–
f′ (x 0)
Ceci est une approximation à la racine.
f(x0)
∴ x1=x0–
f′ (x 0)
Encore selon la série de Taylor, nous avons

(x− x0) 2
f(x) = f(x0) + (x–x0)f′(x0) + f″(x0)
2
(x 1− x0 ) 2
∴ f(x1) =f(x0) + (x1–x0)f′(x0) + f″(x0)
2
Maisf(x1) = 0 asx1est une approximation de la racine.
1
∴ f(x0) + (x1–x0)f′(x0) + (x1–x0)2f″(x0) = 0
2

1 { f ( x 0 )}2f″ (x 0 )
ou f(x0) + (x1–x0)f′(x0) + =0
2 f′ (x 0) 2

f(x0) 1 { f ( x 0 )}2f″ (x 0 )
⇒ x1=x0– –
f′ (x 0) 2 { f ′ (x 0)}3
Cette formule peut être utilisée de manière itérative.

Example 13.The graph of y = 2 sin x and y = log x + c touch each other in the
le voisinage du point x = 8. Trouvez c et les coordonnées du point de contact.
Les graphiques se toucheront si les valeurs de dy/dx à leur point de
le contact est le même.
176 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

dy
Pour y = 2 sin x, = 2 cosx
dx
dy 1
Pour y = logx + c =
dx x
1
∴ 2 cosx = ⇒ xcosx – .5 = 0
x
Laissez f(x) = xcosx - 0.5
∴ f′(x) = cosx - xsinx
∴ La formule itérative de Newton est

xncosxn− 0,5
xn+ 1 =xn–
cosxn− xnsinxn

Forn= 0, x0= 8, première appli. x1= 7.793

Deuxième approximation, x2= 7.789≈ 7,79
Maintenant, y= 2 sin 7.79 = 1.9960
∴ Point de contact→ (7.79, 1.996)
Maintenant, y = logx + c

⇒ 1.996 = log 7.79 + c⇒ c = -0,054.

Exemple 14. En utilisant la valeur de départ x0= i, trouver un zéro de

x4+ x3+ 5x2+ 4x + 4 = 0.

Sol. Par la méthode de Newton

f(i) 3i
x1=i– =i– = 0,486 + 0,919i
f′ (i) 1 + 6i
f(.486+ .919i)
Maintenant, x2= .486 + .919i–
f′ (.486+ .919i)

F
G
− .292+ .174i Je
H1.78+ 6.005je JK
= .486 + .919i– = – 0,499 + 0,866i

− 1 + i3
La racine réelle est x = .
2
Exemple 15. Montrer que la racine carrée de N = AB est donnée par
S N
− 4 + S , où S = A + B.
N ~
 AALGÉBRIQUE ETTTRANSENDANTALEQUATIONS 177

Sol.Let x= N
⇒ x2– N = 0
Laissez f(x) = x2– N
∴ f′(x) = 2x
Par la formule de Newton-Raphson,

f(xn) x2n− Nxn N

xn+ 1=xn– = x n− = +
f′ (x n ) 2x n 2 2x n

A+ B
Laissez xn=
2
A+ B N S N
alors, xn+ 1 = + ~
− + Puisque S = A + B
4 Un+ B 4 S
Exemple 16. Déterminez les valeurs de p et q afin que le taux de convergence de
la méthode itérative
N
xn + 1= pxn+ q
xn 2
pour calculer N1/3devient aussi élevé que possible.
Sol.Nous avons x3= N
∴ f(x) = x3– N
Locationα soyez la racine exacte, nous avons
α3= N
N , nous obtenons
Substituer xn= α +en,xn+ 1= α +en+ 1, N = α3inxn+ 1=pxn+q
x2n

α3
α +en+ 1=p(α +en) +q
(α + en ) 2

α3
= p(α + en ) + q
F e I
2

G
H α JK
α 21 + n

F1 + e I
−2

= p(α + e) + qα HG α KJ
n
n

R2 e + 3 Fe I− ......... |U
= p(α + e)+ qα 1 −|S
2

|T α G Hα JK V
n n

|W
n
178 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

e2n
= p ( α + en ) + qα − 2qen+ 3q – ........
α
⇒ en+ 1 = (p+q– 1)α + (p– 2q)en+ 0(en) + ....... 2

Maintenant pour la méthode afin de devenir de l'ordre le plus élevé possible, c'est-à-dire, de l'ordre
2, nous devons avoir

p+q= 1 etp – 2q= 0

2 1
pour que, p= andq= .
3 3
Exemple 17. Comment la constante devrait-elleα être choisi pour garantir la rapidité maximale possible
convergence avec la formule d'itération?

αxn+ xn+ −2
1
xn+1 = .
α+1
Sol.Depuis limxn= limxn+1 = ξ, nous avons
n→ ∞ n→ ∞

F αξ + 1+ 1Je
G ξ JJ 2

ξ=G
H α+1 K
⇒ (α + 1)ξ3= αξ3+ ξ2 + 1
⇒ ξ3– ξ2– 1 = 0
ξ peut être obtenu en trouvant une racine de l'équation x3–x2– 1 = 0.
Nous avons f(x) = x3–x2– 1
f ′(x) = 3x2– 2x
Puisque f(1.45) = (–)ve et f(1.47) = (+)ve
∴ La racine se situe entre 1,45 et 1,47.
Laissez x0= 1,46
Par la méthode de Newton-Raphson,

Première approximation est

f(x0) Fx − x − 1 Je= 1.465601.

0
3
0
2
x1=x0–
f′ (x 0) = x0 – H3x − 2x JK
G 0
2
0

La deuxième approximation est

f(x1 ) Fx − x − 1 Je
1
3
1
2
x2=x1–
f′ (x 1 ) = x 1– GH3x − 2x JK= 1.46557
1
2
1
 UnLGEBRAIQUE ETTTRANSCENDANTALEQUATIONS 179

D'oùξ = 1,465 arrondi à trois décimales.

Maintenant, nous avons

αxn+ xn+ −2
1
xn+1 = (40)
α+1
Mettren= ξ +enetxn+1= ξ +en+1dans (40), nous obtenons
1
(α + 1)( ξ +en+ 1) = α(ξ +en) + +1
(ξ + en) 2

1 F
1+ I + 1
G
e −2

H ξ JK
n
= α(ξ +en) +
2
ξ
ce qui donne,

2 F I e+ O(e)
(1 +α)en+1 = α − 3
ξ G
H JK
n n
2

2
Pour une convergence rapide, nous devons avoirα =
ξ3

2
∴ α= = 0,636.
(1,465)3
Exemple 18. Méthode de Newton-Raphson pour résoudre l'équation f(x) = c, où
c est une constante à valeur réelle, appliquée à la fonction

SR quand|x|≤ 1U
cos x
f(x) =
Tcos x+ (x− 1) ,quand|x|≥ 1 V
2
W
2

Pour quel c est xn = (– 1) n, quand x0= 1 et les calculs sont effectués

sans erreurs ? Même dans l'arithmétique de haute précision, la convergence est problématique.
Expliquer.
Sol. f(x) – c = 0 (41)
En appliquant la méthode de Newton-Raphson à l'équation (41), nous obtenons

ML
f(x ) − c O

Nf′ (x ) PQ
n
xn+1 = x–n
n

Forn= 0, nous avons

ML
f(x ) − c O

Nf′ (x ) PQ
0
x1= x–0
0
180 CORDINATEUR-BASEDNUMÉRIQUE ETSSTATISTIQUET TECHNIQUES

=1–
L
M
cos1− c OP |∵x0= 1
N− sin 1 Q
Lcos1− c
– 1 = 1 +M sin 1
O
P |∵x1= (– 1)1= – 1
N Q
D'où – 2 sin 1 = cos 1 – c
⇒ c = cos 1 + 2 sin 1
avec cette valeur, nous obtenons

x2= 1, x3= – 1, ...,xn= (– 1)n

Depuisf′(x) = 0 entrex0et les racines et aussi à x=0, la convergence
est problématique malgré l'arithmétique de haute précision.

DEVOIR 3.9
1. En utilisant la méthode de Newton-Raphson, trouvez la racine de x 4–x– 10 = 0 qui est près de tox= 2,
corriger à trois décimales.
2.Calculer une racine positive de 2x - 10x = 7 par la méthode de Newton-Raphson, correct à
log
quatre décimales.
3. (i) Utilisez la méthode de Newton-Raphson pour trouver une racine de l'équation x3– 2x – 5 = 0.
(ii) Utilisez la méthode de Newton-Raphson pour trouver une racine de l'équation x3– 3x – 5 = 0.
4. Trouvez la racine réelle des équations
(i) logx = cosx (ii)x2+ 4 sinx = 0
par la méthode de Newton-Raphson, correct à trois décimales.
5. Utilisez la méthode de Newton-Raphson pour obtenir une racine correcte à trois décimales de la
les équations suivantes :

(i) sinx = 1 - x (ii)x3– 5x + 3 = 0 (iii)x4+x2– 80 = 0

(iv)x3+ 3x2- 3 = 0 (v) 4(x– sinx) = 1 (vi)x– cosx= 0
x
(vii) sinx= (viii)x + logx = 2 (ix) tanx=x.
2
6. Expliquez la méthode de Newton-Raphson pour calculer les racines. Appliquez-la pour trouver x à partir de
x2– 25 = 0. Écrivez un programme en utilisant 'C'.
7. Écrivez un programme informatique en 'C' pour trouver une racine réelle de l'équation f(x) = 0 par la méthode de Newton.
Méthode de Raphson.

8. En utilisant la méthode de Newton-Raphson, obtenez la formuleN pour

et trouver 20 corriger à
2 décimales.
9. Obtenez la racine cubique de 120 en utilisant la méthode de Newton-Raphson, en commençant par x 0= 4.5.
10.Développer un algorithme en utilisant la méthode de Newton-Raphson pour trouver la quatrième racine d'un

nombre positif N, et trouver 4 32 .