麦克劳林公式求sin(x)

无意中发现的一条求sin公式，总结之。

测试平台：Win7 64b + VS2012

克劳林公式：

算法描述：

 1 #include <iostream>
 2 #include <cmath>
 3 using namespace std;
 4 
 5 inline int factorial(int x)
 6 {
 7     if ( x==1 )
 8         return 1;
 9     return factorial(x-1)*x;
10 }
11 
12 inline double square (double x, int n)
13 {
14     double ret=1.0;
15 
16     for (int i = 0; i<n; i++)
17         ret *= x;
18 
19     return ret;
20 }
21 
22 inline double _sin(double x)
23 {
24     double y = x/180*3.1415926;
25     return ( y - 1.0/factorial(3) * square(y, 3) + 1.0/factorial(5)*square(y, 5) );
26 }
27 
28 int main(void)
29 {
30     int x = 30;
31 
32     cout << "sin (30) :" << endl;
33     cout << _sin (30) << endl;
34     cout << sin (30* 3.1415926/180) << endl;
35 
36     cout << "sin (60) :" << endl;
37     cout << _sin(60) << endl;
38     cout << sin (60 * 3.1415926/180) << endl;
39 
40     cout << "sin (90) :" << endl;
41     cout << _sin(90) << endl;
42     cout << sin(90*3.1415923/180) << endl;
43     return 0;
44 }

测试结果：

sin (30) :
0.500002
0.5

sin (60) :
0.866295
0.866025

sin (90) :
1.00452
1