Package findprime

%% Read an integer n and produce a prime number that is
%% about n decimal digits long.

Import "math/prime".
Import "collect/string".
Import "misc/random".

Define nextPrime(?n) =
  head [: x | Let x = each(n _upto_ infinity).
              {prime?(x)} :]
%Define

Execute
  SetRandomSeed "323382873734054838928927378672165667473723837264646878".
  Writeln["How many decimal digits should the prime number have?"].
  Try
    Extract[$(?n:Natural)] from bxStdin.
    {n > 0}
  then
    Writeln[$nextPrime(randomRange(10^(n--1), 10^n))].
  else
    Writeln["I need a positive number."].
  %Try
%Execute

%Package