Метод простых итераций с заданной точностью на С++

#include <cmath>
#include <iomanip>
#include <iostream>
using namespace std;
 
double f(double x)
{
    return pow(2.0,x) - 2*pow(x, 2.0) - 1;
}
 
double g(double x)
{
    return x + 0.5*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;
        cout<<"x    = "<<x   <<endl;
        cout<<"g(x) = "<<g(x)<<endl;
        cout<<"f(x) = "<<f(x)<<endl;
        x = g(x);
    }
    system("pause");
    return 0;
}

#include <iomanip>
#include <iostream>
using namespace std;
 
 
inline double fx(double x, double y)
{
    return (-((x*x - 1)*(exp(y) + y * exp(x))) / (2 * x*exp(y) + 2 * y * exp(x)));
}
 
inline double gy(double x, double y)
{
    return (-((y*y - 1)*(x*exp(y) + exp(x))) / (2 * x*exp(y) + 2 * y * exp(x)));
}
 
int main()
{
    double x, y, xn, yn;
    double eps;
    cout << "Enter initial root value x= ";
    cin >> x;
    cout << "Enter initial root value y= ";
    cin >> y;
    cout << endl;
    cout << "Enter error of calculation : ";
    cin >> eps;
    for (double iter = 1; (eps < fabs(fx(x, y))) && (eps < fabs(gy(x, y))); iter = iter + 1)
    {
        system("cls");
        //*Итераций может быть очень много, поэтому рекомендую забыть
        //о целых а использовать дабл как счётчик, хотя в принципе если 
        //решение не нашли за 10-100 итераций то решения для данного коэффициента
        //при f(x,y) в g(x,y) нет и надо его менять
        cout << "Iteration : " << setprecision(0) << iter << endl;
        cout << "x = " << x << endl;
        cout << "y = " << y << endl;
        cout << "f(x,y) = " << fx(x, y) << endl;
        cout << "g(x,y) = " << gy(x, y) << endl;
        xn = x;
        yn = y;
        x -= fx(xn, yn);
        y -= gy(xn, yn);
 
        //system("pause");
 
    }
    system("pause");
    return 0;
}