你的精确解带入等式后并不成立。matlab的键或fsolve,lsqnonlin都可以求解非线性方程稿源伍组:
f = @(x, y, z)[x^y+y^x-5*x*y*z-85; x^3*y^z*z^y-60; x^z+z^x-y-1];
% 初始值设为裂旁区间中点
fsolve(@(x)f(x(1), x(2), x(3)), [4 3 1.25]) % fsolve求解
lsqnonlin(@(x)f(x(1), x(2), x(3)), [4 3 1.25]) % lsqnonlin求解