clear all close all clc time = [0:0.1:40]; global k1 k2 delta T k u K2 sw w % sw -> 1 w/ dilution; sw->0 w/o dilution delta = 1; T = 100; k = 500; w = 0.5; u = 20; K2 = 5; k1 = 10; k2 = 10; %% ideal integral control sw = 0; [t,x_nd] = ode23s(@Fig2_phospho,time,[0,0]); figure(1) plot(t,x_nd(:,2)/u) hold on %% leaky integral control with various gain sw = 1; for G = [1 10 10^2]; k1 = G*10; k2 = G*10; [t,x] = ode23s(@Fig2_phospho,time,[0,0]); figure(1) plot(t,x(:,2)/u) hold on end legend('\gamma=0,k_i/\gamma=10','\gamma=1,k_i/\gamma=10','\gamma=1,k_i/\gamma=10^2','\gamma=1,k_i/\gamma=10^3')