n = 5;
ps = 0:n-1;
table = mod(ps'*ps,n)
sn = sum(sum(table))
lowerbound = 0.5 - (1/(sqrt(n)))
value = sn/(n*n*n)
upperbound = 0.5