This lecture closes the discussion on numerical accuracy and starts the first section on a numerical method: Direct solvers for linear systems of equations (LSE). After a wrap-up of some mathematical preliminaries we introduce the well-known Gaussian eleminiation procedure to solve LSE.