9 lines
332 B
Text
9 lines
332 B
Text
CubanPrimes := proc(n) local i, cp;
|
|
cp := Array([]);
|
|
for i by 2 while numelems(cp) < n do
|
|
if isprime(3/4*i^2 + 1/4) then
|
|
ArrayTools:-Append(cp, 3/4*i^2 + 1/4);
|
|
end if;
|
|
end do;
|
|
return cp;
|
|
end proc;
|