  Magnetic fields play a crucial role in many astrophysical scenarios and, in
particular, are of paramount importance in the emission mechanism and evolution
of Neutron Stars (NSs). To understand the role of the magnetic field in compact
objects it is important to obtain, as a first step, accurate equilibrium models
for magnetized NSs. Using the conformally flat approximation we solve the
Einstein's equations together with the GRMHD equations in the case of a static
axisymmetryc NS taking into account different types of magnetic configuration.
This allows us to investigate the effect of the magnetic field on global
properties of NSs such as their deformation.

High Energy Astrophysical Phenomena;General Relativity and Quantum Cosmology