function bifur(a0, a1, n)
  figure;
  hold on;
  delta = (a1-a0)/n;
  for i=0:n,
      a=a0+i*delta;
      e=stab_eq(a);
      plot(real(e),imag(e),'o');
  end
