Sum[k^2/(2 k - 1) - k^2/(2 k + 1), {k, 1, 1001}]