くだすれFORTRAN（超初心者用）

くだすれFORTRAN（超初心者用） at TECH

368:2/2
06/06/02 02:32:02
REAL FUNCTION sub_trapez(x, h)
IMPLICIT NONE
REAL, INTENT(IN) :: x, h
sub_trapez = 0.5 * h * ( func(x) + func(x + h) )
RETURN
CONTAINS
REAL FUNCTION func(x)
IMPLICIT NONE
REAL, INTENT(IN) :: x
func = EXP( -x**2 ) / SQRT(pi)
RETURN
END FUNCTION func
END FUNCTION sub_trapez
END MODULE mod_integ
!
PROGRAM oppai
USE mod_integ
IMPLICIT NONE
INTEGER :: i
REAL :: s, s1, s2
pi = 4.0 * ATAN(1.0) ! arctan(1.0) = pi / 4; tan(pi/4) = 1.0
! [-1..1]
s1 = trapez(100, -1.0, 0.0)
! [-2..2]
s2 = trapez(200, -2.0, 2.0)
! [-inf..+inf]
s = trapez(1000, -10.0, 10.0)
PRINT *, s, s1, s2
STOP
END PROGRAM oppai
>>367-368

次ページ

続きを表示

1を表示