Sperimentazioni dell'algoritmo

Proviamo l'algoritmo di risoluzione dei sistemi lineari tramite fattorizzazione $QR$ con la matrice dell'esempio precedente:

» a

a =

     0     7
     5     6

» x=solveQR(a,[1;1])

x =

   0.02857142857143
   0.14285714285714

» inv(a)*[1;1]

ans =

   0.02857142857143
   0.14285714285714

»

ed anche in questo caso abbiamo verificato il risultato ottenuto co nostro algoritmo con il calcolo esplicito del vettore soluzione $x$.