You likely don’t want to be using the normal equations. They lose a lot of numeric stability compared to just solving directly (which will likely use a QR factorization).
Specifically, rather than using Sk'*(Rk\Sk) you can use X0\Yk
4 Likes