The deconvolution problem of the Constituent Quark Model (CQM) in finding the partondistribution function in valon at Q₀² is reduced to an optimum problem and solved successfully with Newton-Raphson method. Moreover, by discretizing the QCD evolution eqs into ordinary simultaneous differential eqs, we obtain the parton distribution functions in valon at arbitary Q² > Q₀², and thus, through convolution, we finally find the Q²-dependent nucleon parton distribution functions in CQM picture.

To avoid the dependence of hadron structure functions upon supposed function forms and to improve the accuracy, a new method to find the hadron structure fundtion at high Q² is discussed in this paper to replace the commonly used moment method. The sturcture functions for gluon and valence quark and sea at high Q² are obtained by discretizing the A-P evolution equations into ordinary simultaneous differential eqs. and then solving them with Runge-Kutta method.