Added comments and the possibility to use a matlab implementation of the perfect foresight model solver.