  We perform MHD simulations of a thin resistive and viscous accretion disk
around a neutron star with the surface dipolar magnetic field of 10$^8$ Gauss.
The system evolution is followed during the interval of 500 millisecond pulsar
rotations. Matter is accreted through a stable accretion column from the disk
onto the star. We also show propagation of the stellar wind through the corona.
Analysis of the mass accretion flux and torques on the star shows that the disk
reaches the quasi-stationary state.
