Массив указателей на функции - C (СИ) - Ответ 9791900

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
 
typedef double(*func_t)(double);
 
int i;
double derivative(double, double, func_t)  ;
double second_derivative(double, double, func_t f) ;
double newton(double, double, func_t);
double func1(double), func2(double), func3(double), func4(double), func5(double), func6(double), func7(double), func8(double);
int main(void)
{
    int maxIteration, i, funcNum;
    double x0,xn;// вычисляемые приближения для корня
    double a, b, c, h, eps, func;// границы отрезка и необходимая точность
    //Объявление массива ссылок на функции
    double (*mas[8])(double) = {func1, func2, func3, func4, func5, func6, func7, func8};
    printf( " Vvedite h " ); 
    scanf( "%lf", &h ); // вводим шаг
    printf( "  Vvedite eps " );
    scanf( "%lf", &eps ); // вводим нужную точность вычислений
    printf( "  Vvedite maximum number of iteration: " );
    scanf( "%lf", &maxIteration );// вводим кол-во итераций
    for(funcNum = 0; funcNum < 8; funcNum++) { 
        printf("funcNum: %d\n", funcNum);
        if (mas[funcNum](a)*second_derivative(x0, h, mas[funcNum]) > 0) 
            x0 = a; // для выбора начальной точки проверяем f(x0)*p2f(x0)>0 ?
        else 
            x0 = b;
        c = (a + b) / 2; 
        if (mas[funcNum](c)*mas[funcNum](a) < 0)
            b=c;
        if (mas[funcNum](c)*mas[funcNum](b) < 0)
            a = c;
            xn = newton(x0, h, mas[funcNum]);
        while (fabs(x0-xn) > eps) {
            x0 = xn;
            xn = newton(x0, h, mas[funcNum]);
        } 
        while( xn < 0 ) {// пока не достигнем необходимой точности, будет продолжать вычислять
            x0+=h;
            xn = newton(x0, h, mas[funcNum]);
            while (fabs(x0-xn) > eps) {
                x0 = xn;
                xn = newton(x0, h, mas[funcNum]);
            }
            printf("xn = %.9f\n",xn); 
            printf("f(x0) = %.9f\n",mas[funcNum](x0)); 
            printf("Number of iteration = %d\n\n",i);
        } 
        printf( "  Exit?=> " );
        scanf( "%lf", &exit );
    }
    getchar();
}
 
//Функции-уравнения, решаемые в задаче
double func1(double x) { return exp(-x) + pow(x, 2) + 2; }
double func2(double x) { return -pow((x-2), 2); }
double func3(double x) { return exp(x)-2*x-2; }
double func4(double x) { return sin(x) + x - 1; }
double func5(double x) { return exp(-pow(x,2));}
double func6(double x) { return 2*exp(-3*x) - x +1; }
double func7(double x) { return log(x) - x + 1; }
double func8(double x) { return log(x) - pow((x-1), 2); }
 
//Первая производная функции
double derivative(double x, double h, func_t f)
{
    return (f(x + h) - f(x - h)/(2*h));
}
//Вторая производная функции
 
double second_derivative(double x, double h, func_t f)
{
    return (derivative(x, h, f) - derivative(x, h, f)/(2*h));
}
double newton(double x, double h, func_t f)
{
    return (x - f(x)/derivative(x, h, f));
}