Salve a tutti,
è il mio post e spero di non violare nessuna regola o di aver sbaglato posto.
Come da titolo la mia domanda riguarda la risoluzione di sistemi lineari con Matlab. Ho creato di mio una funzione che risolve un sistema lineare sfruttando prima la fattorizzazione LU e poi risolvendo con sostituzioni in avanti e indietro. Questo programma che ho scritto dipende, ovviamente, dal condizionamento della matrice in ingresso. Facendo delle prove con matrici hilbertiane e non l'errore sulla soluzione cambia.
Ho provato a fare valutazioni sul peso del condizionamento sfruttando, questa volta, la risoluzione di sistemi già implementata in Matlab (A/b = x). Questo mostra che dopo un certo n (molto piccolo, circa 10) gli errori forniti sono praticamente uguali, tanto da far sì che il grafico si sovrapponga.
E' possibile che sia questo il grafico o c'è un errore che non riesco a trovare nel codice?