there is a polynomial: $$p(x)=1\cdot x^3+bx^2+cx+d$$
And there is a matrix of form - Toeplitz matrix with coeffcients of $p(x)$ on main diagonal: $$P=\pmatrix{1&b&c&d&0&0&0\cr0&1&b&c&d&0&0\cr0&0&1&b&c&d&0\cr0&0&0&1&b&c&d}$$
The output matrix is : $$P'=\pmatrix{1&0&0&0&x_1&x_2&x_3\cr0&1&0&0&x_4&x_5&x_6\cr0&0&1&0&x_7&x_8&x_9\cr0&0&0&1&x_{10}&x_{11}&x_{12}}$$
Is there some way to compute $P'=\frac{P}{p(x)}$. It is easy to compute $P'$ from $P$ by adding rows, but maybe there is a way to show that $P'=\frac{P}{p(x)}$. I have no idea how to mix matrices with polynomials.
