Sindbad~EG File Manager

 int gsturm(POLYnomial, POLYnomial, term *);
 /* generalized Sturm sequence */
 int nroots(POLYnomial P);  /* number of real roots of P, or -1 for error */
 int coste_roy(POLYnomial P, POLYnomial Q, unsigned *cpos, unsigned *cneg);
/* number of real roots of P making Q positive (*cpos) or negative (*cneg) */
 int sylvester(POLYnomial P, POLYnomial Q, unsigned *cpos, unsigned *cneg);
/* same as coste_roy but assuming gcd(P,Q) = 1, so don't have to compute it */
 int nroots_interval(POLYnomial P, term a, term b);
/* number of roots of P in open interval (a,b) */