Update calcPh.m