Otrzymany układ można zapisać równoważnie w postaci (112)
12*2 + al2X2 + ' |
*’ nxn |
:22*2 d" a23X3 |
*• + *2«*n |
n2x2 + «n3*3 + |
’“ + annxn |
*1 = Pi + a — @2 + a
gdzie:
h. ail' |
au = |
ait | |
oraz | |||
«t. = |
0 |
dla i = |
■A |
Przyjmując | |||
«11 |
a\2 |
au. | |
a — |
«21 |
aZ2 - |
a2n |
®nl |
an2 - |
ann |
dla t±j t.j = 1,2,—.n.
oraz
układ 1.1.2 można zapisać w postaci macierzowej jako x = /? 4- a • x
Otrzymaną równość wykorzystuje się do uzyskania kolejnych przybliżeń rozwiązania układu. W tym celu, jako rozwiązanie początkowe, obiera się dowolny wektor (np. wektor zerowy) i oblicza się kolejne iteracje xCO = £+«.*«»
(przybliżenie pierwsze),
= /? + a •
(przybliżenie drugie) itd. Stąd ogólny wzór iteracyjny (1.1.3)
XC*+D = /? + *.*<*>
Kolejne przybliżenia x x , ... , x , ... utworzą ciąg wektorów. Jeżeli istnieje granica tego ciągu lim,,.,,*, x^k) = x% wtedy x jest rozwiązaniem układu równań liniowych. Aby powyższa granica istniała, to ciąg wektorów x(0>, x(1), ... musi być ciągiem zbieżnym. Zbieżność tego ciągu jest równoważna zbieżności metody iteracyjnej.
2