function task2_1()    
    param = 　【ここを埋める】; % パラメータa,b,K,h,Iの値
    x_ini = [1.0,1.0]; % x1,x2 の初期値
    
    time = 0:0.1:10; % シミュレーションを行う時間区間

    for x1_ini = 0.5:0.25:1.5 % x1の初期値を0.5から1.5まで0.25刻みで変動
        x_ini(1) = x1_ini; 
        [t, time_course] = ode45(@(t, x) ODE(t,x,param), time, x_ini);
        figure; 
        hold on;
        plot(t, time_course(:,1), 'r'); % 横軸:時間、縦軸:x1でplot
        plot(t, time_course(:,2), 'b'); % 横軸:時間、縦軸:x2でplot
        xlabel('Time (t)'); % 横軸のラベル名
        ylabel('Delta expression'); % 縦軸のラベル名
        legend({'Cell 1 (x1)','Cell 2 (x2)'}); % 各系列の説明
        title(join(["x1(0) = ",x1_ini])); % タイトルのラベル名
    end

end

function dxdt = ODE(~,x,param) % ODE(t,y,param)という関数の定義


　【ここを埋める】

end