Fix #25: setdbspl and rms can use any ndim

The case where the number of dimensions is larger than 1 is treated as a
special case. If ndim > 1, then we create the necessary number of
`np.newaxis` so that the subtraction or division can be performed.

I'm not sure it's the optimal method, but it works for now.