史蒂芬森加速迭代

迭代公式：

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