
clear				% clear previous definitions

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

xlab = 'Weight (kg)'; 		% X-label
ylab = 'Metabolism (kcal)'; 			% Y-label
mytitle = ''; 	% Title

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

%setnextdir;
%load(['data',nextdir,'kleiber.data']);
xxkleiber = [0.021	0.282	0.41	2.98	1.52	2.46	3.57	4.33 ...	
    5.33	3	4.2	6.6	14.1	24.8	23.6	36	38	46.4	46.8 ...	
    57.2	54.8	57.9	300	435	482	600  ];
yykleiber = [3.6	28.1	35.1	167	83	119	154	191	233	152	207	288	...
    534	875	872	800	1090	1254	1330	1368	1224	1320 ...	
    4221	8166	7754	7877];

x1=0.01;
x2=700;

fig_def1;			% start of defaults and definitions


yy = 66.82*xx.^0.7565;
yy1 = 63.86*xx.^0.7685;


%loglog(xx,yy,'b-','MarkerSize',ms,'LineWidth',lw);
plot(xx,yy,'b-','MarkerSize',ms,'LineWidth',lw);
hold on
plot(xx,yy1,'r-','MarkerSize',ms,'LineWidth',lw);
%loglog(xxkleiber,yykleiber,'bo','MarkerSize',ms,'LineWidth',lw);
plot(xxkleiber,yykleiber,'bo','MarkerSize',ms,'LineWidth',lw);
%plot([0 500],[0 500],'-','MarkerSize',ms,'LineWidth',lw/1.75,'Color',darkgreen);

grid
%legend('Logistic Model','Beverton-Holt','Hassell''s Model','Beetle Data','Identity Map' ,4);
%text(5,80,'h(a) = 6.45a + 73.9','FontSize',fs_text);

%axis square
ax = [0 700 0 10000];
fig_def2;			% end of defaults and definitions

fontlabs = 'Times New Roman';   % Font type used in labels
xlabel(xlab,'FontSize',14,'FontName',fontlabs,'interpreter','latex');  % x-Label size and font
ylabel(ylab,'FontSize',14,'FontName',fontlabs,'interpreter','latex');  % y-Label size and font
title(mytitle,'FontSize',16,'FontName','Times New Roman','interpreter','latex'); % Title size/font
set(gca,'FontSize',12);         % Axis tick font size

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Here we print the figure in postscript
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% print -depsc kleiber_pow.eps		% Color
% print -depsc kleiber_pow2.eps
print -depsc kleiber_log.eps

clear