Current Assay class simulates using forward Euler with fixed step size determined by user. We should use scipy.integrate.odeint for simulating genetic networks. This should improve efficiency and accuracy.