求如何用C语言复数矩阵求特征值和特征向量

请教大家，如何用C语言求出一个复数矩阵的特征值和特征向量呢？
比如4X4的复数矩阵为：
5.8751 5.8774 + 1.5423i 2.0818 + 2.8963i 1.8692 + 1.9634i
5.8774 - 1.5423i 6.8258 2.6390 + 2.8255i 2.2955 + 1.9349i
2.0818 - 2.8963i 2.6390 - 2.8255i 3.9534 2.9155 - 0.3039i
1.8692 - 1.9634i 2.2955 - 1.9349i 2.9155 + 0.3039i 2.2816

他的特征向量为：
-0.0506 + 0.2898i 0.1137 + 0.6431i -0.1009 - 0.3856i 0.3665 + 0.4403i
-0.0942 - 0.2603i -0.3552 - 0.5019i -0.3791 - 0.1255i 0.4927 + 0.3782i
-0.5045 - 0.0249i 0.2804 - 0.2056i 0.6644 - 0.0427i 0.4221 - 0.0503i
0.7627 -0.2731 0.4882 0.3245

特征值为：
0.0412 0 0 0
0 0.1996 0 0
0 0 2.0591 0
0 0 0 16.6361

急求！多谢了！