三角関数 S2 General Cosine
Closed; even; self setting; 34 storage locations;
time:37.6 msec; working spaces 0gL, 2gL.
Accuracy: 2×10-11
Places: 1/8 cosπx in Acc.
where x=C(0oL)
Notes: 1. The result will be correct regardless of overflow in the
argument, since the addition or subtraction is always
correct modulo 2.
2. 1/8 sinπx can be computed by subtracting 1/2 from the
argument to be placed in 0oL. The constant 1/2 can be
obtained from 21st storage loaction in this subroutine.
3. The following is an example of a master routine employing
the 3/4 binary point system, i.e. the binary point being
placed between d8 and d9.
pl x ┐ shift binary point
l 8 ┘
s 21s -1/2
tl o ┐
it │ jump to s
jl s ┘
r 5 adjust binary pont
0s: a 40 ┐ plant link
1 x 20r ┘
2 plo ┐
3 kl 5r │ take modulus
4 nlo ┘
5 s 21r -1/2
6 tlg
7 vlg
8 tl 2g
9 vl 32r ┐
10 sl 30r │
11 vl 2g │
12 al 28r │
13 vl 2g │ summation of
14 sl 26r │ power series
15 vl 2g │
16 al 24r │
17 vl 2g │
18 sl 22r │
19 vlg ┘
20 jl ( ) link
21 | e =1/2
22 | 51471 ┐ a1 ┐
23 | 223880 ┘ │
24 | 84667 ┐ a3 │
25 | 211265 ┘ │
26 | 41781 ┐ a5 │constants of
27 | 229407 ┘ │power series
28 | 9818 ┐ a7 │
29 | 36551 ┘ │
30 | 1344 ┐ a9 │
31 | 71435 ┘ │
32 | 115 ┐ a11 │
33 | 20532 ┘ ┘