-----maxima-bounces at math.utexas.edu wrote: -----
>Oh. Am feeling a bit stupid! Which code is Barton Willis's?
That would be the code in /share/linearalgebra/polynomialp.lisp.
>In which case, do you mean polyp which is in limit.lisp? That works on
>all but your last example. And I really should have found it before
>trying to write my own!
The function polyp might be all that you need (but I haven't looked
at the source)
(%i10) ?polyp(rat(q,x));
(%o10) false
(%i11) ?polyp(x);
(%o11) true
Convert CRE to general form before applying ?polyp
(%i12) ?polyp(?ratdisrep(rat(q,x)));
(%o12) true
BW