In modeling various real processes, an important role is played by methods of solution source identification problem for partial differential equation. The current paper is devoted to approximate of elliptic over determined problem with integral condition for derivatives. In the beginning, inverse problem is reduced to some auxiliary nonlocal boundary value problem with integral boundary condition for derivatives. The parameter of equation is defined after solving that auxiliary nonlocal problem. The second order of accuracy difference scheme for approximately solving abstract elliptic overdetermined problem is proposed. By using operator approach existence of solution difference problem is proved. For solution of constructed difference scheme stability and coercive stability estimates are established. Later, obtained abstract results are applied to get stability estimates for solution Neumann-type overdetermined elliptic multidimensional difference problems with integral conditions. Finally, by using MATLAB program, we present numerical results for two dimensional and three dimensional test examples with short explanation on realization on computer.