The reconstruction problem for steady symmetrical two-dimensional magnetic reconnection is addressed in the frame of a two-fluid approximation with neglected ion current. This approach yields Poisson's equation for the magnetic potential of the in-plane magnetic field, where the right-hand side contains the out-of-plane electron current density with the reversed sign. In the simplest case of uniform electron temperature and number density and neglecting the electron inertia, Poisson's equation turns to the Grad-Shafranov one. With boundary conditions fixed at any unclosed curve (the satellite trajectory), both equations result in an ill-posed problem. Since the magnetic configuration in the reconnection region is highly stretched, one can make use of the boundary layer approximation; hence, the problem becomes well-posed. The described approach is generalized for the case of nonuniform electron temperature and number density. The benchmark reconstruction of the PIC simulations data has shown that the main contribution for inaccuracy arises from replacing Poisson's equation by the equation of Grad-Shafranov. Under this substitution, the reachable cross-size of the reconstructed region is shrinking down to fractions of the proton inertial length. Artificial smoothing, demanded by solving the ill-posed problem, boundary layer approximation represent two alternative methods of problem regularization. In terms of the reconstruction error, they perform nearly the same; the second method benefits from the comparative simplicity and less restrictions imposed on the boundary shape.