1. spl1插值函数：spl1(x0,y0,x,fx,dfx,s,n1,n2,ndim1,ndim2)

2. 高斯积分 gauleg(0.d0,1.d0,arc_gauss,warc_gauss,narc_gauss)

其中：arc_gauss是0-1之间的分段间隔，warc_gauss是l上的分段间隔。

3. 切线和法线的确定

1) 两点a(y,z),b( , )，切线方向为： ，

即： ，

故：已知a、b（切线方向点）两点和他们距离l，则

2） 法线方向为：（-dz，dy）