Numerical Methods for Option Pricing under the Two-Factor Models