In reply to Qui Phan.

It should be straightforward (but less efficient) to use PROC FCMP to define the functions. Each call to the CDF will have to recompute the cumulative PDF from 0 to x. For the quantile, you can just accumulate the PDF until the cumulative sum exceeds p.