N:X:-between(3,N,X).
m:-2000:P,\+ (P:N,P*(P+2)mod\ -N<1;format('~d,~d
',[P,P+2])).