program define pdf_Z
    version 4.0

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

    quietly {
        replace `fx' = exp(-(`x'*`x')/2)/sqrt(2*_pi)
    }

end