This paper approximates a two- or three-dimensional elliptic problem by using a mixed method that approximates the function and its derivative separately, and uses the least squares method to find the solution. Superconvergence results are proved for finite element approximations, and it is shown what regularity conditions and spaces lead to optimal results.