Shadow cast by a thin emission disk around a Schwarzschild black hole¶

[1]:

import astropy.units as u

from einsteinpy.rays import Shadow
from einsteinpy.plotting import ShadowPlotter

Defining the conditions

[2]:

mass = 1 * u.kg
fov = 30 * u.km
# What field of view is the user expecting

[3]:

shadow = Shadow(mass=mass, fov=fov, n_rays=1000)

[4]:

obj = ShadowPlotter(shadow=shadow, is_line_plot=True)
obj.plot()
obj.show()

[5]:

obj = ShadowPlotter(shadow=shadow, is_line_plot=False)
obj.plot()
obj.show()

EinsteinPy