bisekcja

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

double funkcja(double a)
{
    return ((((a-2)*a+3)*a-4)*a+1)*a+0.1;
}

double mz(double xl, double xr, double eps)
{
    double xm, ym, yl, yr;
    if (xl>xr) {
    xm=xl; xl=xr; xr=xm;
    }
    yl=funkcja(xl);
    yr=funkcja(xr);

    while(xr-xl>eps)
    {
        xm=(xl+xr)/2.0;
        ym=funkcja(xm);
        if(ym==0) return xm;
        if (yl*ym<0) {
        xr=xm; yr=ym;
        }
        else {
        xl=xm; yl=ym;
        }
    }
    return xm;
}



int main()
{
    int ilosc;
    scanf("%d", &ilosc);
    double lewyp, prawyp;
    double dokladnosc;
    for (int i=0; i<ilosc; i++)
    {
        printf("Podaj lewy przedzial: ");
        scanf("%lf", &lewyp);
        printf("Podaj prawy przedzial: ");
        scanf("%lf", &prawyp);
        printf("Podaj dokladnosc: ");
        scanf("%lf", &dokladnosc);
        printf("Miejsce zerowe: %.6lf\n", mz(lewyp, prawyp, dokladnosc));
    }

}