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