用 symbolic math toolbox
其实这个问题手算可得到结果。

你看看这个可否？
syms x
a=(1-x)*eye(4)+x*ones(4,1)*ones(4,1)'

a =

[ 1, x, x, x]
[ x, 1, x, x]
[ x, x, 1, x]
[ x, x, x, 1]
inv(a)

ans =

[ -(2*x+1)/(3*x^2-2*x-1),        x/(3*x^2-2*x-1),        x/(3*x^2-2*x-1),        x/(3*x^2-2*x-1)]
[        x/(3*x^2-2*x-1), -(2*x+1)/(3*x^2-2*x-1),        x/(3*x^2-2*x-1),        x/(3*x^2-2*x-1)]
[        x/(3*x^2-2*x-1),        x/(3*x^2-2*x-1), -(2*x+1)/(3*x^2-2*x-1),        x/(3*x^2-2*x-1)]
[        x/(3*x^2-2*x-1),        x/(3*x^2-2*x-1),        x/(3*x^2-2*x-1), -(2*x+1)/(3*x^2-2*x-1)]
我不知道a表达式有没有写错，你看一下

3# ywh19860616

上路了，但是不是他所要求的东西。代码及结果如下：

#macro2.cn

哦，看错了，不好意思，我以为求逆矩阵，原来求特征矩阵,用eig函数，[v,d]=eig(a),v表示特征值，d为特征向量