The simple low-dimensional model proposed by Guzdar and co-workers has been extended and applied to study the nonlinear diamagnetic effect on zonal flow in interchange turbulence. Based on the numerical solution to the derived nonlinear equations, we can draw the conclusion: diamagnetic effect suppresses the enhancement of the normal polarized drift nonlinearity, which could lead to the damping of zonal flow.