clear all close all clc time = [0:0.1:60]; global G delta T Ts k R gamma u w ks sw %% leaky integral control sw = 1; % sw == 0 -> w/o RNA decay % sw == 1 -> w/ RNA decay delta = 3; % half life 13 min T = 1; Ts = 4; ks = 50; R = 50; G = 1; % G = 1/\epsilon gamma = 1; w = 0.5; u = 1; y = T*u*ks/(Ts-T*u); %% leaky integral control sw = 0; % G = 1; k = 10*G; % k = complex removal rate [t,x_nd] = ode23s(@Fig2_sRNA,time,[0,0,0]); figure(1) plot(t,x_nd(:,3)/y); hold on sw = 1; for G = [1 10 100]; k = 10*G; [t,x] = ode23s(@Fig2_sRNA,time,[0,0,0]); figure(1) plot(t,x(:,3)/y); hold on end legend('\delta=0,T=T_s=1,\theta/\beta=10','\delta=3,T=T_s=1,\theta/\beta=10','\delta=3,T=T_s=10,\theta/\beta=10^2','\delta=3,T=T_s=10^2,\theta/\beta=10^3')