the magnetic field, the centrifugal force, and the radiation drag, and so on as seen in Table \[tab:00\]. To consider the screening effect and the formation of the field-aligned electric field due to inertia, the time step also has to be much smaller than the the period of plasma oscillation. This requirement is also satisfied as seen in Table \[tab:00\]. According to Table 2, the Larmor radii and inertial lengths are enough smaller than the typical scale of the system. The ratio of the plasma frequency $\omega_p$ to the gyro frequency $\omega_g$ changes in such a way that $\omega_p / \omega_g \ll 1$ at the magnetic poles, $\omega_p / \omega_g < 1$ at the light cylinder, and $\omega_p / \omega_g > 1$ at the outer boundary in both simulation and actual pulsars. Thus, the magnetization property of plasma is also simulated well in the simulation and actual pulsars.

Results
=======

  ----------------------------- -----------------------------
   (80mm,80mm)[figure1a.eps]{}   (80mm,80mm)[figure1b.eps]{}
   (80mm,80mm)[figure1c.eps]{}   (80mm,80mm)[figure1d.eps]{}
  ----------------------------- -----------------------------

Magnetospheric plasmas are represented by $\sim 1 \times 10^5$ super-particles. After integration of $\sim 10$ periods of time, the amount of particles produced at the stellar surface and the gaps, and the amount of removed particles are balanced, so that the total 