Numerical analysis: Gaussian elimination