function dx = twoOsc(t,x)

%
dx = zeros(4,1);
g = 0.05;
dx(1) = x(2)-g*(x(1)-x(3));
dx(2) = -1*x(1);
dx(3) = x(4)-g*(x(3)-x(1));
dx(4) = -1*x(3);
%}


