<![CDATA[
int main() {
array double a[3][3] = {0.8, 0.2, 0.1,
0.1, 0.7, 0.3,
0.1, 0.1, 0.6};
array double evalues[3],evectors[3][3];
array double p[4];
eigensystem(a, evalues, NULL);
polycoef(p, evalues);
printf("evalues = %f\n", evalues);
printf("characteristic polynomial of matrix a = %f\n", p);
}
]]>
<![CDATA[
evalues = 1.000000 0.600000 0.500000
characteristic polynomial of matrix a = 1.000000 -2.100000 1.400000 -0.300000
]]>