(MAGMA)
//Program finds primes separated by an even number (called gap) which
//have a power of two between them. Program starts with the smallest
//power of two above gap. Primes less than this starting point can be
//checked by inspection.
gap:=16;
start:=Ilog2(gap)+1;
for i:= start to 1000 do
powerof2:=2^i;
for k:=powerof2gap+1 to powerof21 by 2 do
if (IsPrime(k) and IsPrime(k+gap)) then k;
end if;
end for;
end for;
