ENH: add a weighted covariance calculation to numpy.function_base.

Weights can now be included in a covariance calculation in a similar manner as
the existing implementation of a weighted average.
Just a reminder for myself, but ddof may be missing here?