function ozbudak_hysteresis_timeseries
	
	x_init = 0; % TMG (something like IPTG)
	y_init = 0; % Fluorescence of tagged LacY permease
	
	alpha = 0.1;
	beta = 10;
	rho = 25;
	x_out = 5;
	
	
	time = 0:0.5:60;
	param = [ alpha , beta , rho , x_out ]; % Constants
	s0 = [ x_init , y_init ];    % Initial values
	
	[t1,time_course1] = ode15s(@(t,s) ODE(t,s,param),time,s0);
	
	time = 60:0.5:120;
	x_out =
	x_init =
	y_init =
	s0 = [ x_init , y_init ];    % Initial values
	param = [ alpha , beta , rho , x_out ]; % Constants
	
	[t2,time_course2] = ode15s(@(t,s) ODE(t,s,param),time,s0);
	
	figure(1)
	plot([t1; t2],[time_course1(:,2); time_course2(:,2)])
	hold on
end


function dsdt = ODE(t,s,param)
	x = s(1); % TMG (something like IPTG)
	y = s(2); % Fluorescence of tagged LacY permease
	alpha = param(1);
	beta = param(2);
	rho = param(3);
	x_out = param(4);
	
	dsdt(1,:) =
	dsdt(2,:) =
end