A mathematical method of dealing with the imposed boundary conditions with multivariate is presented based upon the linear space theory. It can be applied widely to the finite element method and other numerical variational methods. The merit of this approach is that the main rigid matrix remains unchanged in the process of dealing with the constrained conditions. The sufficient condition of the existance and unigueness of solution is that the original main matrix is positive.