亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频

? 歡迎來到蟲蟲下載站！ | ?? 資源下載 ?? 資源專輯 ?? 關于我們
? 蟲蟲下載站

?? interpolation.inl

?? 主要是關于一些數值計算中所要用到的頭文件,希望對大家有所幫助
?? INL
?? 第 1 頁 / 共 2 頁
字號: 
		{
			u[1]  = (y[kk] - y[kk - 1]) / (x[kk] - x[kk - 1]);
			
			if(kk == (n - 3))
			{
				u[3] = (y[n - 1] - y[n - 2]) / (x[n - 1]- x[n - 2]);
				u[4] = 2.0 * u[3] - u[2];
				
				if(n == 4) u[0] = 2.0 * u[1] - u[2];
				
				else u[0] = (y[kk - 1] - y[kk - 2]) / (x[kk - 1] - x[kk - 2]);
			}
			else
			{
				u[3] = 2.0 * u[2] - u[1];
				u[4] = 2.0 * u[3] - u[2];
				u[0] = (y[kk - 1] - y[kk - 2]) / (x[kk - 1] - x[kk - 2]);
			}
		}
		else
		{
			u[1] = (y[kk] - y[kk - 1]) / (x[kk] - x[kk - 1]);
			u[0] = (y[kk - 1] - y[kk - 2]) / (x[kk - 1] - x[kk - 2]);
			u[3] = (y[kk + 2] - y[kk + 1]) / (x[kk + 2] - x[kk + 1]);
			u[4] = (y[kk + 3] - y[kk + 2]) / (x[kk + 3] - x[kk + 2]);
		}
	}
	
	s[0] =  Abs(u[3] - u[2]);
	s[1] =  Abs(u[0] - u[1]);
	
	if(FloatEqual(s[0],0) && FloatEqual(s[1],0)) p = (u[1] + u[2]) / 2.0;
	
	else p = (s[0] * u[1] + s[1] * u[2]) / (s[0] + s[1]);

	s[0] = Abs(u[3] - u[4]);
	s[1] = Abs(u[2] - u[1]);

	if(FloatEqual(s[0],0) && FloatEqual(s[1],0)) q = (u[2] + u[3]) / 2.0;
	
	else q = (s[0] * u[2] + s[1] * u[3]) / (s[0] + s[1]);
	
	s[0] = y[kk];
	s[1] = p;
	s[3] = x[kk + 1] - x[kk];
	s[2] = (3.0 * u[2] - 2.0 * p - q) / s[3];
	s[3] = (q + p - 2.0 * u[2]) / (s[3] * s[3]);
	
	if(k < 0)
	{ 
		p = t - x[kk];
		s[4] = s[0] + s[1] * p + s[2] * p * p + s[3] * p * p * p;
	}

END: ;
}

//光滑等距插值
template <class _Ty>
void InterpolationSmoothIsometry(_Ty x0, _Ty h, 
					valarray<_Ty>& y, int k, _Ty t, valarray<_Ty>& s)
{
	int kk, m, l;
    _Ty u[5], p, q;

	int n =  y.size();	//給定等距結點的個數
    
	for(m=0; m<5; m++) s[m] = 0.0;
    
	if(n < 1) goto END;

    if(n == 1) 
	{
		s[0] = y[0];
		s[4] = y[0]; 
		goto END;
	}

    if(n == 2)
	{
		s[0] = y[0];
		s[1] = (y[1] - y[0]) / h;
        
		if(k < 0) s[4] = (y[1] * (t - x0) - y[0] * (t - x0 - h)) / h;
        
		goto END;
	}
    
	if(k < 0)
	{
		if(t <= x0 + h) kk = 0;
		else 
			if(t >= x0 + (n - 1) * h) kk = n - 2;
			else
			{
				kk = 1;
				m = n;

				while(((kk - m) != 1) && ((kk - m) != -1))
				{
					l = (kk + m) / 2;
                
					if(t < x0 + (l - 1) * h) m = l;
					else kk = l;
				}
            
				kk = kk - 1;
			}
	}
	else kk = k;
  
	if(kk >= n - 1) kk = n - 2;
    
	u[2] = (y[kk + 1]- y [kk]) / h;
	
	if(n == 3)
	{
		if(kk == 0)
		{
			u[3] = (y[2] - y[1]) / h;
            u[4] = 2.0 * u[3] - u[2];
            u[1] = 2.0 * u[2] - u[3];
            u[0] = 2.0 * u[1] - u[2];
		}
        
		else
		{ 
			u[1] = (y[1] - y[0]) / h;
            u[0] = 2.0 * u[1] - u[2];
            u[3] = 2.0 * u[2] - u[1];
            u[4] = 2.0 * u[3] - u[2];
		}
	}  
	else
	{
		if(kk <= 1)
		{
			u[3] = (y[kk + 2] -y[kk + 1]) / h;
            
			if(kk == 1)
			{
				u[1] = (y[1] - y[0]) / h;
                u[0] = 2.0 * u[1] - u[2];
                
				if(n == 4) u[4] = 2.0 * u[3] - u[2];
				else u[4] = (y[4] - y[3]) / h;
			}
            else
			{ 
				u[1] = 2.0 *u[2] - u[3];
                u[0] = 2.0 *u[1] - u[2];
                u[4] = (y[3] - y[2]) / h;
			}
		}
		else if(kk >= (n - 3))
		{
			u[1] = (y[kk] - y[kk - 1]) / h;
            
			if(kk == (n - 3))
			{
				u[3] = (y[n - 1] - y[n - 2]) / h;
                u[4] = 2.0 * u[3] - u[2];
                
				if(n == 4) u[0] = 2.0 * u[1] - u[2];
                else u[0] = (y[kk - 1] - y[kk - 2]) / h;
			}   
			else
			{
				u[3] = 2.0 * u[2] - u[1];
                u[4] = 2.0 * u[3] - u[2];
                u[0] = (y[kk - 1] - y[kk - 2]) / h;
			}
		}
		else
		{
			u[1] = (y[kk] - y[kk - 1]) / h;
            u[0] = (y[kk - 1] - y[kk - 2]) / h;
            u[3] = (y[kk + 2] - y[kk + 1]) / h;
            u[4] = (y[kk + 3] - y[kk + 2]) / h;
		}
	}
     
	s[0] = Abs(u[3] - u[2]);
    s[1] = Abs(u[0] - u[1]);
    
	if(FloatEqual(s[0],0) && FloatEqual(s[1],0))
		p = (u[1] + u[2]) / 2.0;
	else
		p = (s[0] * u[1] + s[1] * u[2]) / (s[0] + s[1]);
    
	s[0] = Abs(u[3] - u[4]);
    s[1] = Abs(u[2] - u[1]);
    
	if(FloatEqual(s[0],0) && FloatEqual(s[1],0))
		q = (u[2] + u[3]) / 2.0;
	else
		q = (s[0] * u[2] + s[1] * u[3]) / (s[0] + s[1]);
	  
	s[0] = y[kk];
    s[1] = p;
    s[3] = h;
    s[2] = (3.0 * u[2] - 2.0 * p - q) / s[3];
    s[3] = (q+p - 2.0 * u[2]) / (s[3] * s[3]);
    
	if(k < 0)
	{
		p = t - (x0 + kk * h);
        s[4] = s[0] + s[1] * p + s[2] * p * p + s[3] * p * p * p;
	}
    
END: ;
}

//第一種邊界條件的三次樣條函數插值、微商與積分
template <class _Ty>
_Ty Interpolation3Spooling1stBoundary(valarray<_Ty>& x, 
		valarray<_Ty>& y, valarray<_Ty>& dy, valarray<_Ty>& ddy, 
			valarray<_Ty>& t, valarray<_Ty>& z, valarray<_Ty>& dz, 
													valarray<_Ty>& ddz)
{
	int i, j;
    _Ty h0, h1, alpha, beta, g;

	int n =  y.size();	//數組y的長度(元素個數)，給定結點個數
	int m =  t.size();	//指定插值點的個數
    
	valarray<_Ty> s(n);

    s[0] = dy[0];
	dy[0] = 0.0;
    h0 = x[1] - x[0];

    for(j = 1; j < n - 1; j ++)
	{
		h1 = x[j + 1] - x[j];
        alpha = h0 / (h0 + h1);
        
		beta = (1.0 - alpha) * (y[j] - y[j - 1]) / h0;
        beta = 3.0 * (beta + alpha * (y[j + 1] - y[j]) / h1);
        
		dy[j] = -alpha / (2.0 + (1.0 - alpha) * dy[j - 1]);
        s[j] = (beta - (1.0 - alpha) * s[j - 1]);
        s[j] = s[j] / (2.0 + (1.0 - alpha) * dy[j - 1]);
        
		h0 = h1;
	}
    
	for(j = n - 2; j >= 0; j --)
	{
		dy[j] = dy[j] * dy[j + 1] + s[j];
	}

    for(j = 0; j < n - 1; j ++) 
	{
		s[j] = x[j + 1] - x[j];
	}
    
	for(j = 0 ; j < n - 1; j ++)
	{
		h1 = s[j] * s[j];
        ddy[j] = 6.0 * (y[j + 1] - y[j]) / h1 - 2.0 * (2.0 * dy[j] + dy[j + 1]) / s[j];
	}

    h1 = s[n - 2] * s[n - 2];
    ddy[n - 1] = 6.0 * (y[n - 2] - y[n - 1]) / h1 + 2.0 * (2.0 * dy[n - 1] 
		         + dy[n - 2]) / s[n - 2];
    g = 0.0;
    
	for(i = 0;i < n - 1; i ++)
	{
		h1 = 0.5 * s[i] * (y[i] + y[i + 1]);
        h1 = h1 - s[i] * s[i] * s[i] * (ddy[i] + ddy[i + 1]) / 24.0;
        g = g + h1;
	}
    
	for(j = 0; j <= m - 1; j ++)
	{
		if(t[j] >= x[n - 1]) i = n - 2;
        else
		{
			i = 0;
            while(t[j] > x[i + 1]) i ++;
		}
        
		h1 = (x[i + 1] - t[j]) / s[i];
        h0 = h1 * h1;
        
		z[j] = (3.0 * h0 - 2.0 * h0 * h1) * y[i];
        z[j] = z[j] + s[i] * (h0 - h0 * h1) * dy[i];
        
		dz[j] = 6.0 * (h0 - h1) * y[i] / s[i];
        dz[j] = dz[j] + (3.0 * h0 - 2.0 * h1) * dy[i];
        
		ddz[j] = (6.0 - 12.0 * h1) * y[i] / (s[i] * s[i]);
        ddz[j] = ddz[j] + (2.0 - 6.0 * h1) * dy[i] / s[i];
        
		h1 = (t[j] - x[i]) / s[i];
        h0 = h1 * h1;
        
		z[j] = z[j] + (3.0 * h0 -2.0 * h0 * h1) * y[i + 1];
        z[j] = z[j] - s[i] * (h0 - h0 * h1) * dy[i + 1];
        
		dz[j] = dz[j] - 6.0 * (h0 - h1) * y[i + 1] / s[i];
        dz[j] = dz[j] + (3.0 * h0 - 2.0 * h1) * dy[i + 1];
        
		ddz[j] = ddz[j] + (6.0 - 12.0 * h1) * y[i + 1] / (s[i] * s[i]);
        ddz[j] = ddz[j] - (2.0 - 6.0 * h1) * dy[i + 1] / s[i];
	}
    
	return(g);
}

//第二種邊界條件的三次樣條函數插值、微商與積分
template <class _Ty>
_Ty Interpolation3Spooling2ndBoundary(valarray<_Ty>& x, 
		valarray<_Ty>& y, valarray<_Ty>& dy, valarray<_Ty>& ddy, 
			valarray<_Ty>& t, valarray<_Ty>& z, valarray<_Ty>& dz, 
													valarray<_Ty>& ddz)
{
	int i, j;
    _Ty h0, h1, alpha, beta, g;

	int n =  y.size();	//數組y的長度(元素個數)，給定結點個數
	int m =  t.size();	//指定插值點的個數
    
	valarray<_Ty> s(n);
    
	dy[0] = -0.5;
    h0 = x[1] - x[0];
    s[0] = 3.0 * (y[1] - y[0]) / (2.0 * h0) - ddy[0] * h0 / 4.0;
    
	for(j = 1; j < n - 1; j ++)
	{
		h1 = x[j + 1] - x[j];
        alpha = h0 / (h0 + h1);
        
		beta = (1.0 - alpha) * (y[j] - y[j - 1]) / h0;
        beta = 3.0 * (beta + alpha * (y[j + 1] - y[j]) / h1);
        
		dy[j] = -alpha / (2.0 + (1.0 - alpha) * dy[j - 1]);
        s[j] = (beta - (1.0 - alpha) * s[j - 1]);
        s[j] = s[j] / (2.0 + (1.0 - alpha) * dy[j - 1]);
        
		h0 = h1;
	}
    
	dy[n - 1] = (3.0 * (y[n - 1] - y[n - 2]) / h1 + ddy[n - 1] * h1 
		        / 2.0 - s[n - 2]) / (2.0 + dy[n - 2]);
    
	for(j = n - 2; j >= 0; j --)
	{
		dy[j] = dy[j] * dy[j + 1 ] + s[j];
	}
    
	for(j =0; j < n - 1; j ++)
	{
		s[j] = x[j + 1] - x[j];
	}

    for(j =0; j < n - 1; j ++)
	{
		h1 = s[j] * s[j];
        ddy[j] = 6.0 * (y[j + 1] - y [j]) / h1 - 2.0 * (2.0 * dy[j] + dy[j + 1]) / s[j];
	}
    
	h1 = s[n - 2] * s[n - 2];
    ddy[n - 1]= 6.0 * (y[n - 2] - y[n - 1]) / h1 + 2.0 * (2.0 * dy[n - 1] 
		        + dy[n - 2]) / s[n - 2];
    g = 0.0;
    
	for(i = 0; i < n - 1; i ++)
	{
		h1 = 0.5 * s[i] * (y[i] + y[i + 1]);
        h1 = h1 - s[i] * s[i] * s[i] * (ddy[i] + ddy[i + 1]) / 24.0;
        g = g + h1;
	}
    
	for(j = 0; j <= m - 1; j ++)
	{
		if(t[j] >= x[n - 1]) i = n - 2;
        else
		{
			i = 0;
            while(t[j] > x[i + 1]) i = i + 1;
		}
        
		h1 = (x[i + 1] - t[j]) / s[i];
        h0 = h1 * h1;
        
		z[j] = (3.0 * h0 - 2.0 * h0 * h1) * y[i];
        z[j] = z[j] + s[i] * (h0 - h0 * h1) * dy[i];
        
		dz[j] = 6.0 * (h0 - h1) * y[i] / s[i];
        dz[j] = dz[j] + (3.0 * h0 - 2.0 * h1) * dy[i];
        
		ddz[j] = (6.0 - 12.0 * h1) * y[i] / (s[i] * s[i]);
        ddz[j] = ddz[j] + (2.0 - 6.0 * h1) * dy[i] / s[i];
        
		h1 = (t[j] - x[i]) / s[i];
        h0 = h1 * h1;
        
		z[j] = z[j] + (3.0 * h0 - 2.0 * h0 * h1) * y[i + 1];
        z[j] = z[j] - s[i] * (h0 - h0 * h1)* dy[i + 1];
        
		dz[j] = dz[j] - 6.0 * (h0 - h1) * y[i + 1] / s[i];
        dz[j] = dz[j] + (3.0 * h0 - 2.0 * h1) * dy[i + 1];
        
		ddz[j] = ddz[j] + (6.0 - 12.0 * h1) * y[i + 1] / (s[i] * s[i]);
        ddz[j] = ddz[j] - (2.0 - 6.0 * h1) * dy[i + 1] / s[i];
	}
    
    return(g);
}

//第三種邊界條件的三次樣條函數插值、微商與積分
template <class _Ty>
_Ty Interpolation3Spooling3thBoundary(valarray<_Ty>& x, 
		valarray<_Ty>& y, valarray<_Ty>& dy, valarray<_Ty>& ddy, 
			valarray<_Ty>& t, valarray<_Ty>& z, valarray<_Ty>& dz, 
													valarray<_Ty>& ddz)
{
	int i, j;
    _Ty h0, y0, h1, y1, alpha, beta, u, g;

	int n =  y.size();	//數組y的長度(元素個數)，給定結點個數
	int m =  t.size();	//指定插值點的個數
    
	valarray<_Ty> s(n);
    
	h0 = x[n -1 ] - x[n - 2];
    y0 = y[n - 1] - y[n - 2];
    
	dy[0] = 0.0;
	ddy[0] = 0.0;
	ddy[n - 1] = 0.0;

    s[0] = 1.0;
	s[n - 1] = 1.0;
    
	for(j = 1; j < n; j ++)
	{
		h1 = h0;
		y1 = y0;
        h0 = x[j] - x[j - 1];
        y0 = y[j] - y[j - 1];
        
		alpha = h1 / (h1 + h0);
        beta = 3.0 * ((1.0 - alpha) * y1 / h1 + alpha * y0 / h0);
        
		if(j < n - 1)
		{
			u = 2.0 + (1.0 - alpha) * dy[j - 1];
            dy[j] = -alpha / u;
            s[j] = (alpha - 1.0) * s[j - 1] / u;
            ddy[j] = (beta - (1.0 - alpha) * ddy[j - 1]) / u;
		}
	}
    
	for(j = n - 2; j >= 1; j--)
	{
		s[j] = dy[j] * s[j + 1] + s[j];
        ddy[j] = dy[j] * ddy[j + 1] + ddy[j];
	}
    
	dy[n-2] = (beta - alpha * ddy[1] - (1.0 - alpha) * ddy[n - 2])
		      / (alpha * s[1] + (1.0 - alpha)* s[n - 2] + 2.0);
    
	for(j = 2; j < n; j ++)
	{
        dy[j-2]=s[j-1]*dy[n-2]+ddy[j-1];
	}
	
	dy[n - 1] = dy[0];
    
	for(j = 0; j < n - 1; j++)
	{
		s[j] = x[j + 1] - x[j];
	}

    for(j = 0; j < n - 1; j++)
	{
		h1 = s[j] * s[j];
        ddy[j] = 6.0 * (y[j + 1] - y[j]) / h1 - 2.0 
			     * (2.0 * dy[j] + dy[j + 1]) / s[j];
	}

    h1 = s[n - 2] * s[n - 2];
    ddy[n - 1] = 6.0 * (y[n - 2] - y[n - 1]) / h1 
		         + 2.0 * (2.0 * dy[n - 1] + dy[n - 2]) / s[n - 2];
    g = 0.0;

	for(i = 0; i < n - 1; i++)
    {
		h1 = 0.5 * s[i] * (y[i] + y[i + 1]);
        h1 = h1 - s[i] * s[i] * s[i] * (ddy[i] + ddy[i + 1]) / 24.0;
        g = g + h1;
	}
    
	for(j = 0; j <= m - 1; j++)
	{
		h0 = t[j];
        
		while(h0 >= x[n - 1])
		{
			h0 = h0 - (x[n - 1] - x[0]);
		}
        while(h0 < x[0])
		{
			h0 = h0 + (x[n - 1] - x[0]);
		}
       
		i = 0;
        
		while(h0 > x[i + 1]) i++;
        
		u = h0;
        h1 = (x[i + 1] -u ) / s[i];
        h0 = h1 * h1;
        
		z[j] = (3.0 * h0 - 2.0 * h0 * h1) * y[i];
        z[j] = z[j] + s[i] * (h0 - h0 * h1) * dy[i];
         
		dz[j] = 6.0 * (h0 - h1) * y[i] / s[i];
        dz[j] = dz[j] + (3.0 * h0 - 2.0 * h1) * dy[i];
        
		ddz[j] = (6.0 - 12.0 * h1) * y[i] / (s[i] * s[i]);
        ddz[j] = ddz[j] + (2.0 - 6.0 * h1) * dy[i] / s[i];
        
		h1 = (u - x[i]) / s[i];
        h0 = h1 * h1;
        
		z[j] = z[j] + (3.0 * h0 - 2.0 * h0 * h1) * y[i + 1];
        z[j] = z[j] - s[i] * (h0 - h0 * h1) * dy[i + 1];
        
		dz[j] = dz[j] - 6.0 * (h0 - h1) * y[i + 1] / s[i];
        dz[j] = dz[j] + (3.0 * h0 - 2.0 * h1) * dy[i + 1];
        
		ddz[j] = ddz[j] + (6.0 - 12.0 * h1) * y[i + 1] / (s[i] * s[i]);
        ddz[j] = ddz[j] - (2.0 - 6.0 * h1) * dy[i + 1] / s[i];
	}
    
	return(g);
}

//二元三點插值
template <class _Ty>
_Ty Interpolation2Variable3Points(valarray<_Ty>& x, valarray<_Ty>& y, 
											matrix<_Ty> z, _Ty u, _Ty v)
{
	int nn(3),mm,ip,iq,i,j,k,l;
    _Ty b[3],h,w;

	int n =  x.size();	//給定結點X方向上的坐標個數
	int m =  y.size();	//給定結點Y方向上的坐標個數
    
	if(n < 4)
	{
		ip = 0;
		nn = n;
	}
    else if(u <= x[1]) ip = 0;
    else if(u >= x[n - 2]) ip = n - 3;
    else
	{
		i = 1;
		j = n;
        
		while( ((i - j) != 1) && ((i - j) != -1) )
		{
			l = (i + j) / 2;
            
			if(u < x[l - 1]) j = l;
            else i = l;
		}
        
		if(Abs(u - x[i - 1]) < Abs(u - x[j - 1])) ip = i - 2;
        else ip = i-1;
      }

    mm = 3;
    
	if(m < 4)
	{
		iq = 0;
		mm = m;
	}
    else if(v <= y[1]) iq = 0;
    else if(v >= y[m - 2]) iq = m - 3;
    else
	{
		i = 1;
		j = m;
        
		while( ((i - j) != 1) && ((i - j) != -1) )
		{
			l = (i + j) / 2;
            if(v < y[l - 1]) j = l;
            else i = l;
		}
        
		if(Abs(v - y[i - 1]) < Abs(v - y[j - 1])) iq = i - 2;
        else iq = i - 1;
	}
    
	for(i = 0; i < nn; i ++)
	{
		b[i] = 0.0;

        for(j = 0; j < mm; j ++)
		{
			h = z(ip + i,iq +j);

            for(k = 0; k < mm; k ++)
			{
				if(k != j)
					h = h * (v - y[iq + k]) / (y[iq + j] - y[iq + k]);
			}
			
			b[i] = b[i] + h;
		}
	}
    
	w = 0.0;
    
	for(i = 0; i < nn; i ++)
	{
		h = b[i];
        
		for(j = 0; j < nn; j ++)
		{
			if(j != i)
				h = h * (u - x[ip + j]) / (x[ip + i] - x[ip + j]);
		}
		
		w = w + h;
	}
    
	return(w);
}


//二元全區間插值
template <class _Ty>
_Ty Interpolation2VariableWholeInterval(valarray<_Ty>& x, 
						valarray<_Ty>& y, matrix<_Ty> z, _Ty u, _Ty v)
{
	
	int ip, ipp, i, j, l, iq, iqq, k;
    _Ty h, w;
	valarray<_Ty> b(10);

	int n =  x.size();	//給定結點X方向上的坐標個數
	int m =  y.size();	//給定結點Y方向上的坐標個數

	if(u<x[0]||FloatEqual(u,x[0]))
	{
		ip = 1;
		ipp = 4;
	}
    else if(u>x[n-1]||FloatEqual(u,x[n-1]))
	{
		ip = n - 3;
		ipp = n;
	}
    else
	{
		i = 1;
		j = n;
        
		while(((i - j) != 1) && ((i - j) != -1))
		{
			l = (i + j) / 2;
            
			if(u < x[l - 1]) j = l;
            else i = l;
		}
        
		ip = i - 3; 
		ipp = i + 4;
	}
    
	if(ip < 1) ip = 1;
    if(ipp > n) ipp = n;
    if(v < y[0] || FloatEqual(v, y[0]))
	{ 
		iq = 1;
		iqq = 4; 
	}
    else if(v > y[m - 1] || FloatEqual(v, y[m - 1]))
	{
		iq = m - 3;
		iqq = m;
	}
    else
	{
		i = 1;
		j =m;
        
		while(((i - j) != 1) && ((i - j) != -1))
		{
			l = (i + j) / 2;
            
			if(v <y [l - 1]) j = l;
            else i = l;
		}
        
		iq = i - 3;
		iqq = i + 4;
	}
    
	if(iq < 1) iq = 1;
    if(iqq > m) iqq = m;
    
	for(i = ip - 1; i < ipp; i ++)
	{
		b[i -ip + 1] = 0.0;
        
		for(j = iq - 1; j < iqq; j ++)
		{
			h = z(i,j);
            
			for(k = iq - 1; k < iqq; k ++)
			{
				if(k != j)
					h = h * (v - y[k]) / (y[j] - y[k]);
			}
			
			b[i- ip + 1] = b[i - ip + 1] + h;
		}
	}
    
	w = 0.0;
    
	for(i = ip - 1; i < ipp; i ++)
	{
		h = b[i - ip + 1];
        
		for(j = ip - 1; j < ipp; j ++)
		{
			if(j != i)
				h = h * (u - x[j]) / (x[i] - x[j]);
		}
        
		w =w + h;
	}
    
	return(w);
}

//#include "Interpolation.inl"		//類及相關函數的定義頭文件

#endif		//_INTERPOLATION_INL

?? 快捷鍵說明

復制代碼 Ctrl + C
搜索代碼 Ctrl + F
全屏模式 F11
切換主題 Ctrl + Shift + D
顯示快捷鍵 ?
增大字號 Ctrl + =
減小字號 Ctrl + -
亚洲欧美第一页_禁久久精品乱码_粉嫩av一区二区三区免费野_久草精品视频
精品精品国产高清一毛片一天堂| 精品国产制服丝袜高跟| 国产精品久久久久久久蜜臀| 99视频一区二区| 日韩黄色一级片| 日韩精品亚洲一区二区三区免费| 成人丝袜高跟foot| 91久久一区二区| 69久久夜色精品国产69蝌蚪网| 欧美在线视频日韩| 欧美优质美女网站| 日韩亚洲欧美一区| 免费高清成人在线| 国产午夜精品久久久久久久| 自拍偷拍亚洲欧美日韩| 中文字幕制服丝袜成人av| 水野朝阳av一区二区三区| 成人激情视频网站| 日韩美女一区二区三区四区| 中文av一区二区| 国产精品伦一区二区三级视频| 亚洲国产精品精华液ab| 日韩av电影免费观看高清完整版| 久久成人18免费观看| 国产精品538一区二区在线| av在线不卡网| 91片黄在线观看| 久久综合九色综合欧美亚洲| 国产精品网站一区| 欧美国产日产图区| 中文字幕巨乱亚洲| 国产精品国产馆在线真实露脸| 精品视频全国免费看| 亚洲摸摸操操av| 亚洲欧美区自拍先锋| 艳妇臀荡乳欲伦亚洲一区| 欧美日韩情趣电影| 国产盗摄精品一区二区三区在线 | 夜夜亚洲天天久久| 欧美成人免费网站| 在线观看免费一区| 国产**成人网毛片九色| 午夜视频在线观看一区二区三区| 国产亚洲欧美激情| 777色狠狠一区二区三区| 成人免费高清视频在线观看| 日韩av一级片| 亚洲一级电影视频| 日本一区二区高清| 欧美成人精品高清在线播放| eeuss鲁片一区二区三区在线观看| 欧美一区二区三区思思人| 成人毛片在线观看| 日本aⅴ精品一区二区三区| 亚洲色图欧洲色图| 欧美日韩免费电影| 麻豆精品一区二区| 国产精品 欧美精品| 久久综合九色综合97_久久久| 欧美在线观看视频在线| 成人污视频在线观看| 狠狠色综合播放一区二区| 日日夜夜精品视频免费 | 精品国精品国产| 欧美日韩国产天堂| 在线观看亚洲精品| 91黄色免费看| 色94色欧美sute亚洲线路二| 99国产精品国产精品毛片| 国产成人h网站| 国产一区二区三区美女| 久久不见久久见免费视频1| 日韩激情视频网站| 日韩精品91亚洲二区在线观看| 亚洲在线视频网站| 一区二区三区日韩精品| 亚洲欧美国产高清| 亚洲欧美另类久久久精品| 中文字幕色av一区二区三区| 中文字幕av不卡| 国产精品久线在线观看| 亚洲欧洲精品一区二区三区| 国产精品久久久久久一区二区三区| 国产精品人人做人人爽人人添| 国产日韩欧美高清| 国产精品久久看| 亚洲欧洲精品成人久久奇米网| 国产精品传媒入口麻豆| 成人免费在线视频| 一区二区三区 在线观看视频| 一区二区三区中文字幕电影| 亚洲激情图片一区| 亚洲国产精品久久久久秋霞影院| 亚洲成av人在线观看| 首页国产欧美久久| 精品一区二区三区免费毛片爱| 国内精品久久久久影院一蜜桃| 国产精品18久久久久久久久久久久| 国产高清久久久| 成人黄色大片在线观看| 色婷婷国产精品| 欧美日韩国产一级| 精品欧美一区二区在线观看| 国产午夜亚洲精品午夜鲁丝片| 中文字幕的久久| 亚洲在线中文字幕| 懂色av一区二区在线播放| 欧美一区二区二区| 亚洲欧美日韩国产一区二区三区 | 国产+成+人+亚洲欧洲自线| 免费高清在线视频一区·| 一区二区视频在线| 午夜不卡av在线| 国产精品69毛片高清亚洲| 色偷偷一区二区三区| 欧美精品第1页| 日韩午夜电影在线观看| 国产成人精品亚洲777人妖| 97精品国产97久久久久久久久久久久| 91丝袜美腿高跟国产极品老师| 欧美日韩高清一区| 日本一区二区三区四区| 天天av天天翘天天综合网| 国产福利不卡视频| 欧美日韩久久久久久| 国产日韩欧美一区二区三区乱码 | 欧美一区二区三区视频在线| 麻豆freexxxx性91精品| 不卡的av中国片| 国产精品一二三区| 欧美在线免费视屏| 久久伊99综合婷婷久久伊| 亚洲精品国产精华液| 久88久久88久久久| 欧美视频自拍偷拍| 国产精品入口麻豆九色| 老司机精品视频导航| 91国模大尺度私拍在线视频| 欧美精品一区二| 亚洲成av人影院在线观看网| 91小视频在线观看| 2014亚洲片线观看视频免费| 蜜臀精品一区二区三区在线观看 | 99re热这里只有精品视频| 日韩精品一区在线| 亚洲成人免费视| 色哟哟国产精品免费观看| 国产三级一区二区| 激情综合色丁香一区二区| 精品视频资源站| 一区二区三区四区视频精品免费 | 欧美xxxxx裸体时装秀| 亚洲国产精品久久不卡毛片| 国产高清无密码一区二区三区| av高清久久久| 风间由美一区二区av101| 一本大道av一区二区在线播放| 国产精品国产三级国产专播品爱网| 精久久久久久久久久久| 欧美巨大另类极品videosbest | 欧美日本一区二区| 欧美tk—视频vk| 综合久久久久久| 一区二区三区中文字幕电影| 国产麻豆视频一区| 欧美一级艳片视频免费观看| 一区二区三区在线视频观看58| 日本亚洲一区二区| 久久电影国产免费久久电影| 欧美日韩精品一区二区| 亚洲色欲色欲www在线观看| 国产成人啪午夜精品网站男同| 精品国产伦一区二区三区观看体验 | 国产精品成人免费精品自在线观看| 日本麻豆一区二区三区视频| 欧美日韩综合不卡| 国产精品国产精品国产专区不片| 成人97人人超碰人人99| 黄色成人免费在线| 亚洲成人精品一区| 亚洲欧洲美洲综合色网| 中文字幕一区二区三区蜜月| 99久久综合99久久综合网站| 国产视频一区二区三区在线观看| 精品一区二区久久| 精品国产sm最大网站| 国产成人在线视频网站| 精品99999| 久久av资源站| 久久婷婷综合激情| 成人sese在线| 成人欧美一区二区三区白人| 99久久综合国产精品| 337p日本欧洲亚洲大胆色噜噜| 国产91精品免费| 国产精品美女www爽爽爽| 99久久亚洲一区二区三区青草| 国产精品色在线观看| k8久久久一区二区三区| 亚洲色图制服诱惑|