Разные ответы в реализации Метода хорд на С++ и Python

#include <iostream>
#include <cmath>
using namespace std;
double f(double x)
{
    return ((double)((x - sin(x)) / ((double)1 - cos(x)) - ((double)5 / (double)3)));
}
 
double solve(double x0, double x1)
{
    float sol = 0;
    float delta = 0.001;
    while (abs(x1 - x0) > delta)
    {
        sol = x0 - f(x0) * (x1 - x0) / (f(x1) - f(x0));
        cout << " sol " << sol << endl;
        if (f(x0) * f(sol) < 0)
        {
            x1 = sol;
            cout << " x1 " << x1 << endl;
        }
        else
            if (f(x1) * f(sol) < 0)
            {
                x0 = sol;
                cout << " x0 " << endl;
            }
            else
            {
                cout << " break " << endl;
                return sol;
            }
        return sol;
    }
}
int main()
{
    double x0;
    double x1;
    double decision;
    cin >> x0 >> x1;
    decision = solve(x0, x1);
    if (decision == 0)
        cout << "No solution was found" << endl;
    cout << decision << endl;
 
}

/* Secant method (order of convergence = golden ratio = (1 + sqrt(5)) / 2 ≈ 1.618) */
 
#include <iostream>
#include <cmath>
 
double f(double x)
{
    return (x - sin(x)) / (1. - cos(x)) - 5./3.;
}
 
double findRoot(double a, double b, double eps)
{
    while(fabs(b - a) > eps)
    {
        a = b - (b - a) * f(b) / (f(b) - f(a));
        b = a - (a - b) * f(a) / (f(a) - f(b));
    }
    return b;
}
 
int main()
{
    std::cout.precision(16);
    std::cout << findRoot(3., 4., 1e-5);
}