Где m - порядок дифференциальногоуравнения; К - число дискретных точек
отсчета; i - порядок производной;
- дискретные значения искомого решения, записываемые в матрицу "Z"; H - кубическая матрица известных коэффициентов.
Ниже приводится алгоритм решения для случая 3-х диагональной матрицы. Надо сказать,
что методы "стрельбы", "прогонки" являются частными случаями метода Гаусса. Параметры
дифференциальногоуравнения
следующие: m=2,далее см. по тексту.
Иллюстрацию этого метода проведем на решении уравнения:
,с граничными условиями
(здесь
Z
i,k
,
i=1..
m
-1
отсутствуют - есть только
Z
o
,k
. Поэтому у
Z
можно оставить один индекс, например,
k
).