clear all;
close all;




load datataui;

strtop=['{\bf RTC for the linear coupled model}'];
strmed=['   '];
strbot=[...
'\nu=',num2str(nu),...
', \alpha=',num2str(alpha),...
', \beta=',num2str(beta),...
', \tau_{LGN} =',num2str(T),...%', \sigma_{LGN}=',num2str(sigma),...
', \tau_E=',num2str(taue),...
', \sigma_E=',num2str(ae),...
', \tau_I=',num2str(taui),...
', \sigma_I=',num2str(ai),...
', s_{EE}=',num2str(see),...
', s_{EI}=',num2str(sei),...
', s_{IE}=',num2str(sie)...
', s_{II}=',num2str(sii)...
];


paramtE=parmin:dpar:parmin+86*dpar;
for k=1:87
  mineE(k,1)=mine(k,1);
end;

paramtI=parmin:dpar:parmin+57*dpar;
for k=1:58
  miniI(k,1)=mini(k,1);
end;

figure(1);
axes('fontsize',24);
plot(paramtE,mineE(:,1),'k','LineWidth',4);
axis([0 2.75 0.3 1.3]);
hold;
plot(paramtI,miniI(:,1),'-- k','LineWidth',4);
xlabel('\tau_{\it I}','fontname','times', 'fontsize', 34);
ylabel('{\it t}_{min}','fontname','times', 'fontsize', 34);
legend('\fontname{times}\fontsize{34}{\it t}_{{\it E},min}',...
'\fontname{times}\fontsize{34}{\it t}_{{\it I},min}','Location','NorthWest');
%title({strtop,strmed,strbot});