Ray tracing a black hole with solve_ivp: the shadow of M87*

Afterwards you can trace light rays past a black hole with solve_ivp and events, find the critical impact parameter of its shadow, and render the lensed image.