function [ dXdt ] = toggle( t, x )

% toggle dynamics 

% degradation of X 
dx = 0.2;

% degradation of Y 
dy = 0.2;

% Dyamics for LacI mRNA
dXdt(1,:) = 4/(4+x(2)^2) - dx*x(1);

% Dyamics for LacI protein 
dXdt(2,:) = 4/(4+x(1)^2) - dy*x(2);

end