Oups! I forgot to add this routine in commit ac54f7f5. This routine provides a matlab implementation of the perfect foresight model solver.