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;