/* cálculo de una raíz por el método de la bisección */
/* fichero bisec.c */
#include <stdio.h>
#include <math.h> 
#define EPS1 10E-8
#define EPS2 10E-12
void main(void)
{
        double a=1, b=4, c, fa, fb, fc;
        double func(double); 
        c = (a+b)/2.0;
        while (func(c)>EPS1 && (b-a)>EPS2){
                c = (a+b)/2;
                if (func(c)*func(a) > 0)
                        a = c;
                else 
                        b = c;
        }
        printf("El valor de la raíz es: %12.6e\n", c);
}
double func(double x)
{
        return (sin(x)-x);
}