solve(x.F.x+a*x.F.y,a) /*=> a=x.F.x/x.F.y*/ can write as the following using itensor solve(X([],[i])*F([i,j],[])*X([j],[]) + a* X([],[i])*F([i,j],[])*Y([j],[]),a ) /*=> a=X([j], [])/ Y([j], [])*/ But, result is not same. How to solve equations which use itensor?