c program for arriving at a solution using bisection method

c program for arriving at a solution using bisection method

main()
{
 float a,b,c,x1,x2,x,series;
 double d;

 printf("enter a,b,c and x1(pos) & x2(neg)");
 scanf("%f%f%f%f%f", &a, &b, &c, &x1, &x2);

 read:

 x = (x1 + x2) / 2;
 series = a * x * x + b * x + c;
 d = fabs(series);

 if (d > 0.0001)
 {
  if (x * x1 < 0)
   x = x2;
  else
   x = x1;

  goto read;
 }
 else
 {
  printf("ans=%f", x);
 }

 return 0;
}