A Gauss-Galerkin finite element method for a class of singular diffusion equations in two space variables