SUMMARYThe axial singularity inverse method for designing bodies of revolution has been improved by using higher-order doublet elements. The performance of the method for various element orders and other solution parameters is presented in detail. The results indicate that the method is generally more robust, less sensitive to insets and has a better-conditioned coefficient matrix compared with the source method of the same order. The condition number of the matrix is shown to increase with the thickness of the body, the order of the method, the number of elements and the degree of stretching of the node distribution. In general, good performance is attained for most bodies even with f, as low as 2 by using 10-12 second-order doublet elements with insets greater than 0-02L from rounded ends. Increasing the insets to 0.06L appears to improve the accuracy of the method for most bodies but slows its convergence.