Geodesics on a Gaussian bump, deflected by curvature in analogy to gravitational lensing
This program numerically integrates the geodesic equation on a surface defined by a Gaussian function z = a*exp(-(x^2+y^2)/s^2). The bump creates a region of positive curvature at the peak and negative curvature on the flanks, so geodesics that pass over the bump are deflected in a manner analogous to gravitational lensing. The trajectories are computed via the Christoffel symbols of the induced metric.
Controls: Click and drag to orbit, scroll to zoom.
