The paper proposes the numerical method of solution the problems of calculation the stress state in thick-walled cylinders and spheres from physically nonlinear inhomogeneous material. The urgency of solved problem due to the change of mechanical properties of materials under the influence of different physical fields (temperature, humidity, radiation, etc.). The deformation diagram describes the three-parameter formula. The numerical method used the method of successive approximations. The results of numerical calculation are compared with the test analytical solutions obtaining the authors with some restrictions on diagram parameters. The obtained results can be considered quite satisfactory.