Нахождение всех корней на интервале

#include <iostream>
#include <cmath>
 
using namespace std;
const double pi = 3.14;
double Eps = 1e-15;
 
double f(double x)
{
    return(2 * pi * cos(x) - x + 1);
}
 
int main()
{
 
    double a (-6.28),
           b (-3.14),
           c,
           tmp;
 
 
    int n {0};
 
    while ((b - a) > Eps)
    {
        n++;
        tmp = c;
        c = (a + b) / 2;
        if (f(a) * f(c) < 0)
        {
            b = c;
        }
        else
            a = c;
    }
 
}

#include <iostream>
#include <cmath>
 
using namespace std;
 
double sign(double x)
{
    if (x>0) 
       return 1;
    else if (x<0)
       return -1;
    else
       return 0;
}
 
double root(double f (double), double a, double b, double eps=1.0e-8)
{
    double fc,fa=f(a), fb=f(b),c;
    
    while(abs(b-a)>=eps)
    {
        c=0.5*(a+b);
        fc=f(c);
        
        if (sign(fc)<=eps) return c;
        
        if (sign(fa)*sign(fc) < 0)
        {
            b=c;
            fb=fc;
        }
        else
        {
            a=c;
            fa=fc;
        }
    }
    return 0.5*(a+b); 
}
 
void allRoots(double f (double), double a, double b, double h=1.0e-2, double eps=1.0e-8)
{
    double x=a,f1=f(a),f2,r;
    
    while (x<=b)
    {
        x=x+h;
        f2=f(x);
        if (sign(f1)*sign(f2)<0)
        {
            r=root(f,x,x-h,eps);
            cout << r << endl;
        }
        f1=f2;
    }
    
}
 
double g(double x)
{
    return(2 * 3.1415926535 * cos(x) - x + 1);
}
 
int main()
{
 
    allRoots(g,-6.28,7.5);
 
    return 0;
}