Numerical methods for ultraparabolic equations