Программа не доходит до конца

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <locale.h>
double f (double x);
double f1 (double y);
double f2 (double z);
 
 
int main()
{
    int i=0;
    double a, b, x, x1;
    setlocale (LC_ALL, "RUS");
    printf ("Введите отрезок [a; b]\n");
    printf ("Введите а: ");
    scanf_s ("%f", &a);
    printf ("\nВведите b: ");
    scanf_s ("%f", &b);
    if (f(b)*f2(b)>0)
        printf ("\nТочка с абсциссой b является неподвижным концом хорды.\nНачальное приближение равно a.");
    else
    {
        if (f(a)*f2(a)>0)
            printf ("\nТочка с абсциссой а является неподвижным концом хорды.\nНачальное приближение равно b.\n");
    }
    printf ("\nНачальное приближение Х0 = %.6f\n", b);
    x=b;
    do 
    {
        i++;
        x1=x;
        x=x1-((f(x1)*(a-x1))/(f(a)-f(x1)));
    }
    while (fabs(x-x1)>=0.000001);
    printf ("Получили приближенное значение кореня x=%.7f\n", x);
    printf ("Число итераций i=%d\n", i);
    printf ("Невязка = %.8f\n", fabs(f(x)));
    system ("pause");
    return 0;
}
 
double f (double x)
{
    double fun;
    fun = pow(x, 3)+6*x*x-4*x-24;
    return fun;
}
 
double f1 (double y)
{
    double fun1;
    fun1 = 3*y*y+12*y-4;
    return fun1;
}
 
double f2 (double z)
{
    double fun2;
    fun2 = 6*z+12;
    return fun2;
}