[求助]fsolve可不可以解我这个二元非线性方程组

function f=myfunction2(alpha,zta,y,n,r)
z=-y.^alpha;
z=1-exp(z);
h=r*log(alpha*zta)+(alpha-1).*sum(log(y(1:r)))-sum(-z(1:r))+(zta-1).*sum(log(z(1:r)))+(n-r).*log(1-z(r).^zta);
f(1)=diff(h,alpha);
f(2)=diff(h,zta);



n=30;r=20;
alpha0=0.3;zta0=2;
for j=1:n
y(j)=-log(1-rand(1).^(1/zta0));
y(j)=y(j).^(1/alpha0);
end
y=sort(y);
x0=[0.3,2];
syms alpha zta;
[x,fval]=fsolve(@myfunction2,x0,alpha,zta,y,n,r)

是一个参数的最大似然估计，不知道怎么解出alpha,zta.请高手指点，谢谢！！！