matlab求解方程组数值解

求求大神帮帮忙，方程式为
f=@(x)([x(1)-(4*g*u*x(2)/(3*p*x(3)))^(1/3)
x(2)-p*d*x(1)/u
x(3)-16*(1+0.173*x(2)^0.657)/x(2)+0.413/(1+16300*x(2)^(-0.109))])
怎么调初值也得不出结果，而且我没分了，谁能帮帮我

这个方程式没看懂啊，
能写清楚点

functionf=mymagic(x,b,c,e,f,h,m,n) f=[(m/2-n*sin(c/2)+e*cos(f))^2+(h+n*cos(c/2)-e*sin(f))^2-(m/2-n*sin(c/2+x(2))+e*cos(f-b))^2-(h+n*cos(c/2+x(2))-e*sin(f-b))^2 (m/2-n*sin(c/2)+e*cos(f))^2+(h+n*cos(c/2)-e*sin(f))^2-(m/2+e*cos(x(1)+f)-n*sin(c/2-x(2)))^2-(h+n*cos(c/2-x(2))-e*sin(x(1)+f))^2]; fsolve(@(x)mymagic(x,1,2,3,4,5,6,7),[0;0]) 上面是函数，下面是调用的语句，其中最后的[0;0]是迭代的初值，这里使用了fsolve进行数值求解，求解的方法就是牛顿迭代法！ 祝你学习愉快！