Delphi Анализ корреляционной матрицы

[ZHU35]

есть 2 матрицы XY- матрица исходных данных и корреляционная матрица К(квадратная матрица).
кол-во столбцов в матрицах одинаковое: XY.N=K.N
структура К:
x1x1 x1x2 x1x3 x1y
x2x1 x2x2 x2x3 x2y
x3x1 x3x2 x3x3 x3y
yx1 yx2 yx3 yy

структура XY:
x1 x2 x3 y

нужно
в первую очередь: проверить последнюю стоку и если abs(yxn)<0.1 то удалить из XY столбец n

следующее: проверить взаимодействия xnxk, если abs(xnxk)>0,9, то удалить из XY тот столбец для которого взаимодействие с y - меньше (если abs(yxn)<abs(yxk) , то удаляем n)

функции удаления столбцов есть.

проблема в том, что если на первом этапе удалить средний столбец (например x2), то возникает вопрос как не рассматривать взаимодействия x2xn.

XY.N - естественно любое положительное число >2, у - не трогаем;

пример К:
+1.00 - 1.00 +0.93 +0.75
- 1.00 +1.00 +0.60 +0.05
+0.93 +0.60 +1.00 +0.95
+0.75 +0.05 +0.95 +1.00