*!version 1.0.0  6-10-98  Joseph Hilbe   STB-45 dm60
* digamma function: immediate command
* Ex: digami (numnber) ; the digamma value, if valid, will follow

program define digami
  version 5.0
  local dg = (lngamma(`1'+.0001)-lngamma(`1'-.0001))/.0002
  noi di %8.4f `dg'
end