;
function osc3_derivative, t, x
 dxdt  = [-(!gamma_d - !gamma_A)*x[1] - !E*x[0] - !lvo*!E*x[4]/!Cpd, $
          !g*x[0]/x[2] - !E*(x[1]-0.) + !g*x[4], $
          -!gamma_A*x[1], $
          x[1],  $
          !gamma_l*x[1], $
          -!E*(x[5] - !wtot_A)]
 return, dxdt
end