lab7_OOP 0.3 alha, C (clang)

lab7_OOP 0.3 alha

#include  <stdlib.h>
#include  <stdio.h>
#include  <math.h>
#include  <float.h>

double solve (double(*method)(double(*)(double,double),double,double,double,double),double(*f)(double,double),double a,double b,double t,double eps)
{
    return method(f, a, b, t, eps);
}

double f1 (double x, double t)
{
    if (x != 0.0)
    {
        return cos(t/x) - 2.0*sin(1.0/x) + 1.0/x;
    }
    else
    {
        printf("Помилка. Ділення на нуль. f1(%lf, %lf)\n", x, t);
        exit(-1);
    }
}

double f2 (double x, double t)
{
    if (x > 0.0)
    {
        return sin(log(x)) - cos(log(x)) + t*log(x);    }
    else
    {
        printf("Помилка. Функція не визначена на проміжку. f2(%lf, %lf)\n", x, t);
        exit(-1);
    }
}

double bisec (double (*f)(double, double), double a, double b, double t, double eps )
{
    return 0.0;
}

double newton (double (*f)(double, double), double a, double b, double t, double eps )
{
    if (f != NULL)
    {
        double x = b;
        double delta = DBL_MAX;
        int iterations = 0;
        for (; fabs(delta) >= eps && iterations <= 1000; ++iterations)
        {
            double diff = (f(x + 1e-10, t)-f(x, t))/1e-10;
            if (diff != 0.0)
            {
                delta = f(x, t)/diff;
                x -= delta;
            }
            else
            {
                printf("Помилка. Похідна дорівнює нулю. newton(%p, %lf, %lf, %lf, %lf)\n", f, a, b, t, eps);
                exit(-1);
            }
        }
        if (iterations > 1000)
        {
            printf("Помилка. Кількість ітерацій перевищила 1000. newton(%p, %lf, %lf, %lf, %lf)\n", f, a, b, t, eps);
            exit(-1);
        }
        return x;
    }
    else
    {
        printf("Помилка. Нульовий покажчик на функцію. newton(%p, %lf, %lf, %lf, %lf)\n", f, a, b, t, eps);
        exit(-1);
    }
}

int main(void)
{
    printf("%lf", solve(newton,    /*Метод обчислень*/
                         f1,       /*Цільова функція*/
                         -2.0,      /*Перше обмеження, a*/
                         2.0,      /*Друге обмеження, b*/
                         2.0,      /*Параметер, t*/
                         0.0001     /*Точність*/
                        ));
    
    return 0;
}


run \| edit \| history \| help	0