divisor function sigma

Calculate the divisor function sigma_k(n) for positive integer n

sigma(n,0) is equal to divisor_count(n)
sigma(n,1) is equal to sum(divisors(n))
sigma(n,k) is equal to sum([x**k for x in divisors(n)])

Default for k is 1.

References
==========

- http://en.wikipedia.org/wiki/Divisor_function
Why did you add two underscores here? I think you only want one (or even zero), unless there is a good reason for this to be name mangled (and if there is, please note why in a comment).