迭代公式:
/*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;
}