计算方法:拟合曲线代码