CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutSign UpSign In

CoCalc provides the best real-time collaborative environment for Jupyter Notebooks, LaTeX documents, and SageMath, scalable from individual users to large groups and classes!

| Download

Reference manual of the SageMath package kerrgeodesic_gw

Project: Sagittarius
Views: 7899
License: GPL3
Image: ubuntu2004
Tweet
Share
Share
1
from kerrgeodesic_gw import KerrBH

2
a = var('a')

3
M = KerrBH(a)

4
M.boyer_lindquist_coordinates()

5
init_point = M((0, 6, pi/2, 0), name='p0')

6
geod = M.geodesic([0, 300], init_point, mu=1, E=0.883, \

7
                  L=1.982, Q=0.467, a_num=0.998)

8
geod.integrate(step=0.005)

9
graph = geod.plot(coordinates='txy', label_axes=True)

10
sphinx_plot(graph)

11