An accurate method for potential flow about a body of revolution has been developed. Solutions for both axial and transverse motions are obtained by first solving integral equations of the second kind for the surface source distributions. These source distributions are then integrated over the surface to get velocity potentials. A simple differentiation of this potential gives the velocity.
Simpson's rule is adopted to discretize the integral equation. Each interval is subdivided to treat the rapidly varing kernel more accurately. Source values at corresponding points are represented by the values at the primary node points through interpolation.
The results for spheroids are compared with exact solutions and show a marked improvement over the Gaussian quadrature formula.