  In this work we simulate null geodesics for the Bonnor massive dipole metric
by implementing a symbolic-numerical algorithm in Sage and Python. This program
is also capable of visualizing in 3D, in principle, the geodesics for any given
metric. Geodesics are launched from a common point, collectively forming a cone
of light beams, simulating a solid-angle section of a point source in front of
a massive object with a magnetic field. Parallel light beams also were
considered, and their bending due to the curvature of the space-time was
simulated.
