牛頓迭代法求解平方根

#include<stdio.h>
#include<math.h>
int main()
{
double a,x;
double y=0;
a = 17895.4;//測試數據
x = (int)(a/2);//x0
while(fabs(x-y)>=1.0e-5){
y = x;
x = 0.5*(x+a/x);
}
printf("sqrt(%lf)=%lf",a,x);
return 0;
}

 

詳細理論基礎能夠參考數學分析的教材。。