Added function `moment`

This is my first contributed function (i.e. please someone check thoroughly). I've checked the basic functionality against circstat (matlab), and added a number of tests. 

Questions and issues:

1. Currently the arguments `w` and `d` only alter the behaviour of this function if `cent=True`. Should they be doing something additionally in the body of the function?

1. I return the values in a tuple. I notice that most of our other functions return only single values, making them more compatible with e.g. the bootstrap decorator.

1. Further to this: does it make sense to include bootstrapping for the moments, and if so, how should this be done (a simple attempt broke because this function outputs a tuple).

Assuming this function is acceptable, I plan to do skewness and kurtosis over the weekend.
DO NOT MERGE THIS YET. I am writing more tests, and running a test with 2D data and centering over axis=1 breaks this line (AssertionError, dimensions do not match). So this is not yet right, however I'm trying to blind debug because I can't get anything to return around the decorators...