gnuplot_surface, a FreeFem++ code which solves the steady Laplace equation over a disk, and then writes the mesh and solution to a file so that gnuplot() can create a surface plot.
The information on this web page is distributed under the MIT license.