D = polytopes.dodecahedron()
DH = D.Hrepresentation()
X = list(SR.var('x_%i' % i) for i in (0..2))
Ineq=[[el.A()*vector(X) >=el.b()] for el in DH]
show(Ineq[0][0])
sol=[solve(Ineq[i][0],x[0]) for i in range(len(Ineq))]
This is nearly a code send by Slelievre as an answer to a former question. In his case it works nicely. But I need to use indexed variables. And in my case Sagemath complains that 'sage.symbolic.expression.Expression' object is not subscriptable'. Why ?