radau_p Function

private elemental function radau_p(n, x)

Polynomial for Radau-Gauss quadrature

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: n
double precision, intent(in) :: x

Return Value double precision


Contents

Source Code


Source Code

  elemental function radau_p (n, x)
    implicit none
    integer, intent (in) :: n
    double precision, intent (in) :: x
    double precision :: radau_p

    radau_p = (legendre_p(n-1,x) + legendre_p(n,x) ) &
         / (x + dble(1.0))

  end function radau_p