0

次のプログラムにエラーがあり、どれが問題なのかわかりません。

 IMPLICIT NONE
        DOUBLE PRECISION X(100),W(100),lambdainv,g,lambda,alfac
        INTEGER i, N

        N=100
        alfac=939.D0*2.D22

        CALL GAUSS_L(1.D0,100.D0,100,X,W)

        lambdainv=0.D0

        DO i=1,N
          lambdainv=lambdainv+((W(i)*(X(i))**2)/(alfac+ 
     >              (X(i))**2))*g(X(i)))**2

        END DO

          lambda=lambdainv**(-1)

        WRITE(*,*)'lambda=', lambda

        STOP
        END

c       Funcion g
        DOUBLE PRECISION FUNCTION g(X)

        IMPLICIT NONE  
        DOUBLE PRECISION X, mu, pi

c       Inicializamos las variables
        mu=138.d0
        pi=ATAN(1.D0)*4.D0


        g=(2.d0*sqrt(mu**3))/((sqrt(pi))*(mu**2+X**2))

        RETURN
        END

エラーは次のとおりです。

lambdainv=lambdainv+((W(i)*(X(i))**2)/(alfac+                 
          1

Error: Unclassifiable statement at (1)

ありがとう！！

4

1 に答える 1