You can easily model something like this with regular old ordinary differential equations in Mathematica or MATLAB. The hard part would be gathering the data to design realistic equations.