dBdt function

Dynamics from EcoSim