program define pdf_logs
    version 4.0

    local x     = "`1'"
    local fx    =  "`2'"
    local alpha = `3'
    local beta  = `4'

    quietly {
        replace `fx' = exp(-(`x'-`alpha')/`beta')
        replace `fx' = `fx'/(`beta'*(1+`fx')^2)}

end