function selkov()
    % 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'); % 横軸:時間、縦軸:濃度でプロット
    legend({'ADP','F6P'});

end

function dsdt = ODE(t, s, param)
    x = s(1); % ADP
    y = s(2); % F6P
    a = param(1);
    b = param(2);
    
    v1 = 　　　　　　　　　; % 空欄を埋めよ
    v2 = 　　　　　　　　　; % 空欄を埋めよ
    v3 = 　　　　　　　　　; % 空欄を埋めよ
    
    dsdt(1, :) = 　　　　　; % 空欄を埋めよ % ADPの時間変化
    dsdt(2, :) = 　　　　　; % 空欄を埋めよ % F6Pの時間変化
end

