Ray tracing black holes in numpy This is a super simple 100 lines python script to ray-trace black holes. Ray tracing in Schwarzschild metric is equivalent to propagating in euclidian geometry massive particles subject to the force $$F(r)=-\frac{3}{2}h^2\frac{\hat{r}}{r^5}$$ where $h$ is the angular momentum of the particle. We propagate the particles from the camera as follow : Galery