解答非线性方程组 

--------------------------------------------------------------------------------
ri=1.26 ; re=1.31 ; Dw=2.381 ; Ni=0 ; Ne=3000 ; alf0=20; dm=10.7-2.381;Ndeta=0.9;

syms x y vxq vzq
fai=[0 pi/3 2*pi/3 pi 4*pi/3 5*pi/3 ];
B=(ri+re)/Dw-1;
for i=1:6
cosfaiq=cos(fai(i));
axq=B*Dw*sin(20/180*pi)+x;
azq=B*Dw*sin(20/180*pi)+y*cosfaiq;
detaiq=sqrt(vxq^2+ vzq^2)-(ri/Dw-0.5)* Dw;
detaeq=sqrt((axq-vxq)^2+(azq-vzq)^2)-(re/Dw-0.5)*Dw;

cosaiq=vzq/((ri/Dw-0.5)*Dw+detaiq);
sinaiq=vxq/((ri/Dw-0.5)*Dw+detaiq);
coseq=(azq-vzq)/((re/Dw-0.5)*Dw+detaeq);
sineq=(axq-vxq)/((re/Dw-0.5)*Dw+detaeq);
gamaiq=Dw*cosaiq/dm;
gameq=Dw*coseq/dm;

rouaiq=4/Dw+2*gamaiq/(Dw*(1-gamaiq))-1/ri;
roueq=4/Dw-2*gameq/(Dw*(1+gameq))-1/re;

betaq=atan(coseq*sineq/(coseq^2+gameq));%外圈控制
nbq=30/pi*dm/Dw*(Ni+Ne)*(1-gamaiq)*(1+gameq)/((1-gamaiq)*cos(acos(coseq)-betaq)+(1+gameq)*cos(acos(cosaiq)-betaq));
nmq=(30/pi*Ni*(1-gamaiq)*cos(acos(coseq)-betaq)-Ne*(1+gameq)*cos(acos(cosaiq)-betaq))/((1-gamaiq)*cos(acos(coseq)-betaq)+(1+gameq)*cos(acos(cosaiq)-betaq));

fcq=2.26e-11*dm*Dw^3*nmq^2;
mgq=4.5e-12*Dw^5*nbq*nmq*sin(betaq);

ki=2.15e5/(sqrt(rouaiq)*sqrt(Ndeta^3));
ke=2.15e5/(sqrt(roueq)*sqrt(Ndeta^3)); %Ndeta 查表

Qiq=ki*sqrt(detaiq^3);
Qeq=ke*sqrt(detaeq^3);



f1=Qiq*sinaiq-Qeq*sineq+2*mgq/Dw*coseq;
f2=Qiq*sinaiq-Qeq*sineq-2*mgq/Dw*sineq+fcq;

[jie1,jie2]=solve(f1,f2)
jievxq(i)=jie1;
jievzq(i)=jie2;
end

请问各位怎么解不出来，，cpu用100%，好久解不出。
请高手指点。