En este trabajo se simulan geodésicas nulas para la métrica de dipolo masivo de Bonnor, implementando un algoritmo simbólico-numérico en Sage y Python. Este programa es capaz, en principio, de visualizar en 3D las geodésicas de cualquier métrica. Estas geodésicas inician en un punto común, formando colectivamente un cono de rayos de luz, simulando una sección de ángulo sólido de una fuente puntual frente a un objeto masivo con campo magnético. Se consideran también rayos de luz paralelos, y se simuló su cambio de trayectoria debida a la curvatura del espacio-tiempo.