;
function osc2_derivative, t, x
 dxdt  = [-(!gamma_d - !gamma_A)*x[1] - !E*x[0], $
          !g*x[0]/x[2] - !E*(x[1]-0.), $
          -!gamma_A*x[1], $
          x[1]]
 return, dxdt
end