三次样条插值的C语言实现

#include
#include
int main（void）{
	double x[19]={0.523.18.017.9528.6539.6250.6578104.6156.6208.6260.7312.5364.4416.3468494507520};
    double y[19]={5.287949.413.8420.224.928.4431.13536.536.634.631.026.3420.914.87.83.71.50.2};
	double t[12]={24612163060110180280400515};
    double h[19]u[19]l[19]g[19];
	double M[19]n[19];
	void tri（double a[]double d[]double c[]double b[]double M[]）;
	int ij;
	double result;
	for（i=1;i<=18;i++）
		h[i]=x[i]-x[i-1];
	for（i=1;i<=17;i++）
		u[i]=h[i]/（h[i]+h[i+1]）;
	for（i=1;i<=17;i++）
		l[i]=1-u[i];
	for（i=1;i<=17;i++）
		g[i]=6.0/（h[i]+h[i+1]）*（（y[i+1]-y[i]）/h[i+1]-（y[i]-y[i-1]）/h[i]）;
	g[0]=6.0/h[1]*（（y[1]-y[0]）/h[1]