RecipePythonBeginner Draw a phase portrait of a two-variable ODE system if inside and not any(np.allclose(p, q) for q in fixed): # many guesses find the same point
ToolPythonBeginner Eigenvalues with numpy.linalg: normal modes of coupled oscillators print("K @ x0 equals F:", np.allclose(K @ x0, F))print("V[:, 0] equals the textbook vector:", np.allclose(V[:, 0], textbook))print("A X equals X diag(w2): ", np.allclose(A @ X, X @ np.diag(w2_co2)))print("X^T M X equals 1: ", np.allclose(X.T @ np.diag(m_co2) @ X, np.eye(3)))