xxx
n=[150000,200000,140000];
syms x y z
eq1='diff(sum(log(x)+(x-1)*log(n(1)-z)-x*log(y)-((n(2)-z)/y)^x)+sum(-((n(3)-z)/y)^x),x)=0';
eq2='diff(sum(log(x)+(x-1)*log(n(1)-z)-x*log(y)-((n(2)-z)/y)^x)+sum(-((n(3)-z)/y)^x),y)=0';
eq3='diff(sum(log(x)+(x-1)*log(n(1)-z)-x*log(y)-((n(2)-z)/y)^x)+sum(-((n(3)-z)/y)^x),z)=0';
solve(eq1,eq2,eq3)
x=ans.x
y=ans.y
z=ans.z
这样就可以了,
你的维数不一致怎么计算呢,兄弟。
n是数组,x是符号变量,当然报错。
|