牛顿迭代法

#include<stdio.h>
#include<math.h>

double f(double x)
{   
 double a;
 a=x;
    return(a*a*a-3*a-1);
}
 
double f1(double x)
{  
 double b;
 b=x;
 return(b-f(b)/(3*b*b-3));
}

void main()
{
 double x0,x1;
 double e=1e-6;
 int k=1,N=20;
 int flag;
 printf("input x0 and N:/n");
 scanf("%lf %d",&x0,&N);
 while(k!=N)
 {
       if(f(x0)==0)
    {   
     flag=0;
           break;
    }
       x1=f1(x0);
       if(fabs(x1-x0)<e)
    {
    flag=1;
    break;
    }
     k++;
     x0=x1;
 }
     if(flag==1)
   printf("X1:%.3f/n",x1);
  else if(flag==0)
   printf("X0:%.3f/n",x0);
  else
   printf("ERROR!/n");
}