Abstract: The issue of data completion is important for the elliptic type partial differential equation. In the inverse Cauchy problem, we need to complete the boundary data by over-specifying Dirichlet and Neumann data on a portion of the boundary. In this paper, we numerically solve the generalized inverse boundary value problems of Laplace equation in a rectangle with one boundary function and two boundary functions missing, which are more difficult than the inverse Cauchy problem. By using the technique of a boundary integral equation method together with a specially designed Trefftz test function, we can complete the boundary data by requiring minimal extra data. Then solving the Laplace equation with the given data and recovered data by the multiple-scale Trefftz method, we can find the numerical solution in the interior nodal points.