如何用matlab求矩阵的特征多项式和特征值啊,例如我生成的是一个4阶魔术矩阵

magic(4)
ans =
16 2 3 13
5 11 10 8
9 7 6 12
4 14 15 1
>> [d,v]=eig(ans)
d =
-0.5000 -0.8236 0.3764 -0.2236
-0.5000 0.4236 0.0236 -0.6708
-0.5000 0.0236 0.4236 0.6708
-0.5000 0.3764 -0.8236 0.2236
v =
34.0000 0 0 0
0 8.9443 0 0
0 0 -8.9443 0
0 0 0 0.0000