double x_uc[4]={0.000,0,0,0};
double y_uc[4]={0.0000,0,0,0};
//DNA
//double X[10] =
{0.000,7.000,11.197,17.910,28.649,45.826,73.301,117.251,187.551,300.00};
//double Y[10] = {272,3008,5854,10287,18398,32204,56751,95727,151346,232604};
//HCG
//double X[10] =
{0.000,0.500,1.000,5.000,10.000,50.000,100.000,500.000,1000.000,5000.000};
//double Y[10] = {1477,1869,1920,4279,7010,28345,55156,245217,415291,1294682};
//HSV-M
//double X[10] = {0.000,0.500,0.645,0.832,1.073,1.385,1.786,2.304,2.973,3.835};
//double Y[10] = {349,4273,5487,7054,9072,11685,15043,19381,24984,32203};
//Trep
double X[10] = {0.000,0.400,0.761,1.450,2.760,5.254,10.001,19.039,36.245,70.000};
double Y[10] = {988,2620,4889,9686,19423,39364,81962,152041,249416,398202};
int n;
/*string s_a;// 用于 a b c 段线
string s_b;// 用于 b c 段线
string s_c;// 用于 c 段线
string s_avarage;// 用于 b c 段线 求均值
*/
for(n=0;n<9;n++)
{
//computing(X,Y ,n);
int k_a,k_b,k_c;
// 插值函数调用,并把结果输入到文件
if(n == 0)
{
cout<<n<<endl;
for(k_a=0;k_a<4;k_a++)
{
x_ua[k_a] = X[k_a+n];
y_ua[k_a] = Y[k_a+n];
}// 给 x、y 赋初值。
ofstream ofs_all("d:\\spline_all.txt");// 写文件流
if(!ofs_all) return 10000;
tangent(y_ua,dy_ua,4);
for(xa_ua = X[0];xa_ua <= X[1];xa_ua += 0.01)
{
s_ua=spline(x_ua,y_ua,dy_ua,xa_ua,4);