At the end I wanted to substitute qij^2=qij, i,j=1,2,3,4 the error: SyntaxError: keyword can't be an expression
the code: f1=(2*x1)+x2+1;
f2=-x1+(5*x2)-2;
e1=expand(f1^2);
e2=expand(f2^2);
var('q11,q12,q13,q14,q21,q22,q23,q24')
(q11,q12,q13,q14,q21,q22,q23,q24)
E1=e1.substitute(x1=q11+(2q12)-q13-(2q14), x2=q21+(2q22)-q23-(2q24));
E2=e2.substitute(x1=q11+(2q12)-q13-(2q14), x2=q21+(2q22)-q23-(2q24));
E1a=expand(E1);
E2a=expand(E2);
E1aa=E1a.substitute(q11^2=q11,q12^2=q12,q13^2=q13,q14^2=q14,q21^2=q21,q22^2=q22,q23^2=q23,q24^2=q24);
E2aa=E2a.substitute(q11^2=q11,q12^2=q12,q13^2=q13,q14^2=q14,q21^2=q21,q22^2=q22,q23^2=q23,q24^2=q24);
E1aa,E2aa