table类型操作数不支持运算符*

x=x1
y=y1
z=z1

eq1=(18*(2x - 30))/(5((y - z)^2 + (x - 15)^2)^(1/2)) + (36x)/(5((y - 5)^2 + x^2)^(1/2));

eq2=(18*(2y - 10))/(5((y - 5)^2 + x^2)^(1/2)) + (18*(2y - 2z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;

eq3=(141*(2z - 16))/(10((z - 8)^2 + 25)^(1/2)) - (18*(2y - 2z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));

[x,y,z]=solve('eq1','eq2','eq3','x','y','z')
x=subs(x),y=subs(y),z=subs(z)

你好，帮你更正，一般符号运算先定义比如syms x y z：

syms x y z
eq1=(18*(2*x - 30))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + (36*x)/(5*((y - 5)^2 + x^2)^(1/2));
eq2=(18*(2*y - 10))/(5*((y - 5)^2 + x^2)^(1/2)) + (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;
eq3=(141*(2*z - 16))/(10*((z - 8)^2 + 25)^(1/2)) - (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));
[x,y,z]=solve([eq1,eq2,eq3],[x,y,z])

符号计算比较复杂，建议使用数值计算：

function main
x0 = rand(3,1); %初值
x = fsolve(@(x)fun(x), x0)
end

function eq = fun(a)
x = a(1);
y = a(2);
z = a(3);
eq = zeros(3,1);
eq(1)=(18*(2*x - 30))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + (36*x)/(5*((y - 5)^2 + x^2)^(1/2));
eq(2)=(18*(2*y - 10))/(5*((y - 5)^2 + x^2)^(1/2)) + (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2)) + 36/5;
eq(3)=(141*(2*z - 16))/(10*((z - 8)^2 + 25)^(1/2)) - (18*(2*y - 2*z))/(5*((y - z)^2 + (x - 15)^2)^(1/2));
end

最后得到结果：

x =

   5.459266184633963
   1.848091223289034
   7.356436516679399

有帮助望采纳哟