The purpose of study is to develop an algorithm for the analysis of thin shells of revolution based on the hybrid formulation of finite element method in two dimensions using a quadrilateral fragment of the middle surface as a discretization element. Nodal axial forces and moments, as well as components of the nodal displacement vector were selected as unknown variables. The number of unknowns in each node of the four-node discretization element reaches nine: six force variables and three kinematic variables. To obtain the flexibility matrix and the nodal forces vector, a modified Reissner functional was used, in which the total specific work of stresses is represented by the specific work of membrane forces and bending moments of the middle surface on its membrane and bending strains, and the specific additional work is determined by the specific work of membrane forces and bending moments of the middle surface. Bilinear shape functions of local coordinates were used as approximating expressions for both force and displacement unknowns. The dimensions of the flexibility matrix of the four-node discretization element were found to be 36×36. The solution of benchmark problem of analyzing a truncated ellipsoid of revolution loaded with internal pressure showed sufficient accuracy in calculating the strength parameters of the studied shell.