2012-01-14 22:13:02 j15(x)
The 15th spherical Bessel function of the first kind, j15(x).
The shocking thing is that I didn't mess this up and that it worked:

( c )
  1  fxBlock besselExplicit15 = ^(double x) {
2 return (29*(27*(25*(23*(21*(19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
3 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
4 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
5 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
6 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
7 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
8 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
9 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
10 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
11 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
12 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
13 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
14 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
15 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
16 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
17 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
18 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
19 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
20 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
21 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
22 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))))/x -
23 (17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
24 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
25 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
26 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
27 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
28 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
29 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
30 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
31 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
32 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
33 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
34 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
35 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))))/x -
36 (19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
37 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
38 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
39 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
40 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
41 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
42 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
43 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
44 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
45 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
46 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
47 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
48 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
49 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
50 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
51 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
52 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
53 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
54 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
55 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
56 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))))/x -
57 (21*(19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
58 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
59 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
60 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
61 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
62 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
63 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
64 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
65 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
66 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
67 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
68 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
69 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
70 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
71 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
72 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
73 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
74 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
75 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
76 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
77 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))))/x -
78 (17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
79 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
80 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
81 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
82 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
83 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
84 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
85 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
86 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
87 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
88 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
89 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
90 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))))/x -
91 (23*(21*(19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
92 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
93 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
94 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
95 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
96 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
97 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
98 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
99 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
100 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
101 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
102 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
103 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
104 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
105 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
106 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
107 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
108 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
109 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
110 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
111 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))))/x -
112 (17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
113 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
114 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
115 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
116 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
117 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
118 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
119 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
120 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
121 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
122 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
123 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
124 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))))/x -
125 (19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
126 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
127 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
128 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
129 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
130 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
131 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
132 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
133 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
134 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
135 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
136 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
137 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
138 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
139 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
140 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
141 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
142 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
143 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
144 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
145 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))))))/x -
146 (25*(23*(21*(19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
147 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
148 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
149 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
150 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
151 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
152 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
153 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
154 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
155 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
156 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
157 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
158 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
159 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
160 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
161 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
162 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
163 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
164 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
165 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
166 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))))/x -
167 (17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
168 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
169 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
170 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
171 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
172 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
173 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
174 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
175 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
176 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
177 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
178 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
179 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))))/x -
180 (19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
181 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
182 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
183 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
184 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
185 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
186 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
187 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
188 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
189 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
190 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
191 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
192 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
193 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
194 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
195 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
196 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
197 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
198 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
199 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
200 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))))/x -
201 (21*(19*(17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
202 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
203 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
204 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
205 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
206 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
207 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
208 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
209 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
210 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
211 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
212 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
213 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
214 (15*(13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
215 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
216 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))/x -
217 (9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
218 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))/x -
219 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
220 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
221 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))))))/x -
222 (17*(15*((13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
223 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
224 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
225 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
226 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x)))))/x -
227 (11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
228 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
229 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))))/x -
230 (13*(11*(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))
231 )/x-(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))
232 )/x-(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))
233 )/x-(9*(7*(5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))/x-(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x)))/x -
234 (5*(3*(sin(x)/pow(x,2)-cos(x)/x)/x-(sin(x)/x))/x-(sin(x)/pow(x,2)-cos(x)/x))))))));
235 };
  • Jeff (Sun, January 15th, 2012, 3:30pm UTC)
    For those wondering, yes this can be done easily as a recursive function. The goal that drove me to the above was to check on the overhead that occurs at runtime. Below is the C99 version of the Bessel functions of the first kind. The equation above runs about 3 times faster than the recursive version.

    ( c )
      1  __block fxBlock (^bessel)(int) = ^fxBlock (int n) {
    2 if (0==n) return ^(double x) { return sin(x) / x; };
    3 else if (1==n) return ^(double x) { return sin(x) / pow(x,2) — cos(x) / x; };
    4 else return ^(double x) { return (2*(n-1)+1)*(bessel(n-1))(x)/x — (bessel(n-2))(x); };
    5 };

Leave a comment