在多项式中代入矩阵时出现的问题

Question

示例：让

M = Matrix([[1,2],[3,4]]) # and 
p = Poly(x**3 + x + 1)    # then
p.subs(x,M).expand()

给出错误：

类型错误：无法添加 和

这是非常合理的，因为前两项成为矩阵，但最后一项（常数项）不是矩阵而是标量。为了补救这种情况，我将多项式更改为

p = Poly(x**3 + x + x**0)    # then

同样的错误仍然存​​在。我必须手动输入表达式，用 M 替换 x 吗？在这个例子中，多项式只有三个项，但实际上我遇到（多元多项式有）几十个项。