简单迭代

/*f(x)=(3x+1)^(1/3)*/  #include<stdio.h>   #include<math.h>  #define e 1e-5 double f(double x0)  {      double x=x0;       return pow((3*x+1),1.0/3.0);  }  int main()  {      double x0=1.5,x1;      x1=f(x0);      while(fabs(x1-x0)>e)      {          x0=x1;          x1=f(x0);      }      printf("%lf\n",x1);      return 0;  }