F(1)=xa(1)+y(3)/y(1)+1./cf*(xa(3)-xa(1)**2-1./3.*xa(2)**2-2./3.*xa(2)*y(4)/y(1))
  F(2)=xa(2)+2*y(4)/y(1)+1./cf*(xa(4)-2*xa(1)*xa(2))
  F(3)=xa(3)-1./(1+sigm)*((0.5+0.5*noug/nouf+sigm)*(xa(1)**2+1./3.*xa(2)**2)-noug/nouf*cg*y(3)/y(1)+2./3.*(1+sigm)*xa(2)*y(4)/y(1)-z)
  F(4)=xa(4)-1./(1+sigm)*((1+noug/nouf+2*sigm)*xa(1)*xa(2)-(1-noug/nouf)*g*y(1)-noug/nouf*cg*(2*y(4)/y(1)+y(5)/y(1))+3./8.*Cd*xa(2)**2)
  F(5)=xa(5)-1./(1+sigm)*((2+cg/cf+sigm)*xa(1)*xa(2)+cg/cf*(1-noug/nouf)*g*y(1)-cg*(2*y(4)/y(1)+y(5)/y(1))-cg/cf*3./8.*Cd*xa(2)**2)

系数已知，xguess（5）赋的初值都是0，y（5）已知的数组y=[0.3673,0.0,-68.8183,0,-0.507]
求解xa（5）
我用的isml库函数 call neqnf(fcn,errrel,n,itmax,xguess,x,fnorm)
求解的到错误的答案
请问还有没有库函数可以求解？？或者能否有高手给出解答方案！！谢谢各位！！