We developed a new particle simulation code which includes pair creation (magnetic pair creation and photon collision process), propagation of gamma-ray, inertia of particle, interaction of plasma and multi-pole stellar field for steady axisymmetric pulsar magnetosphere. The photon path is solved stochastically by an analytical solution of the mean free path of pair creation processes at the photon position. The superimposed quadrupole magnetic field forms asymmetric electrostatic clouds on tje poloidal plane and the accelerating region is different from the dipole case. Here, we demonstrate some results of a test run for our simulation. We will adopt the code for more complicated cases, such that all above-mentioned effects will be considered together in future work.