splin2.cpp

Visual C++ 常用数值算法集 源代码

void SPLIN2(double X1A[], double X2A[], double YA[][11], double Y2A[][11], int M, int N, double X1, double X2, double& Y)
{ 
	double YTMP[100], Y2TMP[100], YYTMP[100];
	int J,K;
    for (J = 1; J<=M;J++)
	{
        for (K = 1; K<=N;K++)
		{
            YTMP[K] = YA[J][K];
            Y2TMP[K] = Y2A[J][K];
		}
        
        SPLINT(X2A, YTMP, Y2TMP, N, X2, YYTMP[J]);
    }
    SPLINE(X1A, YYTMP, M, 1E+30, 1E+30, Y2TMP);
    SPLINT(X1A, YYTMP, Y2TMP, M, X1, Y);
}