| |
| /* An example that finds all primes between 2 and limit. */ |
| |
| limit = 2047; |
| |
| /* auto prime[], num, p, root, i */ |
| |
| prime[1] = 2; |
| prime[2] = 3; |
| num = 2; |
| if (limit >= 2) print "prime 1 = 2\n" |
| if (limit >= 3) print "prime 2 = 3\n"; |
| scale = 0; |
| for (timeloop = 0; timeloop < 30; timeloop += 1) { |
| num = 2; |
| for ( p=5; p <= limit; p += 2) { |
| root = sqrt(p); |
| isprime = 1; |
| for ( i = 1; i < num && prime[i] <= root; i++ ) { |
| if ( p % prime[i] == 0 ) { |
| isprime = 0; |
| break; |
| } |
| } |
| if (isprime) { |
| num += 1; |
| prime [num] = p; |
| print "prime ", num, " = ", p, "\n" |
| } |
| } |
| } |