#include <numeric.h> 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)