Решение не линейных уравнений за численными методами - C (СИ)

# include <stdio.h>
# include <math.h>
# include <stdlib.h>
# include <time.h>
# include <sys/time.h>
# define eps 0.001
# define E 2.71
# define n 100000
 
Функція знаходження розв’язку рівняння:
 
double Func(double x){
    return (x-1)*(x-1)-0.5*pow(E,x);
}
 
Диференціал функції знаходження розв’язку рівняння:
 
double dFunc(double x){
    return 2*(x-1)-0.5*pow(E,x);
}
 
Ітератор функції знаходження розв’язку рівняння:
 
double iFunc(double x){
    return 0.5*(0.5*pow(E,x)-x*x);
}
 
Метод половинного поділу:
 
double PolovPodil(double a, double b){
    double x;
    int i;
    for (i=0;i<n;i++)
        do{
            x=(a+b)/2;
            if (Func(x)*Func(a)>0) a=x;
            else b=x;
        } while (((a-b)/2)>eps);
    return x;
}
 
Метод хорд:
 
double Hordy(double a, double b){
    double x=a,xn;
    int i;
    for (i=0;i<n;i++)
        do{
            xn=x;
            x=a-((b-a)/(Func(b)-Func(a)))*Func(a);
            if (Func(x)*Func(a)>0) a=x;
            else b=x;
        } while (abs(x-xn)>eps);
    return xn;
}
 
Метод Ньютона:
 
double Newton(double a, double b){
    double x,xn;
    x=(a+b)/2;
    int i;
    for (i=0;i<n;i++)
        do{
            xn=x;
            x=xn-Func(xn)/dFunc(xn);
            if (Func(a)*Func(x)>0) a=x;
            else b = x;
        } while (abs(x-xn)>eps);
    return x;
}
 
Метод Ньютона-Рафсона:
 
double NewtonRafson(double a){
    double x=a;
    int i;
    for (i=0;i<n;i++)
        do{
            x=x-Func(x)/dFunc(x);
        } while (Func(x)>eps);
    return x;
    }
 
 
Метод простої ітерації:
 
double Iteration(double a){
    double x=a,xn;
    int i;
    for (i=0;i<n;i++)
        do{
            xn=x;
            x=iFunc(xn);
        } while (abs(x-xn)>eps);
    return x;
}
 
Головна програма:
 
void main(){
    double a=-100.0,b=0.0,h=0.5;
    clock_t timer;
    while (a<=100){
        b=a+h;
        if (Func(a)*Func(b)<0){
            timer=clock();
            printf("PolovPodil: %g\n",PolovPodil(a,b));
            timer=clock()-timer;
            printf("Processing (100000 times) time is: %.4f seconds\n",(double)timer/CLK_TCK);
            timer=clock();
            printf("Hordy: %g\n",Hordy(a,b));
            timer=clock()-timer;
            printf("Processing (100000 times) time is: %.4f seconds\n",(double)timer/CLK_TCK);
            timer=clock();
            printf("Newton: %g\n",Newton(a,b));
            timer=clock()-timer;
            printf("Processing (100000 times) time is: %.4f seconds\n",(double)timer/CLK_TCK);
            timer=clock();
            printf("NewtonRafson: %g\n",NewtonRafson(a));
            timer=clock()-timer;
            printf("Processing (100000 times) time is: %.4f seconds\n",(double)timer/CLK_TCK);
            timer=clock();
            printf("Iteration: %g\n",Iteration(a));
            timer=clock()-timer;
            printf("Processing (100000 times) time is: %.4f seconds\n",(double)timer/CLK_TCK);
            
        }
        a=b;
    }
}