Нелинейное уравнение методом Ньютона

#include <stdio.h>
#include <math.h>
 
//Объявляем функцию f(x)
double f(double x)
{
    double f;
    f = sqrt(x) + log(x) + sin(4 * x) - 2;
    return f;
}
 
//находим первую производную
double f1(double x)
{
    double f1;
    f1 = 4 * cos(4 * x) + 1/x + 1/(2 * sqrt(x));
    return f1;
}
 
//находим вторую производную
double f2(double x)
{
    double f2;
    f2 = -(16 * sin(4 * x) + pow(x, -2) + 1/(4 * pow(x, 3/2)));
    return f2;
}
 
//
void main() {
    system("chcp 1251");
    //объявляем переменные;
    //a и b - заданный интервал поиска корня;
    //x - искомы корень;
    //e - точность вычисления;
    //k - количество итераций поиска;
    double a, b, e, x;
    int k;
    a = 2.2;
    b = 2.8;
    e = 0.00001;
    printf("a = %f\n", a);
    printf("b = %f\n", b);
    printf("e = %f\n", e);
    if(f(a) * f2(a) > 0) {
        printf("Условие на сходимость выполнено для x = a = %f\n", a);
        x = a;
    }
    else {
        if(f(b) * f2(b) > 0) {
            printf("Условие на сходимость выполнено для x = b = %f\n", a);
            x = b;
        }
        else {
            printf("Условие на сходимость не выполнено\n");
            x = -10E10;
        }
    }
    if(x > -10E10) {
        k = 0;
        while(1) {
            x = x - f(x)/f1(x);
            k+=1;
            if(fabs(f(x)) < e)
            break;
        }
        printf("Корень x = %f\n", x);
        printf("Количество итерация k = %d\n", k);
    }
}