求求帮忙

lifhv

用泰勒公式求sinx和cosx的近似值，精确到10e-5

jackz007

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. /* 本程序利用泰勒级数近似计算 sin(x) 和 cos(x) 的数值，精确到 10e-5 */
   
5. /* 算法参考：http://www.ab126.com/shuxue/3820.html */
   
6. 
   
7. double sinx(double x)
   
8. {
   
9.         long long b , i , j , s                      ;
   
10.         double e , v                                 ;
    
11.         for(e = 1.0 , s = 1 , v = 0.0 , i = 1 ; e >= 0.00001 ; i += 2 , s = - s) {
    
12.                 for(b = 1 , j = i ; j ; j --) b *= j ;
    
13.                 e = pow(x , i) / b                   ;
    
14.                 v += s * e                           ;
    
15.         }
    
16.         return v                                     ;
    
17. }
    
18. 
    
19. double cosx(double x)
    
20. {
    
21.         long long b , i , j , s                      ;
    
22.         double e , v                                 ;
    
23.         for(e = v = 1.0 , s = -1 , i = 2 ; e >= 0.00001 ; i += 2 , s = - s) {
    
24.                 for(b = 1 , j = i ; j ; j --) b *= j ;
    
25.                 e = pow(x , i) / b                   ;
    
26.                 v += s * e                           ;
    
27.         }
    
28.         return v                                     ;
    
29. }
    
30. 
    
31. int main(void)
    
32. {
    
33.         double x = 0.0                               ;
    
34.         printf("sin(%lf) = %lf\n" , x , sinx(x))     ;
    
35.         printf("cos(%lf) = %lf\n" , x , cosx(x))     ;
    
36.         x = M_PI / 6                                 ;
    
37.         printf("sin(%lf) = %lf\n" , x , sinx(x))     ;
    
38.         printf("cos(%lf) = %lf\n" , x , cosx(x))     ;
    
39.         x = M_PI / 3                                 ;
    
40.         printf("sin(%lf) = %lf\n" , x , sinx(x))     ;
    
41.         printf("cos(%lf) = %lf\n" , x , cosx(x))     ;
    
42.         x = M_PI / 2                                 ;
    
43.         printf("sin(%lf) = %lf\n" , x , sinx(x))     ;
    
44.         printf("cos(%lf) = %lf\n" , x , cosx(x))     ;
    
45. }

复制代码

        编译、运行实况：

1. D:\00.Excise\C>g++ -o x x.c
   
2. 
   
3. D:\00.Excise\C>x
   
4. sin(0.000000) = 0.000000
   
5. cos(0.000000) = 1.000000
   
6. sin(0.523599) = 0.500000
   
7. cos(0.523599) = 0.866025
   
8. sin(1.047198) = 0.866025
   
9. cos(1.047198) = 0.500000
   
10. sin(1.570796) = 1.000000
    
11. cos(1.570796) = 0.000000
    
12. 
    
13. D:\00.Excise\C>

复制代码