%
%  g gives the RHS for the system 
%			x'(1) = g1(x,t)
%        x'(2) = g2(x,t)
%
%
function dx=g(t,x)

dx=[3*x(1)-x(2); x(1)-x(2)];

% lorenz equations
% dx = [10*(x(2)-x(1));28*x(1)-x(2)-x(1)*x(3);x(1)*x(2)-8/3*x(3)];

% linear system (stable star)
% dx=[-x(1); -x(2); -x(3)];

% linear system (unstable node)
% dx=[x(1); 2*x(2); 3*x(3)];

% linear system (unstable spiral and node )
% dx=[x(1)-x(2); x(1)+x(2); 2*x(3)];