Source code in the jq programming language

def emit_until(cond; stream): label $out | stream | if cond then break $out else . end;

# For 0 <= $n <= ., factorials[$n] is $n !
def factorials:
   reduce range(1; .+1) as $n ([1];
    .[$n] = $n * .[$n-1]);

def lpad($len): tostring | ($len - length) as $l | (" " * $l)[:$l] + .;

def primes: 2, (range(3; infinite; 2) | select(is_prime));

# Input: the limit of $p
def wilson_primes:
  def sgn: if . % 2 == 0 then 1 else -1 end;

  . as $limit
  | factorials as $facts
  | " n:  Wilson primes\n--------------------",
    (range(1;12) as $n
     | "\($n|lpad(2)) :  \(
       [emit_until( . >= $limit; primes)
        | select(. as $p
            | $p >= $n and
              (($facts[$n - 1] * $facts[$p - $n] - ($n|sgn)) 
               % ($p*$p) == 0 )) ])" );

11000 | wilson_primes