Hi there,
A simple problem appeared and I can't solve it by myself. Plese don't focus on the meaning of the functions and of the values assigned to variables: I suspect the problem is on how I treat aa and bb.
After defining the following functions XX and Hvv that I can call in the code:
XX = -Inverse[aa].(bb.UU );
Hvv = cc.(Inverse[s i - aa].bb);
If I call
X=XX&[aa=4,bb=5,cc=6]
and then call Hvv
Hv=Hvv&[aa=6,bb=7]
I get wrong values for Hv.
But if after calling X=XX&(ecc.) I re-evaluate the definition of Hvv = cc.(Inverse[s i - aa].bb) then everything is fine if I call now Hv=Hvv(ecc.).
I guess I'm misunderstanding the meaning of f(aa)&(aa=a) which I intend as "evaluate f ifor aa=a".
You guys have insight on how to solve this? Thanks.