function task1_pulse3(model_type)
% model_type
%1: Ras
%2: Rap1

time = 1:100; 
pulse_duration_time_vec = 10:10:100;
colormat = jet(length(pulse_duration_time_vec));

figure;
hold on

if model_type ==1 % Ras model
    param = [0.5,5,0.0005,0.005,0.05,100]; % [k1,k2,k3,k4,k5,k6]
    % step刺激
    y0 = [1,0,0,0]; % [S,GEFa,GAPa,Rasa]
    [t,time_course] = ode15s(@(t,y) ODE_step_Ras(t,y,param),time,y0);
    plot(t,time_course(:,4), 'b');
    
    % pulse刺激
    for i = 1 : length(pulse_duration_time_vec)
        % pulse刺激の持続時間を順番に変える
        % 変数pulse_duration_time_vecに変え方が保存されている
        　　　
        　　　
        　　　
        plot(　,　, 'Color', colormat(i,:)); %Rasaの応答の時間変化を図示
        max_response(i) = 　; % Rasaの最大の応答を取得
    end
    
else % Rap1 model
    param = [0.5,0.5,0.2,5];
    % step刺激
    y0 = [1,0,0.005,0]; % [R,GEFa,GAPa,Rap1a]
    [t,time_course] = ode15s(@(t,y) ODE_step_Rap1(t,y,param),time,y0);
    plot(t,time_course(:,4), 'b');
    
    % pulse刺激
    for i = 1 : length(pulse_duration_time_vec)
        % pulse刺激の持続時間を順番に変える
        % 変数pulse_duration_time_vecに変え方が保存されている
        　　　
        　　　
        　　　
        plot(　,　, 'Color', colormat(i,:)); %Rap1aの応答の時間変化を図示
        max_response(i) = 　; % Rap1aの最大の応答を取得
    end
    
end

xlabel('time');
ylabel('response');
if model_type == 1
    title('Ras model')
    ylim([0,0.02])%y軸の範囲を修正
else
    title('Rap1 model')
    ylim([0,1])%y軸の範囲を修正
end

%% durationに対する最大応答強度の関係
figure;
plot(　, 　)
xlabel('duration time of pulse stimulation')
ylabel('max response')

%y軸の範囲を修正
if model_type == 1
    ylim([0,0.02])
else
    ylim([0,1])
end

end

%% Rasモデルでのstep刺激
function dydt = ODE_step_Ras(t,y,param)
S = y(1);
GEFa = y(2);
GAPa = y(3);
Rasa = y(4);

k1 = param(1);
k2 = param(2);
k3 = param(3);
k4 = param(4);
k5 = param(5);
k6 = param(6);

dydt(1,:) = 0;
dydt(2,:) = k1*S*(1-GEFa)-k2*GEFa;
dydt(3,:) = k3*S*(1-GAPa)-k4*GAPa;
dydt(4,:) = k5*GEFa*(1-Rasa)-k6*GAPa*Rasa;
end

%% Rasモデルでのpulse刺激
function dydt = ODE_pulse_Ras(t,y,param, pulse_duration_time)
% task1_pulse1.mのODE_pulse_Rasをコピペ
　　　　
　　　　
　　　　
end

%% Rap1モデルでのstep刺激
function dydt = ODE_step_Rap1(t,y,param)
S = y(1);
GEFa = y(2);
GAPa = y(3);
Rap1a = y(4);

k1 = param(1);
k2 = param(2);
k3 = param(3);
k4 = param(4);

dydt(1,:) = 0;
dydt(2,:) = k1*S*(1-GEFa)-k2*GEFa;
dydt(3,:) = 0;
dydt(4,:) = k3*GEFa*(1-Rap1a)-k4*GAPa*Rap1a;
end

%% Rap1モデルでのpulse刺激
function dydt = ODE_pulse_Rap1(t,y,param, pulse_duration_time)
% task1_pulse1.mのODE_pulse_Rap1をコピペ
　　　　
　　　　
　　　　
end

