function ozbudak_hysteresis_timeseries( )   
    time = 0:0.1:60;

    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;

    param = [ alpha , beta , rho , x_out ]; % Constants
    s0 = [ x_init , y_init ];    % Initial values

    [t,time_course] = ode15s(@(t,s) ODE(t,s,param),time,s0); 

    figure(1)
    plot(t,time_course(:,2))
    hold on

    %
    % Change x_out
    %

    x_init = �@�@�@�@�@�@�@�@
    y_init = �@�@�@�@�@�@�@�@

    s0 = [ x_init , y_init ];    % Initial values

    x_out = �@�@�@�@�@�@�@�@
    param = [ alpha , beta , rho , x_out ]; % Constants

    time = 60:0.1:120;
    [t,time_course] = ode15s(@(t,s) ODE(t,s,param),time,s0); 

    figure(1)
    plot(t,time_course(:,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