function task3_3( )   
	figure(1);
	hold on;
	xlabel('Time (t)'); % 横軸のラベル名
	ylabel('LacY (y)'); % 縦軸のラベル名
	figure(2);
	hold on;
	xlabel('Time (t)'); % 横軸のラベル名
	ylabel('TMGout (Xout)'); % 縦軸のラベル名
    c = colormap(jet); % 表示用のカラーマップ

    【ここを埋める】; % パラメータ
    for 【ここを埋める】 % x_outの前半の値を変化させる
        time = 0:1:60; % シミュレーション前半の時間区間
        x_init = 0; y_init = 0; % 前半部の初期値
        s0 = [ x_init , y_init ];    % 初期値ベクトル
        param = [ alpha , beta , rho , x_out_before ]; % パラメータ・ベクトル
        [t1,time_course1] = ode15s(@(t,s) ODE(t,s,param),time,s0); % ODE走らせる

        time = 60:1:120; % シミュレーション後半の時間区間
        x_init = time_course1(end,1); y_init = time_course1(end,2); % 後半部の初期値
        s0 = [ x_init , y_init ];    % 初期値ベクトル
        【ここを埋める】; % x_outの後半の値（固定）
        param = [ alpha , beta , rho , x_out_after ]; % パラメータ・ベクトル
        [t2,time_course2] = ode15s(@(t,s) ODE(t,s,param),time,s0);  % ODE走らせる
        
        figure(1); % LacYのプロット
    	plot([t1; t2],[time_course1(:,2); time_course2(:,2)], 'linewidth', 2,  'color',c(fix(x_out_before/10*64),:));
        figure(2); % TMGoutのプロット
        plot([t1; t2],[ones(length(t1))*x_out_before; ones(length(t2))*x_out_after], 'linewidth', 2, 'color',c(fix(x_out_before/10*64),:));
    end
end

function dsdt = ODE(~,s,param)
　【ここを埋める】
end