inverse_gammarnd Function

function inverse_gammarnd(shape, scale) result(ans)

Arguments

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

Return Value real(kind=8)


Calls

proc~~inverse_gammarnd~~CallsGraph proc~inverse_gammarnd inverse_gammarnd gammarnd gammarnd proc~inverse_gammarnd->gammarnd

Source Code

function inverse_gammarnd(shape,scale) result(ans)
  real(8) :: shape,scale,ans
  ! If X is gamma(shape, scale) then
  ! 1/Y is inverse gamma(shape, 1/scale)
  ans= 1d0/gammarnd(shape,1d0/scale)
end function inverse_gammarnd