平方根-泰勒展开式求法

平方根求法

Question Description:

已知: $x^2=N$, 求$x$。
问题转化为$f(x)=N-x^2=0$,解$x$。

---

f(x)的泰勒展开式

• $f(x)$在$x=x_0$处展开：
  
  $$f(x)=\sum_{n=0}^{\infty} \frac{f^{(n)}(x_0)}{n!} \times (x-x_0)^n$$
  $$f(x)=f(x_0)+f\prime(x_0)(x-x_0)+\frac{1}{2}f\prime\prime(x_0)(x-x_0)^2+O(x^2)$$
  • 可得：
    $$f(x)=N-x^2=N-x_0^2-2x_0(x-x_0)+O(x)$$
  • 令$f(x)=0$,得：
    $$x=\frac{1}{2}(x_0+N/x_0)$$

---

code

// c++ code
double ans=1, pre=0;
while(abs(ans-pre)>1e-6){
    pre=ans;
    ans=(ans+x/ans)/2;
}