Two meshfree methods are proposed for the numerical solution of boundary value problems (BVPs) for the Laplace equation, coupled with boundary conditions with jump discontinuities. In the first case, the BVP is solved in two steps, using a subtraction of singularity approach. Here, the singular subproblem is solved analytically while the classical method of fundamental solutions (MFS) is applied for the solution of the regular subproblem. In the second case, the total BVP is solved using a variant of the MFS where its approximation basis is enriched with a set of harmonic functions with singular traces on the boundary of the domain. The same singularity-capturing functions, motivated by the boundary element method (BEM), are used for the singular part of the solution in the first method and for augmenting the MFS basis in the second method. Comparative numerical results are presented for 2D problems with discontinuous Dirichlet boundary conditions. In particular, the inappropriate oscillatory behavior of the classical MFS solution, due to the Gibbs phenomenon, is shown to vanish.