#include
int main() {
array double complex a[3][3] = {0.8, 0.2, 0.1,
0.1, 0.7, 0.3,
0.1, 0.1, 0.6};
array double complex a2[3][3] = {complex(0.8,-1), 0.2, 0.1,
0.1, 0.7, 0.3,
0.1, 0.1, 0.6};
array double complex evalues[3];
array double complex evectors[3][3];
eigen(a, evalues, evectors);
printf("eigenvalues of a =\n%f\n", evalues);
printf("eigenvectors of a =\n%f\n", evectors);
eigen(a2, evalues, evectors);
printf("eigenvalues of a =\n%f\n", evalues);
printf("eigenvectors of a =\n%f\n", evectors);
}

eigenvalues of a =
complex(1.000000,0.000000) complex(0.600000,0.000000) complex(0.500000,0.000000)
eigenvectors of a =
complex(0.744845,0.000000) complex(-0.707107,0.000000) complex(-0.408248,0.000000)
complex(0.579324,0.000000) complex(0.707107,0.000000) complex(0.816497,0.000000)
complex(0.331042,0.000000) complex(0.000000,0.000000) complex(-0.408248,0.000000)
eigenvalues of a =
complex(0.796911,-0.968453) complex(0.832118,-0.037024) complex(0.470970,0.005477)
eigenvectors of a =
complex(0.989164,0.000000) complex(0.003912,-0.225655) complex(-0.025122,-0.090331)
complex(-0.019749,0.106886) complex(0.891481,0.000000) complex(0.794571,0.000000)
complex(0.008946,0.098280) complex(0.391299,-0.034802) complex(-0.598227,0.044616)