Вот расчеты для произведения 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)