The paper presents the solution to the geodetic boundary value problem by the finite element method in area of Slovak Republic. Generally, we have made two numerical experiments. In the first one, Neumann BC in the form of gravity disturbances generated from EGM-96 is used and the solution is verified by the quasigeoidal heights generated directly from EGM-96. In the second one, Neumann BC is computed from gravity measurements and the solution is compared to the quasigeoidal heights obtained by GPS/leveling method.