The numerical implementation and stability of first- and second-order absorbing boundary conditions for simulating free surface gravity waves are considered. The free surface waves are solved with a panel method. The stability of the boundary conditions is proved, assuming certain properties of the integral operator. Arguments are given that support these assumptions. The theoretical results are confirmed in the test cases and the boundary conditions give low reflections for plane waves.