Add monadic variants for iterateN, unfoldr, unfoldrN

I needed such functions while working on some toy Monte-Carlo simulations. It turned out that all required stream fusion machniery for these functions was present already but they weren't implemented for some reason

There're no tests. I don't really know what to test here.