Judd and Guu's stochastic perturbation model in Python

I just uploaded a Python version of Judd and Guu's perturbation code to my website (code; notes).

If you happen to be a Python programmer, or a computational economist (even better: both!), then any feedback you can give on the code is highly appreciated. I'm fairly new to Python, so I'm sure I could have programmed some parts more efficiently.

Oh, and here is a picture of a bunch of Canberran kangaroos:

I thought you'd like to know.