exponentialrnd Function

function exponentialrnd(mean) result(c)

Arguments

Type IntentOptional Attributes Name
real(kind=8) :: mean

Return Value real(kind=8)


Calls

proc~~exponentialrnd~~CallsGraph proc~exponentialrnd exponentialrnd mersenne mersenne proc~exponentialrnd->mersenne

Source Code

function exponentialrnd(mean) result(c)
  real(8) :: mean,c
  if (mean <= 0d0) then
     write(*,*) "EXPONENTIALRND: mean must be positive"
  else
     c=-mean*log(mersenne())
  end if
end function exponentialrnd