P-stable boundary value methods for second order IVPs