Verificare cu calculatorul:
sage: A = matrix( QQ, 3,3, [ 1,-1,0, 1,-2,-3, 0,-1,1 ] )
sage: A
[ 1 -1 0]
[ 1 -2 -3]
[ 0 -1 1]
sage: A.det()
-4
sage: A.inverse()
[ 5/4 -1/4 -3/4]
[ 1/4 -1/4 -3/4]
[ 1/4 -1/4 1/4]
sage:
Explicatie:
Peste tot avem blocuri de matrici 3x3 incadrate intre linii,
A, B, C, D, ... si E, F, G, H, ...
care sunt asezate dupa cum urmeaza mai sus:
A E
B F
C G
D H
: :
: :
E X
unde A este matricea cu care am plecat si E este matricea identitate.
Nota: La fiecare pas avem doua matrici in aceeasi linie de matrici 3x3 , sa zicem
U V
cu proprietatea ca U . inversa(V) = A , unde A este matricea cu care plecam.
Dec explicit:
inversa(E) . A = A
inversa(F) . B = A
inversa(G) . C = A
inversa(H) . D = A
: : :
: : :
inversa(X) . E = A .
Cateodata am copiat pur si simplu si am marcat un pivot.
Aplicam algoritmul lui Gauss de eliminare totala de la fiecare pas la urmatorul.
Eliminarea corespunde unei inmultiri din stanga cu o matrice elementara,
acest si numai acest lucru trebuie inteles.
Deci daca plecam cu sau suntem undeva pe drum cu doua matrici
U V
care satisfac inversa(V) . U = A .
atunci (pentru o matrice (inversabila) elementara) si matricile
SU SV
satisfac ce trebuie, anume
inversa( SV ) . SU =
inversa( V ) . inversa( S ) . SU =
inversa( V ) . E . U =
inversa( V ) . U =
A .