平方根原来可以这么算

mikkoyue

1. #include <stdio.h>
   
2. 
   
3. double sqrt(double x)
   
4. {
   
5.   [color=DarkOrange]double[/color] smaller, greater, result;
   
6. 
   
7.   smaller = 1.0;
   
8. 
   
9.   [color=DarkOrange]if[/color] (x > smaller) {
   
10.     greater = x;
    
11.   } [color=DarkOrange]else[/color] {
    
12.     greater = smaller;
    
13.     smaller = x;
    
14.   }
    
15. 
    
16.   result = greater;
    
17. 
    
18.   [color=DarkOrange]if[/color] ( result * result > x ) {
    
19.     result = smaller;
    
20.   }
    
21. 
    
22.   [color=DarkOrange]while[/color] (greater - smaller > .000000001) {     // 此处可以设置为你的精度要求
    
23.    [color=DarkOrange] if[/color] ( result * result < x ) {
    
24.       smaller = result;
    
25.     } [color=DarkOrange]else[/color] {
    
26.       greater = result;
    
27.     }
    
28.     result = (smaller + greater)/2;
    
29.   }
    
30. 
    
31.   [color=DarkOrange]return[/color] result;
    
32. }
    
33. 
    
34. int main()
    
35. {
    
36.   [color=DarkOrange]double[/color] a;
    
37. 
    
38.   printf("Please input a number to a:");
    
39. 
    
40.   scanf("%lf", &a);
    
41. 
    
42.   printf("sqrt(%g) = %g\n", a, sqrt(a));
    
43. 
    
44.   [color=DarkOrange]return[/color] 0;
    
45. }

复制代码

doubleelec

Post by mikkoyue

1. #include <stdio.h>
   
2. 
   
3. double sqrt(double x)
   
4. {
   
5.   [color=DarkOrange]double[/color] smaller, greater, result;
   
6. 
   
7.   smaller = 1.0;
   
8. 
   
9.   [color=DarkOrange]if[/color] (x > smaller) {
   
10.     greater = x;
    
11.   } [color=DarkOrange]else[/color] {
    
12.     greater = smaller;
    
13.     smaller = x;
    
14.   }
    
15. 
    
16. [color=Red] result = greater;
    
17. 
    
18.   if( result * result > x ) {
    
19.     result = smaller;
    
20.   }[/color]
    
21.   [color=DarkOrange]while[/color] (greater - smaller > .000000001) {     // 此处可以设置为你的精度要求
    
22.    [color=DarkOrange] if[/color] ( result * result < x ) {
    
23.       smaller = result;
    
24.     } [color=DarkOrange]else[/color] {
    
25.       greater = result;
    
26.     }
    
27.     result = (smaller + greater)/2;
    
28.   }
    
29. 
    
30.   [color=DarkOrange]return[/color] result;
    
31. }
    
32. 
    
33. int main()
    
34. {
    
35.   [color=DarkOrange]double[/color] a;
    
36. 
    
37.   printf("Please input a number to a:");
    
38. 
    
39.   scanf("%lf", &a);
    
40. 
    
41.   printf("sqrt(%g) = %g\n", a, sqrt(a));
    
42. 
    
43.   [color=DarkOrange]return[/color] 0;
    
44. }

复制代码

红线的地方不就是  result = smaller;   吗？ 还是我看错了？

另外，如果参数是负数，这个程序也能返回个实数值来。

mikkoyue

谢谢楼上的提醒，当参数是非正数时我忘记考虑了，我只考虑了参数大于0的情况，看来得加一句 if (x < 0) exit(-1); 才行，关于2楼提到的红色的部分，大家就想想吧，目前我也不是十分确定，反正这样是一定不会出错，呵。现在改正如下
[PHP]#include <stdio.h>

double sqrt(double x)
{
  double smaller, greater, result;

  if (x <= 0) exit(-1);

  smaller = 1.0;

  if (x > smaller) {
    greater = x;
  } else {
    greater = smaller;
    smaller = x;
  }

  result = greater;

  if( result * result > x ) {
    result = smaller;
  }
  while (greater - smaller > .000000001) {     // 此处可以设置为你的精度要求
    if ( result * result < x ) {
      smaller = result;
    } else {
      greater = result;
    }
    result = (smaller + greater)/2;
  }

  return result;
}

int main()
{
  double a;

  printf("lease input a number to a:");

  scanf("%lf", &a);

  printf("sqrt(%g) = %g\n", a, sqrt(a));

  return 0;
}[/PHP]

在gcc 3.4.1下编译通过　

lucifer

不如写个不用浮点数实现的，这样更有实用价值一些。

mikkoyue

嗯，有道理，努力中。。。

strangeryf

seems reinventing wheels

zouyuchong

笔算的话，我更倾向于连分数逼近。。。

虽然有笔算开方的方法，但连分数逼近方便一点。。

davidliang

不是吧
用乘的啊

我学数值分析的时候好象用什么 牛顿中值定律(?)来逼近的
收敛很快
r*r算法好象差了些....

可惜学得差了些,忘了

ljil

win98+lcc编译成功，呵呵