Although solve() is not really designed to deal with itensor expressions, it
can be used, so long as you avoid reusing the same dummy index in your
equation. I.e., try
solve(X([],[i])*F([i,j],[])*X([],[j]) + a* X([],[u])*F([u,v],[])*Y([],[v]),a
);
(Note also that I corrected the position of the index in the second
occurrence of X in the numerator and in Y in the denominator; as you
originally wrote it, the same index occurred twice in a covariant position,
which is not correct.)
Viktor
-----Original Message-----
From: maxima-bounces at math.utexas.edu [mailto:maxima-bounces at math.utexas.edu]
On Behalf Of Niitsuma Hirotaka
Sent: Thursday, September 02, 2010 12:50 AM
To: maxima at math.utexas.edu
Subject: How to solve equation using itensor
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?
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima