This paper considers the numerical computation of added resistance on ships in the presence of incident waves. As a method of solution, a higher-order Rankine panel method is applied in time domain. The added resistance is evaluated by integrating the second-order pressure on the body surface. Computational results are validated by comparing with experimental data and other computational results on a hemi-sphere, a barge, Wigley hull models, and Series 60 hull, showing very fair agreements. The study is extended to the comparison between Neumann-Kelvin and double-body linearization approaches, and their differences are discussed.