首先是最基本的二分开根号,这个比较容易理解,复杂度比起下面讲的牛顿迭代法要高,更容易理解。
下面给出代码:
#define eps 0.00001
float SqrtByDichotomy(float n)
{
if(n<0)
{
return -1.0;
}
else
{
float low,up,mid,last;
low=0,up=(n>=1?n:1);
mid=(low+up)/2;
do
{
if(mid*mid>n)
up=mid;
else
low=mid;
last=mid;
mid=(up+low)/2;