Метод бисекций

typedef double (*doublefdouble) (double);
double abs (double x)
{
 if (x<0.0)
 {
  return -x;
 }
 return x;
}
double Bisec (double a, double b, double Epsilon, doublefdouble Func)
{
 double c;
 c=(a+b)/2;
 if (abs(b-a)<Epsilon)
 {
  return c;
 }
 if (Func(c)*Func(a)<=0.0)
 { 
  return Bisec(a,c,Epsilon,Func);
 }
 return Bisec(c,b,Epsilon,Func);
}

Bisec(a,b,Epsilon,f); // Вместо f подставить имя функции, вычисляющей левую часть уравнения. Значение Bisec - искомый корень.

double hhklk(double x)
{
 return x*x-1;
}

double x;
x=Bisec(0.0,20.0,0.01, hhklk); 
// Здесь x уже равен корню уравнения

#include "stdafx.h"
#include "conio.h"
#include "locale.h"
#include "iostream"
#include "math.h"
#include "function.h"
using namespace std;
 
int _tmain(int argc,_TCHAR* argv[])
{
    setlocale (LC_CTYPE,"Russian_Rusia.1251");
        double a=0;
        double b=2.0;
        double e=0.1;
        double x;
        for(;;)
        {
            double x1=(a+b)/2-e/2;
            double x2=(a+b)/2+e/2;
        if(f(x2)>f(x1))
            b=x2;
        else
            a=x1;
        double l=b-a;
        if(l<0.2)
            break;
    }
    x=(a+b)/2; 
    double min=f(x); 
    cout<<"\t"<<"Значение x при котором функция имеет минимум "<<x<<endl;
    cout<<"\t"<<"Значение функции в этой точке Y(x)="<<min<<endl;
    _getch();
    return 0;
}

#include "stdafx.h"
#include "cmath"
 
double f (double x)
{   
    double y;
    y=x*x*x*x-4/3*x*x*x+20/3*x*x-112/27*x+10;
 
    return(y);
}

double f (double x);