Shadow cast by an thin emission disk around a Schwarzschild Black Hole

In [1]:
import astropy.units as u

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

Defining the conditions

In [2]:
mass = 1 * u.kg
fov = 30 * u.km 
# What field of view is the user expecting
In [3]:
shadow = Shadow(mass=mass, fov=fov, n_rays=1000)
In [4]:
obj = ShadowPlotter(shadow=shadow, is_line_plot=True)
obj.plot()
obj.show()
In [5]:
obj = ShadowPlotter(shadow=shadow, is_line_plot=False)
obj.plot()
obj.show()