The formulation and development of a numerical algorithm for solving shell theory problems concerning cylindrical shells with variable-stiffness reinforcements are considered. The Hamilton-Ostrogradskii variational principle is used to obtain the equation of motion. The algorithm is based on finite-difference approximation of the initial variational functional. A numerical example is presented.