00001 00012 #include <iostream> 00013 #include <iomanip> 00014 #include <math.h> 00015 #include <stdlib.h> 00016 00017 using namespace std; 00018 00020 #define maxiter 5000 00022 #define DBL_EPSILON 1e-15 00023 00025 int roots(double *a,int n,double *wr,double *wi); 00026 00030 void get_quads(double *a,int n,double *quad,double *x);