C类C第277-278页：数值数学与计算，CHENEY/KINCAID，1985C类C文件：SPL3.FORC类等距点SIN（X）的C自然三次样条函数（SPL3，ZSPL3）C类参数（N=10）尺寸T（N）、Y（N）和U（N）H=1.6875/真实（N-1）DO 2 I=1，NT（I）=真实（I-1）*HY（I）=SIN（T（I））2继续调用ZSPL3（N，T，Y，W，B，U，V，Z）DO 3 I=1,4*N-3X=真实（I-1）*H*0.25D=SIN（X）-SPL3（N、T、Y、Z、X）打印4、I、X、D3继续4格式（2X、I5、F22.14、E10.3）停止结束子程序ZSPL3（N、T、Y、H、B、U、V、Z）尺寸T（N）、Y（N）和H（N）DO 2 I=1，N-1H（I）=T（I+1）-T（I）B（I）=（Y（I+1）-Y（I））/H（I）2继续U（2）=2.0*（H（1）+H（2））V（2）=6.0*（B（2）-B（1））DO 3 I=3，N-1U（I）=2.0*（H（I）+H（I-1））-H（I-1V（I）=6.0*（B（I）-B（I-1））-H（I-13继续Z（N）=0.0DO 4 I=N-1,2，-1Z（I）=（V（I）-H（I）*Z（I+1））/U（I）4继续Z（1）=0.0返回结束功能SPL3（N、T、Y、Z、X）尺寸T（N）、Y（N）和Z（N）DO 2 I=N-1,1，-1差异=X-T（I）如果（DIFF.GE.0.0）转到32继续I=13 H=T（I+1）-T（I）B=（Y（I+1）-Y（I））/H-H*（Z（I+1P=0.5*Z（I）+DIFF*（Z（I+1）-Z（I））/（6.0*H）P=B+DIFF*PSPL3=Y（I）+差异*P返回结束