ConceptPythonIntermediate Polynomial roots are eigenvalues: why Wilkinson's polynomial loses its roots print(f"eigenvalues: {np.sort(np.linalg.eigvals(C3).real)}")diff = np.abs(np.linalg.eigvals(C) - np.roots(coef)).max()
ConceptPythonIntermediate Stiffness: why an explicit solver crawls on a reaction that has long settled ev = np.sort(np.linalg.eigvals(robertson_jac(t, radau.sol(t))).real)lam_fast = np.abs(np.linalg.eigvals(J)).max(axis=1) # eigvals works on a stack of matrices