def METHOD_NAME(self):
    iwork = self._lsoda_solver._integrator.iwork
    rwork = self._lsoda_solver._integrator.rwork
    order = iwork[14]
    h = rwork[11]
    yh = np.reshape(rwork[20:20 + (order + 1) * self.n],
                    (self.n, order + 1), order='F').copy()
    return LsodaDenseOutput(self.t_old, self.t, h, order, yh)