Numerical solutions for unsteady flow in unconfined aquifers