Найти корни уравнения методом Ньютона

#include "stdafx.h"
#include <cmath>
#include <iomanip>
#include <iostream>
using namespace std;
 
double f(double x)
{
     return pow(0.2*(x),3)-cos(x) ;//моя функия
}
 
 
 
double g(double x)
{
    return x - f(x);
}
 
int main()
{
    double x;
    double eps;
    cout<<"Enter initial root value   : ";cin>>x;
    cout<<"Enter error of calculation : ";cin>>eps;
    for(double iter = 1; eps < fabs(f(x)); iter = iter + 1)
    {
        system("cls");
        //*Итераций может быть очень много, поэтому рекомендую забыть
        //о целых а использовать дабл как счётчик, хотя в принципе если 
        //решение не нашли за 10-100 итераций то решения для данного коэффициента
        //при f(x) в g(x) нет и надо его менять
        cout<<"Iteration : "<<setprecision(0)<<iter<<endl;
        if(f(x) == 0)//Чёртовский важный момент(!)
            break;//ведь если df(x) == 0, то будет деление на ноль x - f(x)/df(x)
        cout<<"x    = "<<x    <<endl;
    
     
        cout<<"f(x) = "<<f(x) <<endl;
        x = g(x);
        
    return 0;
 
    }
}

#include <iostream>
#include <cmath>
 
using namespace std;
 
double f(double x)      //возвращает значение функции f(x) = (0.2x)^3-cos(x)
{
    return (0.008*x*x*x-cos(x));
}
double df(double x)     //возвращает значение производной
{
    return (0.024*x*x+sin(x));
}
double d2f(double x)    //возвращает значение второй производной
{
    return (0.048*x+cos(x));
}
 
int main()
{
    int exit = 0, i;    // переменные для выхода из цикла
    double x0, xn;      // вычисляемые приближения для корня
    double a, b, eps;   // границы отрезка и необходимая точность
    do 
    {
        i=0;
        cout<<"Please input [a b]\n=>";
 
        while(1)
        {
            while (!(cin >> a >> b) || (cin.peek() != '\n'))
            {
                cin.clear();
                cin.ignore(20, '\n'); 
                cout << "Input error! Retry input" << endl;
            }
            //Проверяем условия:
            //f(a)*f(b) < 0 (функция принимает значения разных знаков на концах отрезка [a;b]);
            //Производные f'(x) и f"(x) сохранят знак на отрезке [a;b],
            //т.е. функция либо возрастает на отрезке, либо убывает,
            //сохраняя при этом направление выпуклости
            if ((f(a)*f(b)<0) && (df(a)*df(b)>0) && (d2f(a)*d2f(b)>0) )
                break;
            cout << "Input error! Extremum! Retry input" << endl;
        }
 
        // вводим нужную точность вычислений
        cout<<"\nPlease input epsilon\n=>";
        while (!(cin >> eps) || (cin.peek() != '\n'))
        {
            cin.clear();
            cin.ignore(10, '\n'); 
            cout << "Input error! Retry input" << endl;
        }       
 
        //выбираем за начало тот конец отрезка, где касательная к кривой y=f(x) на отрезке [a;b] пересекает ось Ох
        x0 = (df(a)*a-f(a))/df(a);          //абцисса, в которой касательная пересекает ось Ох
        x0 = ((x0>=a) && (x0<=b)) ? a : b;  //Начальная точка из условия a < f'(x)*x-f(x))/df(x) < b
        
        xn = x0-f(x0)/df(x0);               // считаем первое приближение
            
        cout<<++i<<"-th iteration = "<<xn<<"\n";    //выведем
            
        while(fabs(x0-xn) > eps)        // пока не достигнем необходимой точности, будет продолжать вычислять
        {
            x0 = xn;
            xn = x0-f(x0)/df(x0);       // непосредственно формула Ньютона
            
            cout<<++i<<"-th iteration = "<<xn<<"\n";//выведем
        }
            
        cout<<"\nRoot = "<<xn;          // вывод вычисленного корня
 
        cout<<"\nExit?=>";
        cin>>exit;
    } while (exit!=1);  // пока пользователь не ввел exit = 1
    return 0;
}