Magnetized Stars in General Relativity*

While magnetic fields are ubiquitous features of stellar systems, it is only relatively recently that the evolution of compact objects such as neutron stars have been simulated using both general relativity and ideal magnetohydrodynamics. An important aspect of these simulations is beginning with physically reasonable initial data. To this end, we consider the problem of constructing equilibrium configurations of axisymmetric, polytropic stars with strong magnetic fields in general relativity.