figure
time=9; fs=10;
[T1,Y1] = ode23s('logisticCon',[0 time],[.1]);
[T2,Y2] = ode23s('Gompertz',[0 time],[.1, log(14)]);
[T3,Y3] = ode23s('regrowth',[0 time],[.1]);
x = 0:1:9;
y = 0.2*[0.5 3.5 5.5 7 6.5 7.2 6.5 7.3 6.8 7]';
option=odeset('AbsTol',1e-9,'RelTol',1e-9);
plot(T1(:,1),Y1(:,1),'r',T2(:,1),Y2(:,1),'b-.',T3(:,1),Y3(:,1),'k:',x,y,'-.gs','LineWidth',2,...
                'MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize', 4);
text(2,0.5,'logistic model: dx/dt=rx(1-x/K), r=ln(14), K=1.4, x(0)=.1','FontSize',fs);
text(2,0.3,'Gompertz model: dx/dt=ln(K/x(0))e^{-t}x, x(0)=.1','FontSize',fs);
text(2,0.1,'regrowth model: dx/dt=r(1-x/K), r=ln(14), K=1.4, x(0)=.1','FontSize',fs);
legend('logistic','Gompertz','regrowth',2)
title('Comparison of logistic, Gompertz and regrowth models','FontSize',12)
xlabel('time'); ylabel('x'); xlim([0 time]); ylim([0 1.8]);