Sindbad~EG File Manager

MEXPORT_PROVER int gsturm(POLYnomial, POLYnomial, term *);
 /* generalized Sturm sequence */
MEXPORT_PROVER int nroots(POLYnomial P);  /* number of real roots of P, or -1 for error */
MEXPORT_PROVER int coste_roy(POLYnomial P, POLYnomial Q, unsigned *cpos, unsigned *cneg);
/* number of real roots of P making Q positive (*cpos) or negative (*cneg) */
MEXPORT_PROVER 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 */
MEXPORT_PROVER int nroots_interval(POLYnomial P, term a, term b);
/* number of roots of P in open interval (a,b) */