The paper considers the application of the method of direct integration to calculations of annular plates and slabs on a continuous variable elastic base. Ring-shaped plates with variable geometric and mechanical parameters are increasingly used. Not only the elastic base, but also the plate thickness and cylindrical stiffness can be variable parameters here. The need for an analytical method for calculating such structures raises no doubts, since it makes it possible to evaluate the accuracy of finite-element analysis. To date, there are no proposals in the literature regarding a general analytical method for the calculation of annular plates on a variable elastic base. A detailed description of the algorithm of the direct integration method is not given in the paper, and all the calculation formulas for the annular plate are taken from the authors' already published article. The results of numerical implementation of this algorithm for specific examples are considered: a concrete plate, which is rigidly pinch on the inner contour, and its outer contour is free, and a steel plate, which is rigidly pinch on the outer contour, and its inner contour is free. To estimate the results of calculation by the author's method, computer modeling of the considered structures in PC LIRA-SAPR and their calculations by the finite-element method have been executed. The foundation reaction is described by Winkler model with a variable bedding factor. In the first case a bed factor is assumed constant, and in the second case it changes under the linear law. Calculations have shown that discrepancy between deflections calculated by the finite-element method and the author's method does not exceed 1 %, and the results of radial and circumferential moments calculation differ more considerably, amounting to 10 %. The authors explain this difference by the inaccuracy of the numerical analysis associated with the semi-automatic method of constructing a finite-element mesh, which should be made finer. The densification of the mesh in the manual mode of its partitioning significantly reduces the discrepancy between the results of calculating the deflections, radial and circumferential bending moments by the finite-element method and the author's method.