| Download
SageMath notebooks associated to the Black Hole Lectures (https://luth.obspm.fr/~luthier/gourgoulhon/bh16)
Project: BHLectures
Views: 20094Kernel: SageMath 9.1.beta4
Schwarzschild null geodesics: elliptic integrals
This Jupyter/SageMath notebook is relative to the lectures Geometry and physics of black holes.
To run it, you must start SageMath with sage -n jupyter
.
In [1]:
'SageMath version 9.1.beta4, Release Date: 2020-02-13'
In [2]:
Function
We evaluate via the functions elliptic_f and elliptic_kc of SageMath, taking into account that the second argument of elliptic_f
and the argument of elliptic_kc
is and not .
In [3]:
Case
In [4]:
In [5]:
In [6]:
Case
In [7]:
In [8]:
Asymptotic direction from some emission point
In [9]:
In [10]:
In [11]:
In [12]:
In [13]:
In [14]:
In [15]:
In [16]:
verbose 0 (3797: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 800 points.
verbose 0 (3797: plot.py, generate_plot_points) Last error message: ''
In [17]:
verbose 0 (3797: plot.py, generate_plot_points) WARNING: When plotting, failed to evaluate function at 800 points.
verbose 0 (3797: plot.py, generate_plot_points) Last error message: ''
In [18]:
In [19]: