Simulations of Axisymmetric Magnetospheres of Neutron Stars