Найти минимум и максимум заданной функции - C (СИ)

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define MAX_I 100000 //максимальное число итераций: так как не у всякой функции есть минимум или максимум, нужно учесть это, что бы программа была отказоустойчивой
 
double f (double x)
{
    return 0.5*x*x; //f возвращает значение туда, где она вызывается
}
 
int lin_search (double (*f)(double), double a, double b, double eps, double *x)
{
    double h = (b-a)/100; //это шаг линейного поиска, с шагом h мы будем идти по оси х и проверять значения функции
    double t = a;
    double prev = f(t); //это первое значение функции
    double cur = f(t+h); //это следующее
    int i;
    for (i=0; i < MAX_I; i++) //здесь ставим ограничение на тот самый MAX_I, что бы программа не зацикливалась
        {
            while (cur > prev && t+h < b) //просто проверяем значения функции, если cur > prev, то максимум еще не найден и нужно идти дальше
                {
                    t+=h;
                    prev = cur; //здесь продвигаем cur и prev на один шаг h
                    cur = f(t+h);
                }
            //если есть максимум на этом отрезке, то найдется cur < prev и цикл оборвется
            //в этом случае поворачиваем назад и уменьшаем шаг поиска для увеличения точности
            h /= -10;
            if (fabs(h) < eps) break; //когда h станет достаточно малым, то выходим из цикла
        }
    if (i < MAX_I) // проверяем не зациклилась ли у нас функция, если нет, то меняем х и возвращаем положительный результат
        {
            *x = prev;
            return i;
        }
    return -1; //иначе возвращаем отрицательное значение
}
 
int main (void)
{
    double a, b, eps, x;
    printf("Input a: "); scanf("%lf", &a);
    printf("Input b: "); scanf("%lf", &b);
    printf("Input eps: "); scanf("%lf", &eps);
    
    if ( lin_search(f, a, b, eps, &x) > 0 ) printf("x = %lf  -  f(x) = %lf\n", x, f(x));
        else printf("Bad luck, dude\n");
    return 0;
}