非线性方程求解（Newton迭代法通用子程序代码）

#include<math.h>
double f(double x)
{
return(pow(x,6)-x-1);
}
/*******F(X)******/
double g(double x)
{
return(6*pow(x,5)-1);
}
/*******G(X)*****/
main()
{
double x0,x,p;
int i=2,k,N0;
printf("\n enter x0,N0,p:\n");
scanf("%lf,%d,%lf",&x0,&N0,&p);
for(k=0;k<=N0;k++)
  { if(g(x0)==0)
                {       i=1;
                        printf("%d\n",i);
                        break;
                }
                x=x0-f(x0)/g(x0);
printf("\n %d %12.8lf %8.3le",k,x0,f(x0));
if(fabs(x-x0)<p)
{
i=0;;
break;
}