A new numerical scheme for structures of rotating magnetic stars

We have developed a new numerical scheme for obtaining structures of rapidly rotating stars with strong magnetic fields. In our scheme, both poloidal and toroidal magnetic fields can be treated for stars with compressibility and infinite conductivity. By introducing the vector potential and its integral representation, we can treat the boundary condition for the magnetic fields across the surface properly. We show structures and distributions of magnetic fields as well as the distributions of the currents of rotating magnetic polytropic stars with polytropic index N= 1.5. The shapes of magnetic stars are oblate as long as the magnetic vector potential decreases as 1/r when r ->infinity. For extremely strong magnetic fields, equilibrium configurations can be of toroidal shapes.