Post Reply 
(29C) Prime numbers up to 10'000
09-03-2018, 09:26 PM (This post was last modified: 09-04-2018 09:31 AM by Thomas Klemm.)
Post: #6
RE: (29C) Prime numbers up to 10'000
Listing for the HP-67:
Code:
001:         02:  2
002:      33 00:  STO 0     ; ∆n = 2
003:         07:  7
004:      33 01:  STO 1     ; n = 7
005:      33 03:  STO 3     ; π(3) = 7
006:         05:  5
007:      33 02:  STO 2     ; π(2) = 5
008:         04:  4
009:      35 33:  ST I      ; i = 4
010:   31 25 00:  LBL 0     ; for n loop
011:         06:  6         ; 6
012:      34 00:  RCL 0     ; ∆n    6
013:         51:  -         ; 6-∆n
014:      33 00:  STO 0     ; ∆n' = 6-∆n
015:   33 61 01:  STO+ 1    ; n' = n + ∆n
016:      34 01:  RCL 1     ; n
017:         05:  5         ; 5     n
018:         81:  ÷         ; n/5
019:      32 83:  FRAC      ; {n/5}
020:      31 51:  x=0       ; 5 | n ?
021:      22 00:  GTO 0     ; next n
022:      34 01:  RCL 1     ; n
023:         07:  7         ; 7     n
024:         81:  ÷         ; n/7
025:      32 83:  FRAC      ; {n/7}
026:      31 51:  x=0       ; 7 | n ?
027:      22 00:  GTO 0     ; next n
028:      34 01:  RCL 1     ; n
029:      33 24:  STO (i)   ; π(i) = n
030:      35 72:  PAUSE     ; is prime
031:      31 34:  ISZ       ; i = i+1
032:      35 34:  RC I      ; i
033:         02:  2
034:         05:  5         ; 25    i
035:      32 81:  x>y       ; 25 > i ?
036:      22 00:  GTO 0     ; next n
037:   31 25 01:  LBL 1     ; for n loop
038:         06:  6         ; 6
039:      34 00:  RCL 0     ; ∆n    6
040:         51:  -         ; 6-∆n
041:      33 00:  STO 0     ; ∆n' = 6-∆n
042:   33 61 01:  STO+ 1    ; n' = n + ∆n
043:         02:  2
044:      35 33:  ST I      ; i = 2
045:   31 25 02:  LBL 2     ; for i loop
046:      34 24:  RCL (i)   ; π(i)
047:      34 01:  RCL 1     ; n         π(i)
048:      34 24:  RCL (i)   ; π(i)      n     π(i)
049:         81:  ÷         ; n/π(i)    π(i)
050:      32 71:  x≤y       ; n ≤ π(i)² ?
051:      32 51:  x=y       ; n = π(i)² ?
052:      22 03:  GTO 3     ; n ≥ π(i)²
053:      34 01:  RCL 1     ; n
054:      35 72:  PAUSE     ; is prime
055:      22 01:  GTO 1     ; next n
056:   31 25 03:  LBL 3     ; n/π(i)
057:      32 83:  FRAC      ; {n/π(i)}
058:      31 51:  x=0       ; π(i) | n
059:      22 01:  GTO 1     ; next n
060:      31 34:  ISZ       ; i = i+1
061:      22 02:  GTO 2     ; next i

Registers:
Code:
00: ∆n
01: n
02: 5
03: 7
04: 11
05: 13
06: 17
07: 19
08: 23
09: 29
10: 31
11: 37
12: 41
13: 43
14: 47
15: 53
16: 59
17: 61
18: 67
19: 71
20: 73
21: 79
22: 83
23: 89
24: 97
25: i
Find all posts by this user
Quote this message in a reply
Post Reply 


Messages In This Thread
RE: (29C) Prime numbers up to 10'000 - Thomas Klemm - 09-03-2018 09:26 PM



User(s) browsing this thread: 1 Guest(s)