Parabolic orbits are not working

[nielsw2]

It looks like parabolic orbits (e=1) are not working:

orbit = KerrGeoOrbit[0, 11, 1, 1]

returns a bunch of 1/0 errors.

[MvdMeent]

Is it just Schwarzschild?

[MvdMeent]

Actually, I don't think we ever implemented parabolic orbits, did we?

[MvdMeent]

Apparently, I started doing this at some point. I'll investigate.

[MvdMeent]

After looking at it, I remember why I previously abandoned this. Taking the e->1 limit can be done relatively easily for many parts of the solution. For example, the polar solution reverts to the trivial Schwarzschild solution. The however, the limit of the t solution appears non-trivial. My first attempt yielded a divergent answer. Probably the best route is to take the limit of the EOMs and then solve the resulting elliptic integrals.