Merge from vendor branch LIBSTDC++:
[dragonfly.git] / contrib / bc / Examples / primes.b
1
2 /* An example that finds all primes between 2 and limit. */
3
4 define primes (limit) {
5     auto num, p, root, i
6
7     prime[1] = 2;
8     prime[2] = 3;
9     num = 2;
10     if (limit >= 2) print "prime 1 = 2\n"
11     if (limit >= 3) print "prime 2 = 3\n";
12     scale = 0;
13
14     for ( p=5; p <= limit; p += 2)  {
15         root = sqrt(p);
16         isprime = 1;
17         for ( i = 1;  i < num && prime[i] <= root; i++ ) {
18             if ( p % prime[i] == 0 ) {
19                 isprime = 0;
20                 break;
21             }
22         }
23         if (isprime) {
24             num += 1;
25             prime [num] = p;
26             print "prime ", num, " = ", p, "\n"
27         }
28      }
29 }
30
31
32 print "\ntyping 'primes (10)' will print all primes less than 10.\n"