bessely1 Function

public function bessely1(x)

Bessel function of the second kind of order one .

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x

Return Value real(kind=wp)


Calls

proc~~bessely1~~CallsGraph proc~bessely1 bessely1 proc~caljy1 caljy1 proc~bessely1->proc~caljy1 proc~nan nan proc~bessely1->proc~nan proc~ninf ninf proc~bessely1->proc~ninf

Source Code

  real(wp) function bessely1(x)
    !! Bessel function of the second kind of order one \(Y_1(x)\).
    !
    !! \(\lbrace x \in \mathbb{R} \mid x \gt 0 \rbrace\)

    real(wp), intent(in) :: x

    if (x < 0) then
      bessely1 = nan()
    else
      call caljy1(x, bessely1, 1)
      if (bessely1 <= -huge(0.0_wp)) then
        bessely1 = ninf()
      end if
    end if
  end function bessely1