Fortran函数返回NaN

时间:2017-10-08 18:07:09

标签: fortran

我试图通过迭代找到功能值。这是我的代码。

   x(3)=g1
   x(4)=1.0
   x(5)=g2
   x(8)=1.0
   x(15)=1.0



    a11=x(10)**2!x(7)**2+x(8)**2+x(9)**2+x(10)**2
    a12=x(7)*x(12)!+x(8)*x(13)+x(9)*x(14)+x(10)*x(15)
    a21=a12
    a22=x(12)**2+x(13)**2+x(14)**2+x(15)**2
    d1=-(x(2)*x(7)+x(3)*x(8)+x(4)*x(9)+x(5)*x(10))
    d2=-(x(2)*x(12)+x(3)*x(13)+x(4)*x(14)+x(5)*x(15))

   det0=a11*a22-a12*a21
   det1=d1*a22-d2*a12
   det2=d2*a11-d1*a21

   dx=det1/det0
   dy=det2/det0

     p1=an+3.
     p2=an+1.
     f(1)=x(2)
     f(2)=x(3)
     f(3)=-p1*x(1)*x(3)+2.*p2*x(2)**2-x(4)
     f(4)=x(5)
     f(5)=-pr*p1*x(1)*x(5)+4.*an*x(2)*x(4)

     do k=1,2
     i(1)=5*k+1
     do j =2,5
      i(j)=i(j-1)+1
     end do
     f(i(1))=x(i(2))
     f(i(2))=x(i(3))
     f(i(3))=-p1*(x(1)*x(i(3))+x(3)*x(i(1)))+4.*p2*x(2)*x(i(2))-x(i(4))
     f(i(4))=x(i(5))
     f(i(5))=-pr*p1*(x(1)*x(i(4))+x(5)*x(i(1)))+4.*an*pr*   (x(2)*x(i(4))+x(4)*x(i(2)))
     end do

我期待一些数值但是一些变量给出NaN值a11,a12,a21,a22,det0,det1,det2。我猜这些功能正在产生问题。

0 个答案:

没有答案