added diff eqs based on Steve's code

What's up with this tr indexing --- it appears to be hard coded, but I don't see what's going on with it conceptually or how I should code this into a model.