sj1 Function

private elemental function sj1(x)

Uses

FIXME : Add documentation

Arguments

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

Return Value real(kind=kind_rs)


Contents

Source Code

sj1

Source Code

  elemental function sj1 (x)
    use warning_helpers, only: is_zero
    implicit none
    real (kind=kind_rs), intent (in) :: x
    real (kind=kind_rs) :: sj1
    if (is_zero(x)) then
       sj1 = 0.5_kind_rs
    else
       sj1 = bessel_j1(x) / x
    end if
  end function sj1