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

#include <iostream>
#include <iomanip>
#include <cmath>
#include <stdexcept>
 
using namespace std;
 
double f(double x) {
    return x * x * x - 4.0 * x + 3;
}
 
typedef double (*pFunc)(double);
 
class FuncPoint { 
    double x, y;
    pFunc pf;
public:
    FuncPoint(pFunc _pf = nullptr) : x(0.0), y(0.0), pf(_pf) {
        if (pf)
            y = pf(x);
    }
    FuncPoint(double _x, double _y, pFunc _pf = nullptr) : x(_x), y(_y), pf(_pf) {
        if (pf)
            y = pf(x);
    }
    double getX() const { return x; }
    void setX(double _x) {
        x = _x;
        if (pf)
            y = pf(x);
    }
    double getY() const { return y; }
    void setY(double _y) {
        if (pf)
            throw logic_error("attempt to set dependent parameter");
        else
            y = _y;
    }
    pFunc getF() const { return pf; }
    void setF(pFunc _pf) {
        pf = _pf;
        if (pf)
            y = pf(x);
    }
};
 
double fixedPointIteration(pFunc pf, double a, double b, double eps) {
    FuncPoint p1(f), p2(f);
    p1.setX(a);
    p2.setX(b);
    if (p2.getY() * p1.getY() > 0.0)
        throw invalid_argument("root absent");
    double k = (p2.getY() - p1.getY()) / (p2.getX() - p1.getX());
    double xi = p2.getX();
    while (fabs(xi - p1.getX()) > eps) {
        xi = p1.getX();
        p1.setX(p1.getX() - p1.getY() / k);
        //cout << "Xi = " << xi << "    Xp = " << p1.getX() << "    delta = " << fabs(xi - p1.getX()) << endl;
    }
    return p1.getX();
}
 
int main()
{
    double x1, x2, eps, x;
    cout << "Enter x1:  ";
    cin >> x1;
    cout << "Enter x2:  ";
    cin >> x2;
    cout << "Enter eps: ";
    cin >> eps;
 
    try {
        cout << setiosflags(ios_base::scientific);
        x = fixedPointIteration(f, x1, x2, eps);
        cout << "\nroot = " << x << " f(x) = " << f(x) << endl;
    }
    catch (exception & e) {
        cout << "*** Exception: " << e.what() << endl;
    }
 
    return 0;
}