#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]