% Orbit diagram

global alfa
figure; hold;
for alfa=2.8:0.01:4
   y0=[0.1]; n=300; [t,y]=eulere('quadr',0,n,y0,n);
   plot(alfa,y(200:300),'linewidth',2);
end
axis([2.8 4 0 1]);
