Подсчёт корней на промежутке изоляции

#include<iostream>
#include<cmath>
#include<iomanip>
 
using namespace std;
 
double f(double x)
{
    return sqrt(x) - pow(cos(x), 2) - 2;
}
 
int main()
{
    double a, b, n, h, x, y, result;
    int k;
    setlocale(LC_ALL, "rus");
    a = 4; b = 8; n = 5;
    h = (b - a)/n;
    cout <<"a = " << a << endl;
    cout <<"b = " << b << endl;
    cout <<"n = " << n << endl;
    cout <<"h = " << h << endl;
    k = 0;
    cout <<setw(15)<<"x"<<setw(15)<<"f(x)"<<endl;
    for(x = a; x <= b; x = x + h)
    {
        cout <<setw(15)<<x<<setw(15)<<f(x);
        if(x < b && f(x)*f(x + h) < 0)
        {
            cout <<" Промежуток изоляции "<<x<<"; "<<x + h<<endl;
            k = k + 1;
        }
        else
        {
            cout <<endl;
        }
    
    }
    cout <<"Количество корней k = "<<k<<endl;
 
    do  
    {
        y = result;
        result = a - (f(a) * (b - a))/(f(b) - f(a));
        a = result;
        cout << "Roots:\n" << result << endl;
    }
    while (fabs(y - result) >= h);
    return 0;
}

for(x = a; x <= b; x = x + h)
{
    cout <<setw(15)<<x<<setw(15)<<f(x);
    if(x < b && f(x)*f(x + h) < 0)
    {
        cout <<" Промежуток изоляции "<<x<<"; "<<x + h<<endl;
        k = k + 1;
 
        double a1 = x;
        double b1 = x + h;
        double result = f(a1);
        do  
        {
            y = result;
            result = a1 - (f(a1) * (b1 - a1))/(f(b1) - f(a1));
            a1 = result;
        }
        while (fabs(y - result) >= h);
        cout << "Root #" << k << ": " << result << endl;
    }
    else
    {
        cout <<endl;
    }
   
}