简单迭代

/*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;  
}