Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=8), | intent(in) | :: | x | |||
real(kind=8), | intent(in), | optional | :: | d |
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