I am new to SymPy. These are my definitions of symbols :
n = symbols(\'n\', integer=True, positive=True) rho = symbols(\'rho\', integer=True) c = IndexedBase(