
clear				% clear previous definitions

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

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

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

x1=0;
x2=10;

fig_def1;			% start of defaults and definitions

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

y1 = sqrt(xx);
y2 = -sqrt(xx);
y3 = 0*xx;
y4 = 0*nxx;

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

grid on;

plot_arrow(5,-4,5,-2.5);
plot_arrow(5,4,5,2.5);
plot_arrow(5,0.5,5,2);
plot_arrow(5,-0.5,5,-2);

plot_arrow(-2.5,-2.5,-2.5,-1);
plot_arrow(-2.5,2.5,-2.5,1);

ax=[-5 10 -5 5];

%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 pitch.eps		% Color
%print -deps pitch.ps		% BW

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Here we preview the figure
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%!gs cricket.ps
%!gv cricket.ps &

clear