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.