一个开n次方的函数sqrt()

Code:

1. #define PCS 1e-8   
2. //x为底数，y为指数   
3. double pow(double x,double y)   
4. {   
5.     double pow=1;   
6.     if(y==0) return 1.0;   
7.     else if(y>0){   
8.             while(y!=0){   
9.                 pow*=x;   
10.                 --y;   
11.             }   
12.             return pow;   
13.         }   
14.     else{   
15.         y=-y;   
16.         while(y!=0){   
17.             pow*=x;   
18.             --y;   
19.         }   
20.         pow=1.0/pow;   
21.         return pow;   
22.     }   
23. }   
24.   
25. //x为被开方数，n为开方次数   
26. double sqrt(double x,unsigned n){   
27.        
28.     double lower=0,higher=x;   
29.     double value=(lower+higher)/2;   
30.     if(n==1) goto end;   
31.     //利用二分法求解   
32.     while(pow(value,n)-x>PCS||pow(value,n)-x<-PCS){   
33.         if(pow(value,n)-x>0) higher=value ;   
34.         else lower=value;   
35.         value=(higher+lower)/2;   
36.     }   
37.     return value;   
38. end:   
39.     return x;   
40. }