function selkov_bendixson( )
    syms x y a b ; 
    f =  ; % ODEを線形化 
    g =  ; % ODEを線形化 
    J =  ; % f，gからヤコビ行列Jをつくる 
      
    J = subs(J, [a,b] ,[0.06 , 0.6]) ; % パラメータa，bに値を代入 
    B =  ; % 対角要素の和 

    figure(1) 
    hold on ; 
    for i=0:0.5:3 % xの値を振る 
        for j=0:1:10 % yの値を振る 
            B_value = subs( ) ; % 式Bのx，yにi，jを代入 
            if( B_value > 0 ) 
                plot(i,j,'ko','MarkerFaceColor','w') ; 
            else 
                plot(i,j,'ko','MarkerFaceColor','k') ; 
            end 
        end 
    end 
end 
