def max_abs_eigs(q, d): """ Returns the largest of the absolute values of the eigenvalues of system matrix at q in direction d """ J = jacobian(q, d) return max(abs(eigvals(J)))