QR a comme argument une matrice carrée numérique A
d’ordre n et deux noms de variables : var1 et var2.
QR factorise numériquement cette matrice sous la forme Q*R où
Q est une matrice orthogonale (tQ*Q=I) et R est une matrice triangulaire supérieure.
QR(A) renvoie R, met Q dans var1
et met R dans var2.
On a Q=A*inv(R).
On tape :
On obtient la matrice R :