Упрощение нахождения передаточной функции при помощи Maxima

[Hitmanmix]

Уважаемы форумчане, столкнулся с такой проблемой. Иногда по структурной схеме нужно определить передаточную функцию, используя стандартные преобразования. Пусть, для примера, имеется такая незаурядная структурная схема.



Пусть для примера будут вот такие передаточные функции в звеньях.

$W1(p)=\frac{2 \cdot p}{5\cdot p^2+7\cdot p+21}$

$W2(p)=\frac{5\cdot p}{9\cdot p+2}$

$W3(p)=\frac{8 \cdot p^2+1}{2\cdot p^2+3\cdot p}$

Очевидно, что итоговая передаточная функция такой структурной схемы будет иметь вид.

$W(p)=\frac{W1(p)\cdot W2(p)}{1-W1(p)\cdot W2(p)\cdot W3(p)}$

Так вот вопрос. Можно ли про помощи пакета Maxima найти передаточную функцию. В чем проблема, как известно, при выполнении преобразований p сокращать нельзя, из-за этого понижается порядок дифференциального уравнения, стандартная упрощалка Maxima та же самая ratsimp прекрасно сократит p и выражение будет неверным. То же эффект даст и collect из MathCad. Вопрос, может как-то можно Максиме запретить сокращать p, может есть какая-то директива или типа того, подскажите, если не трудно

[Hitmanmix]

То есть я хочу сказать, что если сделать вот так



То получается неверный результат из-за того, что Максима сокращает p. А вот какая передаточная функция должна быть у такой схемы, это я в ручную выполнил преобразования, не сокращая p



Вот и вопрос, может есть какая функция в Максима или код какой, чтобы сразу получить результат, как на последней картинке!