N = 1000000000

numbers = ones(N, 1);

for i = 2 : floor(sqrt(N))
  factor = i * (i + 1);
  multiples = factor : factor : N;
  numbers(multiples) = 0;
end

N - sum(numbers)