// file:    testPrimes.C
// author:  Robert Keller
// purpose: demonstrate use of Primes class
// compile: g++ testPrimes.C primes.o polya.o

#include "Primes.H"

Poly square(Poly x)
  {
  return x*x;
  }

main()
  {
  integer N = 1000;
  cout << "The first " << N << " primes: " << endl;
  Polylist primes = Primes::primes;
  cout << primes.prefix(N) << endl;

  cout << "Mapping over an infinite list: "
       << primes.map(square).prefix(10) << endl;
  }