//InterpolationHermiteIsometry.cpp  
//埃米爾特等距插值

#include <iostream.h>		//輸入輸出流頭文件
#include "Interpolation.h"	//插值頭文件

void main()
{
	double x, h, t, s;

	double y[10] = 
	{
		0.904837, 0.818731, 0.740818, 0.670320, 0.606531,
		0.548812, 0.496585, 0.449329, 0.406570, 0.367879
	};
	
	valarray<double> yy(y,10);
	valarray<double> dy(10);

	for(int i = 0; i < 10; i ++) dy[i] = -y[i];
	
	x = 0.1;
	h = 0.1;

	cout << " InterpolationHermiteIsometry()" << endl << endl;
	
	t = 0.752;
	cout << " t = " << t;
	s = InterpolationHermiteIsometry(x,h,yy,dy,t);
    cout << "\t s = " << s << endl;
}