+- HP Forums (https://www.hpmuseum.org/forum)
+-- Forum: HP Software Libraries (/forum-10.html)
+--- Forum: HP-65/67/97 Software Library (/forum-12.html)
+--- Thread: Euler's Constant (/thread-148.html)
Euler's Constant - Namir - 12-18-2013 05:58 AM
Euler's Constant
Gamma = H(n) - ln(n) + Poly(n)
Poly(n) = –1/(2*n) + 1/(12*n^2) –
1/(120*n^4) + 1/ (252*n^6) –
1/(240*n^8) + 1/(132*n^10)
HP-67 Implementation
Memory Map
R0 = n
R1 = i = 1 to n
R2 = sum of 1/i
R3 = Poly(n)
Implementation
Code:
1 LBL A
2 STO 0
3 STO 1
4 0
5 STO 2
6 LBL 0
7 RCL 1
8 1/X
9 STO+ 2
10 RCL 1
11 1
12 -
13 STO 1
14 X#0?
15 GTO 0
16 RCL 0
17 X^2
18 1/X
19 ENTER
20 ENTER
21 ENTER
22 132
23 1/X
24 *
25 2
26 4
27 0
28 CHS
29 GSB 2
30 252
31 GSB 2
32 120
33 CHS
34 GSB 2
35 1
36 2
37 GSB 2
38 RCL 0
39 1/X
40 2
41 /
42 -
43 RCL 0
44 LN
45 -
46 RCL 2
47 +
48 RTN
49 LBL 2
50 1/X
51 +
52 *
53 RTN