Based on classical finite volume method for solving two dimensional elliptic equation with discontinuous efficient, a modified finite volume method is presented in this paper. The finite difference method is cooperated with for dealing with some grid point around boundary, and more Taylor expansion items of flux function are intercepted in order to obtain higher numerical accuracy. There is first order accuracy at interface and second order accuracy in the whole domain. Numerical results indicate that the modified finite volume method is effective.