//InterpolationAitkenNotIsometry.cpp	
//埃特金不等距逐步法插值 

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

void main()
{
	double eps = FLOATERROR, t, z;
	double x[11] = 
	{
		-1.0, -0.8, -0.65, -0.4, -0.3, 
		 0.0,  0.2,  0.4,   0.6,  0.8, 1.0
	};
	
	double y[11] = 
	{
		0.0384615, 0.0588236, 0.0864865, 0.2, 0.307692,
		1.0,       0.5,       0.2,       0.1, 0.0588236, 0.0384615
	};
    
	valarray<double> yy(y,11), xx(x,11);

	cout << " InterpolationAitkenNotIsometry()" << endl << endl;
	
	t = 0.15;
	z = InterpolationAitkenNotIsometry(xx, yy, t, eps);
	cout << " t = " << t << "\t" << "z = " << z << endl;

	t = 0.05;
	z=InterpolationAitkenNotIsometry(xx, yy, t, eps);
	cout << " t = " << t << "\t" << "z = " << z << endl;
}