dens_bethe Function

elemental function dens_bethe(x, d)

Arguments

Type IntentOptional Attributes Name
real(kind=8), intent(in) :: x
real(kind=8), intent(in), optional :: d

Return Value real(kind=8)


Source Code

elemental function dens_bethe(x,D)
  real(8),intent(in)          :: x
  real(8),intent(in),optional :: d
  real(8)                     :: dens_bethe,d_
  complex(8)                  :: root,d2
  d_=1.d0;if(present(d))d_=d
  d2=dcmplx(d_,0.d0)
  root=dcmplx((1.d0-1.d0*((x/d_))**2),0.d0)
  root=sqrt(root)
  dens_bethe=(2.d0/(3.141592653589793238d0*d_))*root
end function dens_bethe