program define pdf_par
    version 4.0

    local x  = "`1'"
    local fx = "`2'"
    local a  = `3'
    local k  = 1

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

end