Primality by trial division: Difference between revisions

m
Line 766:
 
<lang SNOBOL4>define('isprime(n)i,max') :(isprime_end)
isprime isprime = n; max = sqrt(n); i = 1
le(n,1) :s(freturn)
eq(n,2) :s(return)
eq(remdr(n,2),0) :s(freturn)
max = sqrt(n); i = 1
isp1 i = le(i + 2,max) i + 2 :f(return)
eq(remdr(n,i),0) :s(freturn)f(isp1)
Anonymous user