INPUT:
Either:
- R,vR,vT,z,vz[,phi]:
- floats: phase-space value for single object (phi is optional)
- numpy.ndarray: [N] phase-space values for N objects
- Orbit instance: initial condition used if that’s it, orbit(t) if there is a time given as well as the second argument