program define pdf_weib
    version 4.0

    local x  = "`1'"
    local fx = "`2'"
    local c  = `3'

    quietly{
        replace `fx'=exp(log(`c')+(`c'-1)*log(`x')-(`x')^`c')}

end