#include
#define Cero 0.00001
float f(float x){
return ((x*x*x)-2);
}
float Regual_falsi (float a,float b, float tol){
float c = a - ( (b-a)/ (f(b)-f(a)) )*f(a);
int i=0;
while ( (f(c) != Cero) && ( fabs(b-a) > tol )&& i<10 ){
if (f(a) * f(c) < Cero)
b=c;
else
a=c;
c = a - ( (b-a)/ (f(b)-f(a)) )*f(a);
i++;}
return c;
}
int main(){
printf("el valor de la raiz es: %f \n", Regual_falsi (-2.0, 2.0, 0.000001 ) );
return 0;
}
No hay comentarios:
Publicar un comentario