integration.f90

Go to the documentation of this file.

MODULE integration

CONTAINS

! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
!  This subroutine calculates integral of a function y(x)
!  from X=0
SUBROUTINE integral_value(N,X,Y,INTY)

    use itm_types

    IMPLICIT NONE

    INTEGER :: n                                          ! number of radial points (input)
    INTEGER :: i

    REAL (R8) :: x(n), &                                  ! argument array (input)
         y(n), &                                  ! function array (input)
         inty(n)                                  ! function integral array (output)


    inty(1)=y(1)*x(1)/2.0_r8
    DO i=2,n
       inty(i)=inty(i-1)+(y(i)+y(i-1))*(x(i)-x(i-1))/2.0_r8
    ENDDO

    RETURN
END SUBROUTINE integral_value
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  
! + + + + + + + + + + + + + + + + + + + + + + + + + + + +  

END MODULE integration