A B-spline based high order panel method was developed for the motion of bodies in an ideal fluid, either of infinite extent or with a free boundarysurface. In this method, both the geometry and the potential are represented by the B-spline, which guarantees more accurate results than most potential based low order methods. In the present work, we applied this B-spline based high order method to the radiation problem of floating bodies. The boundary condition on the free surface was satisfied by adopting a Kelvin-type Green function and irregular frequencies were removed by placing additional control points on the free surface surrounding the body. The numerical results were validated by comparison with existing numerical and experimental results.