A theoretical analysis of the L-2 least-squares finite element method (LSFEM) for solving the Stokes equations in the velocity-vorticity pressure (VVP) first-order system formulation with the Dirichlet velocity boundary condition is given. The least-squar