前提
fortranで次の非線形方程式 f(x)=tanh(x)+0.2x+0.3= 0 の数値解をNewton-Raphson 法を用いて求めようとしたところ、コンパイルの段階で次のようなエラーが出てしまいました。何か解決策はありますでしょうか。
発生している問題・エラーメッセージ
gfortran keisan.f90 f951.exe:Error:Unexpected end of file in 'keisan.f90'
該当のソースコード
fortran
1PROGRAM keisan 2IMPLICIT REAL(A-H,O-Z) 3EXTERNAL FUN1,DIVFUN1 4DIMENSION X_D(50,3) 5 6DATA EPS/1.E-7/ 7 8AA=0.2 9BB=0.3 10 11WRITE(*,*) ' ' 12write(*,*) ' N X2 F(X2) REL_ERROR THEO_ERROR' 13 14AA=0.2 15BB=0.3 16XX=1.1 17X1=XX 18DO N=1,100 19 X2=X1-FUN1(X1,AA,BB)/DIVFUN1(X1,AA) 20 21 EPSD=ABS(X2-X1) 22 EPSX=ABS(EPSD/X1) 23 24 WRITE(*,120) N,X2,FUN1(X2,AA,BB),EPSX 25120 FORMAT(' ',I3,5E17.7) 26 27 IF(EPSX.LT.EPS) EXIT 28 X1=X2 29ENDDO 30WRITE(*,*) ' ' 31 32END 33 34REAL FUNCTION FUN1(X,AA,BB) 35 36IMPLICIT REAL(A-H,O-Z) 37 38FUN1=tanh(X)+AA*X+BB 39 40RETURN 41END 42 43REAL FUNCTION DIVFUN1(X,AA) 44 45IMPLICIT REAL(A-H,O-Z) 46 47DIVFUN1=(1/cosh(X))**2+AA 48 49RETURN 50END 51 52SUBROUTINE BISECTION(NSTEP,EPS,X1,X2,AA,BB,XD,NEND) 53
試したこと
拡張子をf95に変えてみたり、文字コードを色々変えてみたりしましたが上手くいきませんでした。
補足情報(FW/ツールのバージョンなど)
PCはWindowsを使用しています。
回答1件
あなたの回答
tips
プレビュー