OFFSET
1,1
MAPLE
P:= select(isprime, [seq(i, i=3..20000, 2)]):
G:= P[2..-1]-P[1..-2]:
i0:= ListTools:-SelectFirst(j -> G[j]>=42, [$1..nops(G)]):
$ P[i0]+1 .. P[i0]+43; # Robert Israel, Jul 05 2017
MATHEMATICA
rPrimeNext[n_]:=Module[{k}, k=n+1; While[PrimeQ[k]||GCD[n, k]!=1, k++ ]; k]; a=1; lst={a}; Do[AppendTo[lst, a=rPrimeNext[a]], {n, 0, 2*7!}]; lst1=lst; q=41; lst={}; Do[If[lst1[[n+q]]-lst1[[n]]==q, AppendTo[lst, lst1[[n]]]], {n, 0, Length[lst1]-q}]; lst; a=FromDigits[lst]; Table[n, {n, a, a+q}]
CROSSREFS
KEYWORD
nonn,fini,full
AUTHOR
Vladimir Joseph Stephan Orlovsky, May 15 2010
EXTENSIONS
Edited by Robert Israel, Jul 05 2017
STATUS
approved