  In this work, we study the effects of different magnetic field configurations
in neutron stars described by a many-body forces formalism (MBF model). The MBF
model is a relativistic mean field formalism that takes into account many-body
forces by means of a meson field dependence of the nuclear interaction coupling
constants. We choose the best parametrization of the model that reproduces
nuclear matter properties at saturation and also describes massive neutron
stars. We assume matter to be in beta-equilibrium, charge neutral and at zero
temperature. Magnetic fields are taken into account both in the equation of
state and in the structure of the stars by the self-consistent solution of the
Einstein-Maxwell equations. We assume a poloidal magnetic field distribution
and calculate its effects on neutron stars, showing its influence on the
gravitational mass and deformation of the stars.
