A theoretical analysis of the combined finite volume/least squares approximations to boundary value problems for second-order elliptic equations in mixed first-order system formulation with variable coefficients in two- and three-dimensional bounded domai