ENH: Replace exponentiation with cumulative product in np.vander

Speeds calculation up by ~3x for 100x100 matrices, and by ~45x for
1000x1000
Merge pull request #4568 from jaimefrio/cumprod-vander.

ENH: Replace exponentiation with cumulative product in np.vander