[口袋里的long]

最近在求解一个非线性方程组，用到fsolve函数，但是不会用啊，百度下根据人家语句直接把数换了，但是解不出结果来特来。
function y = myfun( x )
a=x(1);
b=x(2);
y(1)=exp(-(5775500/a)^b)-163/10340;
y(2)=a*gamma(1+1/b)-5745500/163;

h=fsolve('myfun',[88821,2],optimset('display','off'))
求解结果为
h =
1.0e+04 *
8.8821 0.0002
这完全就是我给的初值嘛，不知道为什么会出现这个结果，烦请大神们指教下如何求解这个非线性方程组啊？？？？