In this thesis we have sought to solve Richards’ equation. In time we applied a backward Euler discretization and in space we used a continuous Galerkin finite element discretization. The focal point was how to efficiently ...