The method of QCD analysis of singlet and nonsinglet structure functions, based of their Jacobi polynomial reconstruction from perturbative QCD predictions for the Mellin moments, is described. The accuracy and stability of the method are demonstrated with the help of BCDMS hydrogen data. 38 refs.; 5 figs.; 2 tabs.