for d in xrange(m.ndim + 1, m.ndim * 3):
        x = mu.expand_dims(np.copy(m), d)
        assert_equals(x.ndim, d)
        assert_equals(x.shape[x.ndim - m.ndim:], m.shape)
        assert_equals(x.shape[:x.ndim - m.ndim], (1,) * (x.ndim - m.ndim))