I cannot see any differences using your latest code (with rectform).
1/(x^4+1) doesn't work, and, in general 1/(x^n+1) doesn't seem to work iff
n>=4.
B. r., Z.
2008/4/4, Kov?cs Zolt?n <kovzol at matek.hu>:
>
> For me, your previous code works correctly on Maxima 5.13. And it also
> works for 1/(x^2+1), 1/(x^2+3), 1/(x^4-1), but not for 1/(x^4+1). I'm trying
> your latest suggestion. Z.
>
> 2008/4/4, Robert Marik <marik at mendelu.cz>:
> >
> > 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?)
> > >
> > >
> >
>
>
> --
> matek.hu --- mert matekb?l ?retts?gizni mindenkinek kell...
>
--
matek.hu --- mert matekb?l ?retts?gizni mindenkinek kell...