Project Euler/28

function sumDiagonals(n):
    k = (n - 1) / 2
    sum = 1  # Center
    for m from 1 to k:
        sum += 16 * m * m + 4 * m + 4
    return sum

n = 1001
result = sumDiagonals(n)
print(result)