A numerical method is proposed for solving multi-dimensional hyperbolic-parabolic differential equations with the nonlocal boundary condition in t and Dirichlet and Neumann conditions in space variables. The first and second order of accuracy difference schemes are presented. The stability estimates for the solution and its first and second orders difference derivatives are established. A procedure of modified Gauss elimination method is used for solving these difference schemes in the case of a one-dimensional hyperbolic-parabolic differential equations with variable coefficients in x and two-dimensional hyperbolic-parabolic equation.