高斯-勒让德求积工具

chiangtp

kerb 发表于 2017-9-17 21:38
第二版
第二版会生成一个例子程序，例子是sin(x)在区间[0.0,3.1415926]的积分，区间分成5段（程序中整形变 ...

謝謝分享

雞蛋裡挑個骨頭, 還請包含

[Fortran] 纯文本查看 复制代码

PROGRAM test
  INTEGER, PARAMETER :: dp=KIND(1.0D0)
  REAL(KIND=dp) :: x=1.2, y=1.2D0     ! (KIND=8),  literal "8" is compiler dependent
  WRITE(*,*) x, y                     ! x, default precision only
END PROGRAM test

WRITE(*,'(4X,A)',ADVANCE='NO') 'A = ' ! [), 'A = '] ---> [) 'A = '],  no ","
READ(*,*) A                           ! (F15.8), formatted READ, dangerous

FUNCTION GUALEGINTG12(YOUR_FUNC, OTHERS, A, B)
  IMPLICIT NONE
  INTEGER, PARAMETER :: dp=KIND(1.0D0)
  REAL(KIND=dp)             :: GUALEGINTG12
  REAL(KIND=dp), EXTERNAL   :: YOUR_FUNC
  REAL(KIND=dp), INTENT(IN) :: OTHERS, A, B