Thanks for the pointer to diff, we still need to do 
```
powers = np.diff(powers_, axis=1)
powers -= 1
powers[0, :] += 1
```