Вот расчеты для произведения A, B и C на s.
Итак, А. с * Б. с - С. s:
t = [-88.0, 592.666, -1063.777, 805.833, -294.777, 51.5, -3.444]
Результат получается расширением в виде полиномов и оперированием с ними, а не оперированием матрицами, как это делается здесь.
Z = [24, -50, 35, -10, 1]
Этот результат здесь.
h = t / Z = [-3.666, 17.055, -3.444]
Мне удалось получить этот результат в виде серии Тейлора на виджете Wolfram Alpha.
Для получения результата вы можете использовать эти значения:
(-(3+4/9)*x⁶+51.5*x⁵-(294+7/9)*x⁴+(805+5/6)*x³-(1063+7/9)*x²+(592+2/3)*x-88)/((x-1)(x-2)(x-3)(x-4))
Or:
(-(3+4/9)*x⁶+51.5*x⁵-(294+7/9)*x⁴+(805+5/6)*x³-(1063+7/9)*x²+(592+2/3)*x-88)/(x⁴-10x³+35x²-50x+24)
Or:
P(x) = (-(3+4/9)*x⁶+51.5*x⁵-(294+7/9)*x⁴+(805+5/6)*x³-(1063+7/9)*x²+(592+2/3)*x-88)
Q(x) = (x⁴-10x³+35x²-50x+24)