integrate with scipy.odeint package, shorten alpha/betas