Fix the enumeration process of polynomial features

This is a fix for issue #3191 
I have to admit that I could not manage running nosetests in my virtualenv of my fork.

I tested the code in an external way, check it out [here](http://nbviewer.ipython.org/gist/aboSamoor/fb015e88c5a0152abd76)