Binary_interval_method

#include<stdio.h>
#define N 100
float cal(float p)
{
float sum=p*p*p+p-4;
return sum; 
}
float main(){
float a;
float b;
float flag;
int n;
float p;
float FP;
float FA;
flag=0.001f;
n=0;
a=1;
b=3; 
if((FA=cal(a))==0)
{
return a;
}
if((FP=cal(b))==0)
{
return b;
}
while(n<N||n==N)
{
p=(a+b)/2.0;
FP=cal(p);
n++;
if(FP==0||(b-a)/2.0<flag)
{
printf("x: %f  N:%d ok \n",p,n);

return p;
}
else 
if(FA*FP>0)
{
a=p;
FA=FP;
}
else 
b=p;
}

printf("error\n");
return 0;  
}