Метод деления отрезка пополам

//---------------------------------------------------------------------------
#include <stdio.h>  // printf, scanf
#include <conio.h>  // getch
#include <math.h> // Математические функции (sin cos atan exp pow и т.д.) и константа Пи - M_PI
#include <windows.h>
#include <stdlib.h>
 
 
//---------------------------------------------------------------------------
 
 
double my_f1(double x){  // функция вариант 1 (один параметр типа double и результат double)
    return (pow(4,1.0/3)-pow(sin(x/10),2))/pow(x,1.0/2)-x;
}
 
double my_f4(double x){  // функция вариант 4 (один параметр типа double и результат double)
    return 1.5-(pow(x,1.0/2)+pow(x,1.0/3))/exp(3.0/2)-x;
}
 
double check_root(double a, double b, double (*f)(double))
{
    double fa = f(a), fb = f(b);
    if (fa*fb > 0) { printf("Корня на отрезке AB нет\n"); return 0; }
    else return 1;
}
 
double root_half(double a, double b, double (*f)(double), double eps, double &c, int &iter )
{
    double fa = f(a);
    do{
        c=(a+b)/2.0;
        double fc = f(c);
        if (fc*fa > 0) { a = c; } else { b = c; }
        iter++;
    } while (fabs(b-a)>=eps && iter <=2000);
    return 0;
}
double simple_iter(double a, double b, double (*f)(double), double eps, double &c, int &iter )
{
    double x0,x1,delta;
    x0=(a+b)/2.0;  //стартовое значение - середина отрезка ab
 
    do {
        x1=x0+0.5*f(x0);
        delta=fabs(x1-x0);
        x0=x1;
        iter++;
    } while (delta>=eps && iter<=2000);
    c=x1;
    return 0;
}
 
 
void print_header(int NFunc)
{
int i;
 
for (i=1; i < 80; i++) printf("="); printf("\n");
printf("| Функция №%d  |            Метод деления        |         Метод простых       |\n",NFunc);
printf("| Погрешность |           отрезка пополам       |            итераций         |\n");
for (i=1; i < 80; i++) printf("="); printf("\n");
printf("|             |    Корень      |   Итераций     |    Корень     |   Итераций  |\n");
for (i=1; i < 80; i++) printf("="); printf("\n");
}
 
 
 
void main()
{
    SetConsoleOutputCP(1251);
    int i,n,z, iter1 = 0 , iter2=0, iter3 = 0,iter4= 0;
    double a, b, root1=0, root2=0; // вещ.число с плавающей точкой двойной точности
    double eps = 0.1; //стартовая погрешность
    //вводим исходые данные
 
    printf("N (от 1 до 10)=? ");
    scanf("%d", &n);  // ввод десятичного(%d) n - целое число
    if (n<1 || n>10) {
        printf("Недопустимое значение N\n");
        printf("\nPress any key to exit");  getch(); return;
    }
    printf("Введите значение A (0<A<2) ");
    scanf("%lf", &a);  // ввод A - вещественное число;
    if (a<=0 || a>=2) {
        printf("Недопустимое значение A\n");
        printf("\nPress any key to exit");  getch(); return;
    }
    printf("Введите значение B (0<A<B<2) ");
    scanf("%lf", &b);
    if (b<=0 || b>=2 || b<=a) {
        printf("Недопустимое значение B\n");
        printf("\nPress any key to exit");  getch(); return;
    }
 
 
    print_header(1);
    if (!check_root(a,b,my_f1)==0) {
      for (i=1; i <= n; i++) {
        root_half(a,b,my_f1,eps, root1, iter1);
        simple_iter(a,b,my_f1,eps, root2, iter2);
        z = ceil(fabs(log(eps)/log(10.0)));
        printf("| %-12.*lf|  %-12.*lf  |     %4d       |  %-12.*f |    %4d     |\n",z, eps,z, root1,iter1,z, root2,iter2);
        eps = eps/10.0;
      }
    }
    //сбрасываем погрешность и итерации перед вторым вариантом
    eps = 0.1;
    print_header(4);
    if (!check_root(a,b,my_f4)==0) {
      for (i=1; i <= n; i++) {
        root_half(a,b,my_f4,eps, root1, iter3);
        simple_iter(a,b,my_f4,eps, root2, iter4);
        z = ceil(fabs(log(eps)/log(10.0)));
        printf("| %-12.*lf|  %-12.*lf  |     %4d       |  %-12.*lf |    %4d     |\n",z,eps,z, root1, iter3 ,z,root2,iter4);
        eps = eps/10;
      }
    }
 
 
    printf("\nPress any key to exit");  getch(); // задержка консольного экрана до нажатия любой клавиши
    return; //выход из main()
}
//---------------------------------------------------------------------------

double my_f1(double x){  // функция вариант 1 (один параметр типа double и результат double)
    return (pow(4,1.0/3)-pow(sin(x/10),2))/pow(x,1.0/2)-x;
    //return (log(44.8)-sin(pow(x,1.0/2)))/2-x;