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
原文链接: https://www.cnblogs.com/williamjie/p/11131314.html
欢迎关注
微信关注下方公众号,第一时间获取干货硬货;公众号内回复【pdf】免费获取数百本计算机经典书籍
原创文章受到原创版权保护。转载请注明出处:https://www.ccppcoding.com/archives/298262
非原创文章文中已经注明原地址,如有侵权,联系删除
关注公众号【高性能架构探索】,第一时间获取最新文章
转载文章受原作者版权保护。转载请注明原作者出处!