temp=-2:.1:40;

q10_eppley=1.88;
q10_anderson=1.46;

b_eppley=log(q10_eppley)/10;
b_anderson=log(q10_anderson)/10;

Ae=-4000;

tempo=20;
tkel=273.15;

onet=ones(size(temp));

eppley=onet.*exp(b_eppley*(temp-tempo));
ander=onet.*exp(b_anderson*(temp-tempo));
arren=onet.*exp(Ae*((onet./(temp+tkel))-(onet./(tempo+tkel))));

figure
plot(temp,eppley,'k--','LineWidth',2);
hold on
plot(temp,ander,'b-','LineWidth',2);
plot(temp,arren,'r-','LineWidth',2);
bolder
axis([-2 35 0 2.5]);
xlabel('temperature (degC)'); ylabel('growth function (unitless)')
