In this paper, wave resistance for a submerged body is calculated by a higher order panel method. The Neumann-Kelvin problem is solved by the source or normal dipole distribution method. The body surface is represented by a bicubic B-spline and the singularity strengths are approximated by a bilinear form. The results calculated by the higher order panel method are compared with those by the lowest order panel method developed by Hess & Smith. The convergence rate of the higher order panel method is much better than the lowest order panel method. But the wave resistance calculated by the higher order panel method still shows discrepancy with an analytic solution at low Froude number like that by the lowest order panel method.