AYUDE Yen vedi eiii a
peu moin eti
et 8^, ..., r* par les éléments de premiere espéce de E". On voit que dans les matrices E,
seuls les éléments de premiére espéce sont nécessaires.
26. Variante : calcul par accroissements.
A partir du deuxiéme stade d'itération, on peut calculer les accroissements des diverses
quantités, au lieu de leurs valeurs mémes; cependant, ce procédé ne doit pas s'étendre au
calcul des matrices D, qui sont obtenues par multiplication; l'opération serait beaucoup plus
longue. On calculera donc, au deuxiéme stade : D' comme au n* 24, puis
AD = D! — [Dv
dont les éléments de seconde espéce fournissent
AM; = ADS, =D, ;
AM — AD —AD, =DL Di ;
AM = ADS, FAD), =D + Dh ,
AA = ADS, Di,
ensuite AC! par les équations linéaires (27), et, au lieu de (36), les éléments de première
espèce de
AE! = Et — Et = C8 DS C++ AC — D! = AC* + Et — D'
qui fournissent, par (30), les accroissements
| 2]
Api AELN, oA (AE 1
> : ; 7 1
F AES)
mais, à ce moment, il faut repasser aux valeurs elles-mémes par
Bic BY AB}, .., EB =711 + Ar}
Au stade n, on calcule D"' comme au n° 25, puis
AD"? — pr pn?
dont les éléments de seconde espèce fournissent AMY, AMY, ANTI, AN ces valeurs,
introduites dans le système linéaire (27), donnent AC”—* ; l’accroissement correspondant de
E est
AE" = E" — EM = Ch — D
Or, le stade précédent a fourni
AE" x CM Om
On a donc, en soustrayant ces deux égalités membre à membre,
AE" — AE"? — AC — AD"?
ou encore, au stade (n + 1)
AE" = AE" + AC" — AD" (38)
dont les éléments de première espèce fournissent AB” , …, Ar”, et enlin
Bit BE HART D ax, TETE + Af] (39)
/ >
des
la 1
(3)
d'al
mu
lem
gen
à |
ray
soi
por
tên
axe
où