Derivatives of symbolic matrix power

[Upabjojr]

We need support for derivatives of matrix expressions with symbolic powers.

[nileshpatra]

can I take up this issue?

[Upabjojr]

can I take up this issue?

if you can find a way to solve this, feel free to open a PR.

[oscarbenjamin]

I guess that this is not the ideal result:

In [83]: n = symbols('n')

In [84]: A = MatrixSymbol('A', 3, 3)

In [85]: diff(A**n, A)
Out[85]: ArrayDiagonal(ArrayTensorProduct(I, I, n*HadamardPower(A, n - 1)), (1, 4), (3, 5))