曲率计算代码