Nightingale deviation of order p down an axis

[galenseilis]

Just like most functions in NumPy that act on a vector, it is desirable for the Nightingale deviation to be computed along a fiber of a multidimensional array (i.e. informally, a tensor). This can be accomplished by computing the function along an axis. Thus adding an axis parameter that is passed into subsequent functions is appropriate.