#include<stdio.h>
#include<math.h>
#define x0 1.5
#define MAXREPT 1000
#define EPS 0.5E-5
#define G(x) x-(x*x*x-x-1)/(3*x*x-1)
void main()
{
double x_k=x0,x_k1=x0;
int i;
printf(" k xk\n");
for(i=0;i<MAXREPT;i++)
{
printf(" %d %g\n",i,x_k1);
x_k1=G(x_k);
if(fabs(x_k1-x_k)<EPS)
{
printf("The Root is x=%g,k=%d\n",x_k1,i);
return;
}
x_k=x_k1;
}
printf("After %d repeate,no solve.\n",MAXREPT);
}
#include<stdio.h>
#include<math.h>
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余1页未读,立即下载