Abstract. This paper presents the numerical solution for the option on the maximum of two assets using Improving Modified Gauss-Seidel (IMGS) iterative method. Actually, this option can be governed by two-dimensional BlackScholes partial differential equation (PDE). The Crank-Nicolson scheme is applied to discretize the Black-Scholes PDE in order to derive a linear system. Then, the IMGS iterative method is formulated to solve the linear system. Numerical experiments involving Gauss-Seidel (GS) and Modified Gauss-Seidel (MGS) iterative methods are implemented as control methods to test the computational efficiency of the IMGS iterative method.