ENH: replace exponentiation with cummulative product in  np.vander

Speeds calculation up by ~3x for 100x100 matrices, and by ~45x for
1000x1000