
clear				% clear previous definitions

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% here you override any default:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

xlab = '\mu'; 		% X-label
ylab = 'y (Equilibria)'; 			% Y-label
mytitle = 'Transcritical Bifurcation'; 	% Title

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% here you plot your function(s):
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

x1=0;
x2=10;

nx1 = -10;
nx2 = 0;
nxx = linspace(nx1,nx2,500);

fig_def1;			% start of defaults and definitions

y1 = xx;
y2 = 0*xx;
y3 = 0*nxx;
y4 = nxx;

plot(xx,y1,'b-','MarkerSize',ms,'LineWidth',lw);
hold on
plot(xx,y2,'b--','MarkerSize',ms,'LineWidth',lw);
plot(nxx,y3,'b-','MarkerSize',ms,'LineWidth',lw);
plot(nxx,y4,'b--','MarkerSize',ms,'LineWidth',lw);

grid on;

plot_arrow(5,2,5,4);
plot_arrow(5,8,5,6);
plot_arrow(5,-2,5,-4);

plot_arrow(-5,4,-5,2);
plot_arrow(-5,-6,-5,-8);
plot_arrow(-5,-4,-5,-2);

ax=[-10 10 -10 10];

%legend('Susceptibles','Infected','Removed',2);
%text(25,0.85,'Susceptibles','FontSize',fs_text);

fig_def2;			% end of defaults and definitions

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Here we print the figure in postscript
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
print -depsc trans.eps		% Color
%print -deps trans.ps		% BW

clear