该多元非线性方程组不适宜用solve()函数命令来求带扒茄，可以用fsolve()函数命令来求其在某区间内蠢察的此棚数值解。求解过程如下：
func=@(x)[exp((1193*x(3))/20-(543*x(3)^2)/5-(3139*x(3)^(1/2))/100+2867/1000)-x(1),exp((1193*x(4))/20-(543*x(4)^2)/5-(3139*x(4)^(1/2))/100+2867/1000)-x(2),-(x(1)-751/5000)/(x(2)-751/5000)-8994/625,-(x(3)-57/625)/(x(4)-57/625)-8994/625]; %自定义函数
x0=[0.40181 0.075967 0.23992 0.12332]; %初值
x= fsolve(@(x)func(x),x0); %方程组的解
w=x(1),y=x(3),z=x(4),x=x(2)
运行结果
