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
It would also be useful to add a latex representation of the formula, so that it renders nicely in Sphinx. Look at some of the special functions in the functions module to see how this is done.