estimatePhi function

estimatePhi