A numerical method is proposed for solving Bitsadze-Samarskii-type nonlocal boundary value problem for multidimensional elliptic partial differential equations. The second and fourth orders of accuracy stable difference schemes are presented. The stability and almost coercive stability of these difference schemes are established. The method is illustrated by numerical examples.