#include <iostream>
#include <cstdlib>
#include <math.h>

using namespace std;

double f(double x);
double biseccion ( double a, double b, double tol, int n);
int main()
{
    double a, b, tol, raiz;
    int n;
    cout<< "Ingrese el n\xA3mero de iteraciones:  ";
    cin>>n;
    cout << "Ingrese a:  ";
    cin>>a;
    cout<< "Ingrese b:  ";
    cin>>b;
    cout<< "Ingrese la tolerancia:  ";
    cin>>tol;
    raiz=biseccion(a,b,tol,n);
    cout<<"La ra\xA1z es: "<< raiz <<endl;
    system("pause");
    return 0;
}
 
 double f(double x)
 {
        return x*x*x+4*x*x-10;
 }
 double biseccion(double a, double b, double tol, int n)
 {
        double c;
        int nolter=0;
        do
        {
            c=(a+b)/2;
            if(f(a)*f(c)<0)
            {
               b=c;
            }
            else
            {
               a=c;
            }
            cout<<nolter<<"\t"<<a<<"\t"<<b<<"\t"<<c<<"\t"<<f(c)<<endl;
            nolter++;
         }
         while((fabs(f(c))>tol)&&(nolter<n));
         return c;
 }