dx(t)/dt = x(t) -beta*y(t); dy(t)/dt = gamma*x;