Здравствуйте.
Пишу программу для интерполяции сеточной функции сплайнами.
Использую следующие уравнение:

Получаю уравнения для расчета коэффициента c:



Получаю трех диагональную матрицу:

где

Собственно, решить такую матрицу методом прогонки не могу. Потому что по формулам количество получаемых коэффициентов A, B, C, D не равно, в силу того что некоторые из них выходят за пределы массива точек. Где D матрица правых частей уравнения.
Тогда я нашёл такое решение:
https://imageup.ru/img82/3558408/slide_13.jpgМатрица дополняется нулями и единицами. И это решение работает
Собственно вопрос в том откуда нули и единицы берутся? Я уже понял, что вроде как из краевых условий, но как именно не понял.