There are 2 repositories under poincaresurfaceofsection topic.
This mathematica code numerically solves the Hamilton and geodesic equations in Slowly-rotating Kerr