Расширьте произведение сумм в Sympy, введя новые фиктивные индексы ⇐ Python

[Anonymous]

Я хочу использовать Sympy для упрощения некоторых выражений, содержащих произведения сумм. Для этого потребуется разложить произведения и отменить равные члены.

Код: Выделить всё

i, j, k, l = symbols('i, j, k, l', cls=Dummy)
A = IndexedBase('A')
N = symbols('N', integer=True, positive=True)
s = Sum(A[i]*A[j], (i, 1, N), (j, 1, N))
expr = s**2

Выполнение expr.doit(), expr.expand() или expr.simplify() оставит выражение без изменений.
Выражение отображается как


Как расширить продукт и ввести фиктивные индексы, чтобы преобразовать выражение в следующий вид?



Подробнее здесь: https://stackoverflow.com/questions/782 ... my-indices