حل معادله نیوتون رافورس

#include
#include
#include
void main(void)
 {

  float h,x1,f,f1,g,a,b,x;
  int i,conter;

 top: clrscr();
      printf(" \n******************************in file**************************");
      printf("\ninput [a,b]");
      printf("\ninput a=");
      scanf("\n%f",&a);
      printf("input b=");
      scanf("\n%f",&b);
      printf("\nbaseh f(x)=>    [%f,%f]",a,b);
      x=a;
/****************   f(x)   ****************/
      f=x*pow(2,x)-1;
      g=x*x*2;
/*******************************************/
      printf("\nf(a)=%f",f);
      x=b;
/*****************   f(x)   ***************/
      f1=x*pow(2,x)-1;
/*******************************************/
      printf("\nf(b)=%f",f1);
      if ((f1*f)<0)
        printf("\nIF ONE TRUE");
      else
       {
   printf("\n IF ONE FALSE --F(a)F(B)<0--");
   goto top;
       }
//****************************************************************
      printf("\ninput x1:=");
      scanf("%f",&x1);
      printf("\nnumber repeter :=");
      scanf("%d",&conter);
      printf("\n************************out file****************************");

      for(i=1;i<=conter;i++)
    {
        h=x1-(f/g);
        printf("\n%f ",h);
        x1=h;
    }

      getch();
}