/* roots.hh */
#ifndef ROOTS_HH
#define ROOTS_HH

/* Function for computing the roots of a quadratic. */
void compute_roots( double A, double B, double C, double &root1, double &root2);

#endif /* end of ROOTS_HH */