模拟开平方根

linha0

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. void main()
   
5. {
   
6. #if(1)
   
7.     double Sqrtl( double dn);
   
8.     double n;
   
9.    
   
10.     scanf("%lf",&n);
    
11. 
    
12.     printf("平方根为：%.4lf\n",Sqrtl(n));
    
13.     printf("实际平方根：%lf\n",sqrt(n));
    
14. #endif
    
15. }
    
16. //开平方根
    
17. double Sqrtl( double dn)
    
18. {
    
19.     double n=0,sn=0;
    
20.     n=dn/2;
    
21. 
    
22.     while((sn=n*n)>dn)//获取粗略接近真实平方根的值n
    
23.     {
    
24.         n/=2;
    
25.     }
    
26.     while((dn-sn)>0)//获取最接近真实平方根的值n
    
27.     {
    
28.         if((dn-sn)>n*2) //减少循环次数
    
29.             n++;
    
30.         else
    
31.             n=n+1e-4;//控制精度
    
32.         sn=n*n;
    
33.     }
    
34.     return n;
    
35. }

复制代码

输入一个数：
49
输出结果：
平方根为：7.0001

有时，结果是正确的（7.0000），但大多数时间都是（7.0001）。并且，测试其它结果时是正确的。比如：9，16，10000

不知道这是为什么呢？还有，我想知道系统的Sqrt（）是怎么做的，不知道有什么办法没有？

仰望天上的光

n=n+1e-4;//控制精度
写为n=n+1e-5;或者更小就可以得到7了。
要想sqrt函数效率更好，LZ可以百度 牛顿迭代  平方根

linha0

即使用：n=n+1e-5;
也会出现7.00001的；之前我是用n=n+1e-3;的，不过他出现7.001……:@

linha0

而且，通过单步调试，按照逻辑来说，也是不应该出现那点小尾巴的！！:funk:

linha0

1. //牛顿迭代法求平方根,效率特别高
   
2. double SqrtN(double n)
   
3. {
   
4.     double xn=n/2,yn=0;
   
5.     int i=0;//计算循环次数
   
6. 
   
7.     while(yn!=xn)//7的平方根=2.6457513110646，编译器的最大精度；如果能自主控制那就更好了（如：yn=2.645和xn=3.432比较之类）
   
8.     {
   
9.         yn=xn;
   
10.         xn=(xn+n/xn)/2;
    
11.         i++;
    
12.     }
    
13. 
    
14.     return xn;
    
15. }

复制代码