MS>Насколько я помню, данный алгоритм принимает матрицы след размерностей
MS>

MS>A * B = C, где A = dim(N * N), B = dim(1 * N), C = dim(1 * N)
MS>

MS>Вопрос. Может ли данный алгоритм принять матрицу B (а это матрица неизвестных переменных) размерность N*N?

Метод трехдиагональной прогонки решает лишь ленточные (трехдиагональные) системы.
то бишь те, для которых

Уравнение A*B=C

для любого i, для любого j, если |i — j| > 1, то A_i_j = 0

Зато решает хорошо и резво