An improved numerical scheme for the approximate solution of the Parabolic Wave model