Решение нелинейного уравнения методом простой итерации

# include <iostream.h>
# include <math.h>
# include <windows.h>
# include <iomanip.h>
 
int counter=0;
double Root (double,double,double);
double fun (double);
 
int main ()
{
    SetConsoleCP (1251);
    SetConsoleOutputCP (1251);
 
    double x,a,b,eps,tmp;
    cout <<"Введите a="; cin >>a; 
    cout <<"Введите b="; cin >>b; 
 
    if (a>b)
    {
        tmp=a; a=b; b=tmp;
    }
 
    cout <<"Введите eps="; cin >>eps; cout <<endl;
 
    x=Root (a,b,eps);
    cout <<"Корень уравнения x="<<x<<endl;
    cout <<"Количество итераций="<<counter<<endl;
 
    return 0;
}
 
double Root (double a,double b,double eps) 
{
    if (fun(a)*fun(b)>0)
    {
        cout <<"На отрезке ["<<a<<";"<<b<<"] нет корня."<<endl;
        exit (1);
    }
 
    double x0, x1=a;
    do
    {
        counter++;
        x0=x1;
        x1=fun(x0);
    }
    while (fabs(x1-x0)>eps);
 
    return (x1+x0)/2;
}
 
double fun (double x)
{
    return 2*atan(x)-0.5*pow(x,3);
}