bialko.m

function v = bialko(~,x)
   global alfa;
   global beta;
   global gamma;
   global miu;
   global h;

   v = zeros(2,1);
   v(1) = (alfa/(1+x(2)**h))-(miu*x(1));
   v(2) = (beta*x(1))-(gamma*x(2));
end