lab7, C (clang) - rextester

lab7

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

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 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", newton(f1,       /*Цільова функція*/
                         0.41,     /*Перше обмеження, a*/
                         1,      /*Друге обмеження, b*/
                         110.0,      /*Параметер, t*/
                         0.0001     /*Точність*/
                        )); 
    return 0;
}


run \| edit \| history \| help	0