MATLAB fsolve：串联的矩阵的维度不一致，如何解决？

Question

a=9;
d=0.6;
xr1=1;
xg1=1;
xb1=1.2;

x0=[2,2,2,2,2,2,0.6];      % 初始参考点

fun=@(x)[

(x(1)+x(2)+x(1)+2*x(7))-a ,...

(x(2)*3^0.5+d+(x(2)*3^0.5+d)/2)-a   ,...
   
(3^0.5*(x(1)*x(2)+x(2)*x(2)/2)/a/a)/(3^0.5*(x(3)*x(4)+x(4)*x(4)/2)/a/a)-xr1/xg1 ,...
(3^0.5*(x(3)*x(4)+x(4)*x(4)/2)/a/a)/(3^0.5*(x(5)*x(6)+x(6)*x(6)/2)/a/a)-xg1/xb1 ,...
(3^0.5*(x(1)*x(2)+x(2)*x(2)/2)/a/a)/(3^0.5*(x(5)*x(6)+x(6)*x(6)/2)/a/a)-xr1/xb1 ,...
    ];
x=fsolve(fun,x0);   %调用函数

Lr1=x(1)
Lr2=x(2)
Lg1=x(3)
Lg2=x(4)
Lb1=x(5)
Lb2=x(6)
Ld =x(7)

错误使用 vertcat
串联的矩阵的维度不一致。

出错 fsolve (line 230)
fuser = feval(funfcn{3},x,varargin{:});

出错 SIX (line 41)
x=fsolve(fun,x0); %调用函数

原因:
Failure in initial objective function evaluation.
FSOLVE cannot continue.