Polynomial for Radau-Gauss quadrature
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | n | |||
double precision, | intent(in) | :: | x |
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