function repressilator()
    % Repressilatorモデル
    % ODEを数値的に解く
    s0 = [1, 0.5, 0.0, 1, 1, 1]; % 各mRNA、タンパク質の初期値
    param = [100, 1, 5, 2]; % パラメータ alpha, alpha0, beta, n
    time = 0.01:0.05:50; %シミュレーションを行う時間
    %%
    [t, time_course] = ode45(@(t, s) ODE(t, s, param), time, s0); % ODEを解く
    figure(1);
    plot(　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　); %空欄を埋めよ
    legend({'TetR', 'LacI', 'lambda cI'})
    % 横軸:時間、縦軸:変数でプロット（3種類のタンパク質）
    figure(2);
    plot3(　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　); %空欄を埋めよ
    grid on;
    xlabel('TetR');
    ylabel('LacI');
    zlabel('lambda cI');
    % 横軸:TetR、縦軸:LacI、高さ軸: λcI
    
end

function dsdt = ODE(t, s, param)
    x = s(1); % tetR-lite
    y = s(2); % lacI-lite
    z = s(3); % λcI-lite
    X = s(4); % TetR
    Y = s(5); % LacI
    Z = s(6); % λcI
    
    alpha = param(1);
    alpha0 = param(2);
    beta = param(3);
    n = param(4);
    
    dsdt(1, :) = 　　　　　　　　　　　　　　; %空欄を埋めよ % tetR-lite
    dsdt(2, :) = 　　　　　　　　　　　　　　; %空欄を埋めよ % lacI-lite
    dsdt(3, :) = 　　　　　　　　　　　　　　; %空欄を埋めよ % λcI-lite
    dsdt(4, :) = 　　　　　　　　　　　　　　; %空欄を埋めよ % TetR
    dsdt(5, :) = 　　　　　　　　　　　　　　; %空欄を埋めよ % LacI
    dsdt(6, :) = 　　　　　　　　　　　　　　; %空欄を埋めよ % λcI
end
