INPUT:
Either:
- R,vR,vT,z,vz:
- floats: phase-space value for single object
- numpy.ndarray: [N] phase-space values for N objects
- numpy.ndarray: [N,M] phase-space values for N objects at M times
- Orbit instance or list thereof; can be integrated already