ENH: interpolate: allow polyint evaluation on arbitrary shape arrays

Directly from ticket http://projects.scipy.org/scipy/ticket/1021

I didn't read this through or really test it, but the tests seem to pass, and it's from the original author of the module.