  We analyze a divergence based first order system least squares method applied
to a second order elliptic model problem with homogeneous boundary conditions.
We prove optimal convergence in the $L^2(\Omega)$ norm for the scalar variable.
Numerical results confirm our findings.
