各位帮帮忙，用matlab语言进行如下计算

a = rand(6); % 用你的矩阵替换，应该是m*m的实对称矩阵
a1 = mean(a,1); % 各列的平均值，也就是(1/m)∑aik
a2 = mean(a,2); % 各行的平均值，也就是(1/m)∑akj

[m,n] = size(a); % 行数列数
c1 = repmat(a1,m,1); % 行向量a1重复m行
c2 = repmat(a2,1,n); % 列向量a2重复n列

c = exp(c1 + c2); % 求出c

[v,d] = eig(c) % 特征值和特征向量

[d_max,ind] = max(abs(diag(d))) % 最大特征值和最大特征值所在的列数
v_max = v(:,ind) % 最大特征值对应的列