OFFSET
1,1
MAPLE
A094748 := proc(nmax) local a, diffs, d; a := [3] ; diffs := [] ; while nops(a) < nmax do d := 2 ; while true do if isprime(op(-1, a)-d) and not op(-1, a)-d in a and not d in diffs then a := [op(a), op(-1, a)-d] ; diffs := [op(diffs), d] ; break ; elif isprime(op(-1, a)+d) and not op(-1, a)+d in a and not d in diffs then a := [op(a), op(-1, a)+d] ; diffs := [op(diffs), d] ; break ; else d := d+2 ; fi ; od ; od; RETURN(diffs) ; end: A094748(100) ; # R. J. Mathar, Jul 29 2007
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, May 24 2004
EXTENSIONS
More terms from R. J. Mathar, Jul 29 2007
STATUS
approved