Wegen der Symmetrie und positiven Definitheit von
im Falle des
Vollranges von
bietet sich zur Lösung des Normalgleichungssystems das
CHOLESKY-Verfahren an.
Leider handelt es sich dabei um einen numerisch instabilen Algorithmus, der sich jedoch
bei Problemen mit ,,großem`` Residuum
und ,,kleiner`` Lösung
numerisch gutartig verhält.