weighted_moments function

Compute weighted moments