.ssc code for the HR 8799 system, orbits are set up to be face-on and circular (this is the likely configuration according to Marois et al. 2008) and mean anomalies are calculated using the positions given for 2008 Aug 12. Probably could do slightly better by doing a best fit to all of the observed positions but I'm leaving that for now. Gives quite a good match to the images shown in the press releases.
EDIT: the .ssc file is in the add-on in my next post.
Fomalhaut b would be trickier to model because the system is not being observed face-on, the orbit is likely eccentric and there doesn't seem to be a convenient table of values for the positions available.
