A non-linear boundary value problem for Laplace's equation on a polygonal domain is solved with a Nystrom method applied to a boundary integral equation reformulation. The method uses meshes which are graded toward the corners. Numerical examples are included.