(PARI) is(n)=if(kronecker(17, n)>0 && kronecker(-17, n)>0 && n>2 && isprime(n), kronecker(lift((1+sqrt(Mod(17, n)))/2), n)>0, n==17) \\ Charles R Greathouse IV, Nov 11 2012
(Magma) /* By first comment: */ [17] cat [p: p in PrimesInInterval(3, 2200) | LegendreSymbol(-17, p) eq 1 and exists{x: x in ResidueClassRing(p) | x^4-x^2 eq 4}]; // Bruno Berselli, Nov 11 2012