PT计算公式C程序