We consider the Cauchy problem for Laplacian. Using the single layer representation, we obtain an equivalent system of boundary integral equations. We show the singular values of the ill-posed Cauchy operator decay exponentially, which means that a small error is exponentially amplified in the solution of the Cauchy problem. We show the decaying rate is dependent on the geometry of he domain, which provides the information on the choice of numerically meaningful modes. We suggest a pseudo-inverse regularization method based on singular value decomposition and present various numerical simulations.