A method of computing Laplacian potentials for rotationally symmetric Dirichlet boundary-value problems is developed. The algorithm is valid forNsurfaces of revolution of arbitrary cross-section, and gives the potentials both inside the surfaces of revolution. The solution is obtained by solving the integral equation formulation using a step approximation to the source and point matching for testing. Sample computations are given, and a general computer program is available.