function ozbudak( )

    time = 0:0.1:30;

    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)
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