
clear				% clear previous definitions

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

xlab = 't'; 		% X-label
ylab = 'Probability'; 			% Y-label
mytitle = 'Stochastic Birth Model'; 	% Title

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

x1=0;
x2=10;

fig_def1;			% start of defaults and definitions

p1 = exp(-0.5*xx);
p2 = (50*exp(0.01*xx)-50).*exp(-0.51*xx);
p3 = (1275*exp(0.01*xx).^2-2550*exp(0.01*xx)+1275).*exp(-0.52*xx);
p4 = (22100*exp(0.01*xx).^3-66300*exp(0.01*xx).^2+66300*exp(0.01*xx)-22100).*exp(-0.53*xx);

plot(xx,p1,'b-');
hold on;
plot(xx,p2,'r-');
plot(xx,p3,'g-');
plot(xx,p4,'k-');

grid

ax=[0 10 0 1];
legend('P_{50}(t)','P_{51}(t)','P_{52}(t)','P_{53}(t)',1);
%text(20,0.3,'Removed','FontSize',fs_text);

fig_def2;			% end of defaults and definitions

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Here we print the figure in postscript
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
print -depsc birth_prob.eps		% Color

clear