ToolPythonBeginner Eigenvalues with numpy.linalg: normal modes of coupled oscillators return np.diag(s[:-1] + s[1:]) - np.diag(s[1:-1], 1) - np.diag(s[1:-1], -1)residual = np.abs(K / m @ V - V @ np.diag(w2)).max()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)))
RecipePythonBeginner Fit a curve to data with error bars and draw a confidence band perr = np.sqrt(np.diag(pcov))
ToolPythonBeginner Numerical integration with scipy.integrate: the area under a measured peak d_slope, d_intercept = np.sqrt(np.diag(C))