function f(n)
    x = randn(n)
    norm(x-mean(x))^2
end

n=11
mean([f(n) for i=1:1_000_000])

n=5
mean([f(n) for i=1:1_000_000])

# example 
n = 4
P = eye(Int,n) - 1//n