Zoltan, my code has still problems with the original function from your
question (returns imaginary unit in the answer on my computer). Add two
rectform commands on the line starting with fac: ..., like this
fac:
fac*(ratsimp(expand((n-rectform(rhs(sol[i])))*(n-conjugate(rectform(rhs(sol[i])))))))^multiplicities[i])
Best regards
Robert
Kov?cs Zolt?n wrote:
> Dear Robert, thank you for your help. WMI2 (development version) is
> already using your code. ;-)
>
> What I changed is
>
> if expand(expr)#expand(fac) then error(),
>
> to
>
> if (not (expand(expr) = expand(fac))) then error(),
>
> because WMI2 uses # for internal operations. (Probably this should be
> changed someday to another symbol, not used by Maxima. Any suggestions?)
>