function fitzhugh_nagumo()
    % FitzHugh-Nagumoモデル
    % ODEを数値的に解く
    s0 = [1, 0]; % 膜電位(v)の初期値と不活性化変数(w)の初期値 
    param = [0.0 , 0, 0, 10]; % パラメータ a, b, Input, tau
    time = 0.01:0.05:200; %シミュレーションを行う時間
    
    [t, time_course] = ode45(@(t, s) ODE(t, s, param), time, s0); % ODEを解く
    figure(1);
    plot(　　　　　　　　　　　　　　　　　　　　　　　　　　　　); % 横軸:時間、縦軸:変数でプロット
    figure(2);
    plot_nullcline(param);
    hold on;
    plot(　　　　　　　　　　　　　　　　　　　　　　　　　　　　); % 横軸:時間、縦軸:変数でプロット
    hold off;
    
end

function dsdt = ODE(t, s, param)
    v = s(1); % 膜電位
    w = s(2); % 不活性化変数
    a = param(1);
    b = param(2);
    Input = param(3);
    tau = param(4);
    
    dsdt(1, :) = 　　　　　　　　　　　　　; % 膜電位の時間変化
    dsdt(2, :) = 　　　　　　　　　　　　　; % 不活性化変数の時間変化
end


function plot_nullcline(param)
    a = param(1);
    b = param(2);
    Input = param(3);
    tau = param(4);
    v = -2:0.05:2;
    w1 = 　　　　　　　　　;
    w2 = 　　　　　　　　　;
    plot(v, w1, 'r', v, w2, 'm');
end