比如我们要求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. }