An approach to model the Martian magnetization can be done using classical source models. Classical models, such as uniformly magnetized spheres and cylinders, allow for the introduction of additional constraints related to the available information of the magnetic field and its sources. The use of a suitable conservative numerical scheme in Cartesian coordinates was carried out for numerical studies. In this work the motion of different charged particles under the influence of different magnetized sources have been analyzed by using the proposed numerical scheme. For that purpose, electron, proton and alpha particles were used. In addition, the impact of the gravitational effect on the particles motion was also studied.