Sub SPLIN2(X1A(), X2A(), YA(), Y2A(), M, N, X1, X2, Y)
    Dim YTMP(100), Y2TMP(100), YYTMP(100)
    For J = 1 To M
        For K = 1 To N
            YTMP(K) = YA(J, K)
            Y2TMP(K) = Y2A(J, K)
        Next K
        Call SPLINT(X2A, YTMP, Y2TMP, N, X2, YYTMP(J))
    Next J
    Call SPLINE(X1A, YYTMP, M, 1E+30, 1E+30, Y2TMP)
    Call SPLINT(X1A, YYTMP, Y2TMP, M, X1, Y)
End Sub