We know that there are bugs that show up when gcd : subres
that disappear when gcd : spmod. I'm not sure if the following
is a bug or if it shows that some calculations can be
(very) slow with gcd : subres and fast with gcd : spmod. I've
never let the gcd : subres calculation run to completion.
gcd : spmod;
/* gcd : subres; */
showtime : all$
phi : rat(n0/(x-x0) + n1/(x-x1) - (n0 + n1)/(x-x2));
psi : rat(m0/(x-x0) + m1/(x-x1) + m2/(x-x2));
eq2 : phi * (diff(psi,x) + psi^2) - psi * (diff(phi,x) + phi^2)$
okay;
Should we change the default for gcd from subres to spmod?
Barton