Application of the least square method to the solving linear differential-algebraic equations

We consider application of the least square method to the numerical solution of a linear system of ordinary differential equations (ODEs) with an identically singular matrix multiplied a higher derivative by the desired vector-function. We discuss the behavior of the gradient method for minimizing the functional of the residual square in the Sobolev space and some other issues. The results of the numerical experiments are given.