Fourth-order difference method for quasilinear Poisson equation in cylindrical symmetry