用牛顿迭代法求浮点数的平方根

比如我们要求a的平方根，首先随便猜一个近似值x，然后不断令x等于x和a/x的平均数，迭代几次后x的值就已经相当精确了。

看下面(假设a=2,我们求2的平方根)：

1.先随便猜个数，比如我猜2的平方根为7

( 7  + 2/7 ) / 2 = 3.64287514
( 3.64287514  + 2/3.64287514 ) / 2 = 2.095946017

......

......

( n+ 2/n) / 2 = 1.414....

下面是代码：

1. #define ABS(VAL) (((VAL)>0)?(VAL):(-(VAL)))   
2. //用牛顿迭代法求浮点数的平方根   
3. float mysqrt(float x) {   
4.     float g0,g1;   
5.     if(x==0)   
6.         return 0;   
7.     g0=x/2;   
8.     g1=(g0+x/g0)/2;   
9.     while(ABS(g1-g0)>0.01)   
10.     {   
11.         g0=g1;   
12.         g1=(g0+(x/g0))/2;   
13.     }   
14.     return g1;   
15. }