In this paper, we consider the pricing problem for the extremum options by constructing a double nonaffine stochastic volatility model. The joint characteristic function of the logarithm of two asset prices is derived by using the Feynman–Kac theorem and one-order Taylor approximation expansion. The semiclosed analytical pricing formulas of the European extremum options including option on maximum and option on minimum of two underlying assets are derived by using measure change technique and Fourier transform approach. Some numerical examples are provided to analyze the pricing results of extremum options under affine model, nonaffine model, Black–Scholes model, and the influences of some model parameters on the option. Numerical results show that the analytical pricing formulas have higher computational efficiency and accuracy than those of Monte Carlo simulation method. Also, results of sensitivity analysis report that the nonaffine models are more effective than other existing models in capturing the effect of volatility on option pricing.