Fourth order finite difference methods for the system of 2-d nonlinear elliptic equations with variable coefficients