sqrt()平方根计算函数的实现1——二分法

C语言标准库：

http://www.cplusplus.com/reference/cmath/

在一个区间中，每次拿中间数的平方来试验，如果大了，就再试左区间的中间数；如果小了，就再拿右区间的中间数来试。比如求sqrt(16)的结果，你先试（0+16）/2=8，8*8=64，64比16大，然后就向左移，试（0+8）/2=4，4*4=16刚好，得到了正确的结果sqrt(16)=4。

实现时我第一次没有进行精度控制，导致一些不能完整开方的数运行不出来，也忘记了对非整数的判断。代码1为改正后的代码，后面代码2为第一次写的代码，贴在这里以作对比提醒。

代码1：

esp常在C++中用来控制迭代精度，是函数程序事先声明的常量，在微积分中是无限小值。

float esp=0.000000;

double sqrt(double x)

{

double down,up,n,last;

down=0;

up=x;

n=(down+up)/2;

if(x<=0)

return x;

while((n*n!=x)&&abs(last-n)>esp)

{

if(n*n<x)

{

down=n;

last=n;

n=(down+up)/2;

}

else

{

up=n;

last=n;

n=(down+up)/2;

}

}

return n;

}

代码2：

double sqrt(double x)

{

double down,up,n;

down=0;

up=x;

n=(down+up)/2;

while(n*n!=x)

{

if(n*n<x)

{

down=n;

n=(down+up)/2;

}

else

{

up=n;

n=(down+up)/2;

}

}

return n;

}

参考：http://www.2cto.com/kf/201206/137256.html