plend function

Compute Product Limit Estimate for Non-detects