Численные методы решения уравнений (Метод Ньютона)

#include <iostream>
using namespace std;
double find(double x, double eps)
{
    double func(double x) {  //функция, описывающая уравнение
        return (x - 1 / (3 + sin(3.6 * x)));
    }
    double diff(double x) {  //функция, описывающая производную уравнения
        double h = 0.01;
        return (func(x + h) - func(x - h)) / (2.0 * h);
    }
    double answer() {  //функция для вычисления корня уравнения методом Ньютона
        double a = 0.1, b = 0.85, e = 0.0001, x[100];
        int i = 0;
        x[i] = 0 - func(0) / diff(0);
        while (a <= b or x[i + 1] - x[i] <= e) {  //цикл для нахождения корня с заданной точностью
            x[i + 1] = a - func(a) / diff(a);
            a += 0.01;
        }
        return x[i];
    }
    int main()
    {
    cout << "Задание \n";
    cout << "Точное значение = " << 0.2624 << '\n';
    cout << "Корень, найденый методом Ньютона = " << answer() << '\n';  //вывод на экран результата функции, которая вычисляет корень уравнения методом Ньютона
    cout << "Погрешность вычислений = " << abs(0.2624 - abs(answer())) << '\n'; //сравнение точного значения и найденного 
    system("pause");
}

while (a <= b or x[i + 1] - x[i] <= e) {  //цикл для нахождения корня с заданной точностью
            x[i + 1] = a - func(a) / diff(a);
            a += 0.01;
        }

while (fabs(x[i + 1] - x[i]) >= e) {  //цикл для нахождения корня с заданной точностью
            x[i + 1] = x[i] - func(x[i]) / diff(x[i]);
        }

#include <iostream>
using namespace std;
 
int sum(int a, int b)
{
    return a + b;
}
 
int product(int a, int b)
{
    return a * b;
}
 
int foo(int a, int b, int (*function)(int, int))
{
    return function(a, b);
}
 
int main()
{
    cout << foo(2, 5, &sum) << endl;
    cout << foo(2, 5, &product);
}

#include <iostream>
using namespace std;
 
int sum(int a, int b)
{
    return a + b;
}
 
int product(int a, int b)
{
    return a * b;
}
 
typedef int (*Function)(int, int);
 
int foo(int a, int b, Function f)
{
    return f(a, b);
}
 
int main()
{
    cout << foo(2, 5, &sum) << endl;
    cout << foo(2, 5, &product);
}

x-= func(x) / diff(x);