#define n 3

/* Calculation of eigenvalues (D) and eigenvectors (V) of 3x3 matrix A */
void eigen_decomposition(double A[n][n], double V[n][n], double d[n]);