prodestLP function

Estimate productivity - Levinsohn-Petrin method