Newton

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


/*#define F0 ((3.0*x0+1)/(x0*x0)) 
#define f0 (-(3.0*x0+2)/(x0*x0*x0))
*/
/*#define F0 (x0*x0*x0-1)/3.0
#define f0 x0*x0*1.0
*/
/*#define F0 x0*x0*x0-3*x0-1
#define f0 3*x0*x0-3
*/
  /*
#define F0 (x0*x0*x0-x0-1)
#define f0 (3*x0*x0-1)
*/
double calF(double x){
 return 4+2/3.0*cos(x);
}
double calf(double x){
return -2/3.0*sin(x);
}
int main()
{
    int n=1,cal=0;
    int N;
    long double x0,flag,x1;
    printf("请分别输入x0(初始值),N(控制迭代次数),C(控制迭代次数)\n");
scanf("%lf %d %lf",&x0,&N,&flag);
    while(n<N)
    {
               printf("迭代%d次数:  近似解x%d: %lf   F%d: %lf F'%d: %lf\n",cal,cal,x0,cal,calF(x0),cal,calf(x0));
               if(calf==0)
               {
                   printf("公式不符合条件\n");
                   return  0;
  }
               x1=x0-calF(x0)/calf(x0)*1.0;
  if(fabs(x0-x1)<flag)
               {
                   printf("方程的近似解为:%lf\n",x1);
                   ++cal;
  ++n;
                   return 0;
               }
               ++cal;
  ++n;
               x0=x1;
    }
    printf("方程的近似解为:%8f\n",x1);
    return 0;
}