Ray tracing black holes in Schwarzschild metric.
The equations of motion of the light beams are integrated using the Runge Kutta method.
The accretion disk is purely generative, no textures are used.
Rendering may take a while since we need to solve the EOM in Schwarzschild for each pixel of the image.