Заданы коэффициенты полинома 3 степени и точность. Найти корень уравнения методом ньютона с введенной точностью?

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include <locale.h>
 
 
int main()
{
   setlocale(LC_ALL,"Rus");
   double y,a0,a1,a2,a3,x,f,f1;
   printf("Введите значения коэфициэнтов a0, a1, a2, a3:");
   scanf("%lf %lf %lf %lf", &a0, &a1, &a2, &a3);
   double Lb3 (double x)
   {
    return a0+a1*x+a2*pow(x,2)+a3*pow(x,3);
   }
   
}

#include <iostream>
#include <cmath>
 
double f(double x, double a0, double a1, double a2, double a3)
{
    return a0 + a1*x + a2*x*x + a3*x*x*x;
}
 
double df(double x, double a1, double a2, double a3)
{
    return a1 + 2.*a2*x + 3.*a3*x*x;
}
 
double findRoot(double x, double eps, double a0, double a1, double a2, double a3)
{
    double h = f(x, a0, a1, a2, a3) / df(x, a1, a2, a3);
    while (fabs(h) >= eps)
    {
        h = f(x, a0, a1, a2, a3) / df(x, a1, a2, a3);
        x = x - h;
    }
    return x;
}
 
int main()
{
    double a0, a1, a2, a3, x0, eps;
    std::cout << "a0 = "; std::cin >> a0;
    std::cout << "a1 = "; std::cin >> a1;
    std::cout << "a2 = "; std::cin >> a2;
    std::cout << "a3 = "; std::cin >> a3;
    std::cout << "x0 = "; std::cin >> x0;
    std::cout << "eps = "; std::cin >> eps;
    std::cout.precision(16);
    std::cout << findRoot(x0, eps, a0, a1, a2, a3);
}

#include <iostream>
#include <cmath>
double a0, a1, a2, a3;
double f(double x)
{
    return a0 + a1 * x + a2 * x * x + a3 * x * x * x;
}
double df(double x)
{
    return a1 + 2. * a2 * x + 3. * a3 * x * x;
}
double findRoot(double x, double eps)
{
    double h = f(x) / df(x);
    while (fabs(h) >= eps)
    {
        h = f(x) / df(x);
        x -= h;
    }
    return x;
}
int main()
{
    double x0, eps;
    std::cout << "a0 = "; std::cin >> a0;
    std::cout << "a1 = "; std::cin >> a1;
    std::cout << "a2 = "; std::cin >> a2;
    std::cout << "a3 = "; std::cin >> a3;
    std::cout << "x0 = "; std::cin >> x0;
    std::cout << "eps = "; std::cin >> eps;
    std::cout.precision(16);
    std::cout << findRoot(x0, eps);
    return 0;
}

while (fabs(h) >= eps)
{
    h = f(x) / df(x);
    x -= h;
}