Normal() again

[Christian Bauer]

Hi!

On Fri, 19 Nov 1999, Do Hoang Son wrote:
> > normal(q10^2-q10^4*(q10^2-q21^2)^(-1)+q21^2*q10^2*(q10^2-q21^2)^(-1));
> power::eval(): division by zero

Fixed.

> Maybe the next bug is related to the above one:

No.

> normal(q10^2*(I*Pi^2*(q10^2-q21^2)^(-1/2)*q10^2*(-(1/2)*(q10^2-q21^2)^(-1/2)*q10^2+(1/2)*q21^2*(q10^2-q21^2)^(-1/2))+1/2*I*Pi^2*q10^2));
> q10^2*(1/2*I*Pi^2*q21^2*(q10^2-q21^2)^(-1)*q10^2-1/2*I*Pi^2*(q10^2-q21^2)^(-1)*q10^4+1/2*I*Pi^2*q10^2)
> > normal(");
> 0

This is caused by (x^(1/2))^2 being treated (during normal()) as foobar^2
and not being further evaluated until foobar gets re-inserted in the very
end. A workaround is to call expand() before normal().

Bye,
Christian

-- 
  / Coding on PowerPC and proud of it
\/ http://www.uni-mainz.de/~bauec002/