function task4_1()
    % selcovモデル
    % ODEを数値的に解く
    % パラメータの指定
    s0 = [1, 1]; % ADP(x)の初期値とF6P(y)の初期値 
    param = [0.06, 0.6]; % パラメータ a, b
    time = 0.01:0.1:100; %シミュレーションを行う時間
    
    % 数値計算
    [t, time_course] = ode45(@(t, s) ODE(t, s, param), time, s0); % ODEを解く
    
    % figure 1 に描く
    figure(1);
    plot(t, time_course(:,1), 'r', t, time_course(:,2), 'b'); % 横軸:時間、縦軸:濃度でプロット
end

function dsdt = ODE(~, s, param)
    x = s(1); % ADP
    y = s(2); % F6P
    a = param(1);
    b = param(2);
    
    v1 = b;
    v2 = a .* y + x.^2 * y;
    v3 = x;
    
    dsdt(1, :) = v2 - v3; % ADPの時間変化
    dsdt(2, :) = v1 - v2; % F6Pの時間変化
end