def METHOD_NAME(self, w0, t, w):
    """ """
    if w.shape[-1] == 1:
        w = w[..., 0]
    pos_unit = self._func_units["length"]
    t_unit = self._func_units["time"]
    vel_unit = pos_unit / t_unit
    from ..dynamics import Orbit
    orbit = Orbit(
        pos=w[:self.ndim] * pos_unit,
        vel=w[self.ndim:] * vel_unit,
        t=t * t_unit,
    )
    return orbit