In[1], the exact analytic method for the solution of differential equation with variable coefficients was suggested and an analytic expression of solution was given by initial parameter algorithm. But to some problems such as the bending, free vibration and buckling of nonhomogeneous long cylinders, it is difficult to obtain their solutions by the initial parameter algorithm on computer. In this paper, the substructure computational algorithm for the exact analytic method is presented through the bending of non-homogeneous long cylindrical shell. This substructure algorithm can he applied to solve the problems which can not he calculated by the initial parameter algorithm on computer. Finally, the problems can he reduced to solving a low order system of algehraic equations like the initial parameter algorithm Numerical examples are given and compared with the initial para-algorithm at the end of the paper, which confirms the correctness of the substructure computational algorithm.