cauchyrnd Function

function cauchyrnd(median, scale) result(ans)

Arguments

Type IntentOptional Attributes Name
real(kind=8) :: median
real(kind=8) :: scale

Return Value real(kind=8)


Calls

proc~~cauchyrnd~~CallsGraph proc~cauchyrnd cauchyrnd mersenne mersenne proc~cauchyrnd->mersenne

Source Code

function cauchyrnd(median,scale) result(ans)
  real(8) :: ans,median,scale,p
  if (scale <= 0d0) then
     write(*,*) "CAUCHYRND: Scale parameter must be positive"
  end if
  ans = median + scale*tan(PI*(mersenne() - 0.5))
end function cauchyrnd