Le cercle noir de rayon \(R_s\) représente l’horizon des évènements du trou noir.
Les photons arrivent de la droite (\(x=+\infty\) et \(y=b\)) avec une valeur du paramètre d’impact \(b\) qui donne une déviation de \(4\pi\) (2 tours complets).
Après une intégration numérique suivant \(\varphi\), la figure est tracée avec les coordonnées cartésiennes \(x=r\cos\varphi\) et \(y=r\sin\varphi\), en prenant un hypothétique trou noir de la masse du soleil \(M\odot\)
soit \(R_s=2 \ 953\ m\), \(b_{crit}=7\ 672,73\ m\) et avec \(b_{4\pi} = 7\ 672,75\ m\).
La valeur de \(b\) peut être obtenue en choisissant la valeur de la variation finale \(\Delta\varphi\) = déviation + \(\pi\) (valeur cible dans un tableur Excel ou donnée d’entrée dans un script Python).
Ici, elle vaut \(b_{crit} + 2\ cm\).
Une déviation de \(8\pi\) (4 tours complets) est obtenue pour \(b_{8\pi} = b_{crit} + 7\ 10^{-8}\ m\).
Pour un trou noir d’une autre masse \(M\), les tracés sont conservés en appliquant le facteur d’échelle \(\frac{M}{M\odot}\).