function task2_2()    
    param = 【ここを埋める】; % パラメータa,b,K,h,Iの値
    x_ini = [1.0,1.0]; % x1,x2 の初期値
    
    time = 0:1:200; % シミュレーションを行う時間区間
    figure;
    hold on;
    for 【ここを埋める】 % Iを変動させる
        param(5)= I;
        for x1_ini = 0.5:0.1:1.5 % x1の初期値を変動させる
            x_ini(1) = x1_ini; 
            [~, time_course] = ode45(@(t, x) ODE(t,x,param), time, x_ini);
            x1_end = time_course(end,1);
            plot(I, x1_end, 'r+', 'markersize', 10, 'linewidth', 2);
        end
    end
	xlabel('Cell-cell interaction (I)'); % 横軸のラベル名
	ylabel('Final state of cell 1 (x1)'); % 縦軸のラベル名

end

function dxdt = ODE(~,x,param) % ODE(t,y,param)という関数の定義
【ここを埋める】
end