Quadrat special primes: Difference between revisions

=={{header|FreeBASIC}}==

<lang freebasic>

#include "isprime.bas"

dim as integer p = 2, j = 1

print 2;" ";

do

do

if isprime(p + j*j) then exit do

j += 1

loop

p += j*j

if p > 16000 then exit do

print p;" ";

j = 1

loop

print

</lang>

{{out}}

<pre>

2 3 7 11 47 83 227 263 587 911 947 983 1019 1163 1307 1451 1487 1523 1559 2459 3359 4259 4583 5483 5519 5843 5879 6203 6779 7103 7247 7283 7607 7643 8219 8363 10667 11243 11279 11423 12323 12647 12791 13367 13691 14591 14627 14771 15671

</pre>