function bifurcation1D()
    rmin = -0.05; % rの最小値
    rmax = 0.05; % rの最大値
    dr = (rmax-rmin)/10; % rの刻み幅
    xmin = -0.05; % xの最小値
    xmax = 0.05; % xの最大値
    dx = (xmax-xmin)/10; % xの刻み幅
    [r,x] = meshgrid(rmin:dr:rmax, xmin:dx:xmax); % rとxの格子点を作成
    Vr = zeros(size(r)); % ベクトルの横成分はゼロ
    Vx = r.*x; % ベクトルの縦成分 ＝ 微分方程式の右辺
    figure; % 図を出す
    hold on; % 重ね書きOK
    quiver(r,x,Vr,Vx); % ベクトル場の作成
    fimplicit(@(r,x) r.*x, [rmin, rmax, xmin, xmax]); % 平衡点の計算
    xlabel('Parameter (r)'); % x軸のラベル名
    ylabel('State variable (x)'); % y軸のラベル名
end