牛顿迭代法源程序

#include
#include


int main()
{ int i=0; /*定义迭代次数*/
FILE *fp;
float a,x,X; /*定义a开平方，迭代前的值为x,迭代后的值为X。*/
fp=fopen("牛顿法迭代求a的平方根数据分析.doc","w");
scanf("%f%f",&a,&X); /*输入a的值以及x的迭代初值*/
fprintf(fp,"牛顿迭代格式为：X=0.5*(x+a/x)\n以a=%f,初始值x=%f为例\n",a,X);
do
{ x=X;
i++;
X=0.5*(x+a/x); /*迭代格式*/
printf("迭代次数:i=%d\nx=%f,X=%f\n",i,x,X);

fprintf(fp,"迭代次数:i=%d时\n迭代前的值x=%f,迭代后的值X=%f\n",i,x,X);

}while(fabs(x-X)>0.0000001);
printf("%f的平方根是%f\n",a,X);
fprintf(fp,"%f的平方根是%f\n",a,X);
fclose(fp);
return 0;
}