Closed-Null-Geodesics-2D. This code computes elliptic LCSs in 2D, and can be used to compute general null-geodesics with minor modifications.