Actually, I can use the Leibnitz rule for higher-order derivatives of

the product term for expansion, and then reverse the order of two sums and get

.

Then write the matrix form of this equation. However, I was looking for a different way.

For instance, having the coloumns of the unknown both involving and , if possible.