pcens_cdf function

Compute primary event censored CDF