加速迭代.cpp

数值计算方法的程序： 加速迭代 牛顿迭代 Gauss 杜里特尔分解法 克洛特 追赶法 平方根法 改进平方根法 雅可比迭代 高斯 赛德尔迭代 拉格朗日插值法 分段线性插值 分段

#include<stdio.h>
#include<math.h>
float y(float x0)
{return pow(1+x0,1.0/3);
}

float yd(float x0)
{return -1.0/3*pow(1+x0,-2.0/3);
}
float x(float x0,float c)
{return 1.0/(1-c)*(y(x0)-c*x0);
}

void main()
{float x0=1.5,x1,c;
printf("求x(3)-x-1=0在[1,2]的根,误差10(-3).\n");
c=yd(x0);
x1=x(x0,c);
for(;1;)
{if(fabs(x1-x0)>=1e-3)
{x0=x1;
x1=x(x0,c);}
else break;

}

printf("%f",x1);


}