This article is devoted to the algorithm of the forecast of stress and strain state nonho-mogeneous layered physical nonlinear rock massif. The modeling based on applying the whole package of computational methods: variational method, discrete continuation on numeric parameter method, quasi-linearization of nonlinear value boundary problem, finite difference method, Thomas algorithm and iterative process.

Development of computer engineering in mathematical modeling gives a possibility to use the numerical methods parallel with the analytical methods. Among them finite difference method, finite elements method, boundary element method and others. FDM is the effective one, which can be used in geomechanical problems extensively.