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

float epsi();
float funcao(float x);

int main(){

    float x1, x2, pontoMedio;

    printf("x1 e x2: ");
    scanf("%f%f", &x1, &x2);

    if(funcao(x1) * funcao(x2) < 0){
        pontoMedio = (x1+x2)/2;
        while(funcao(pontoMedio) > epsi()){
            if(funcao(x1) * funcao(pontoMedio) < 0)
                x2 = pontoMedio;
            else
                x1 = pontoMedio;
            pontoMedio = (x1+x2)/2;
        }
        printf("A raiz eh %.5f", pontoMedio);
    }
    else
        printf("Nao ha raizes no intervalo dado :c !");


    return 0;

}

float epsi(){

    double aux = 1, soma;

    soma = 1 + aux;

    while(soma > 1){

        aux /= 2;

        soma = aux + 1;

    }

    return aux * 2;

}

/*
float funcao(float x){

    return pow(x ,2) - pow(2, x);

}
*/

float funcao(float x){

    return pow(x, 3) - 9 * x + 3;

}