function prob28

n1 = 1;
sum = 0;

for (i=1:1001)
  sum = sum + n1;
  df = floor(i/2) * 2;
  n2 = n1 + df;
  n1 = n1 + 2*i;
  sum = sum + n2;
end

sum - 1