close all
clear all
global tau1 tau0 arn AptoA

x0=[0 0 0];
tspan=[0,500];
[t,x]=ode45(@laser,tspan,x0)
plot(t,x)

%figure;
%plot(t,x(1),'k',t,x(2),'r',t,x(3),'g')
xlabel('\tau','FontSize',16)
ylabel('\nu_1, \nu_2, \rho','FontSize',16)
set(gca,'fontsize',16)
title('\nu_1, \nu_2, \rho vs. \tau','FontSize',16)
text(100,1.10,'\nu_1','Color','b','FontSize',16)      % 0.65
text(150,1.10,'\nu_2','Color','r','FontSize',16)      % 1.35
text(200,1.10,'\rho','Color',[1,0.5,0],'FontSize',16) % 0.45
annotation('textbox',...
    [0.6 0.54 0.25 0.25],... 
    'String',{['\tau_1 = ' num2str(tau1)], ['\tau_0 = ' num2str(tau0)], ['r_N = ' num2str(arn)], ['AptoA = ' num2str(AptoA)]},...
    'FontSize',12,...
    'FontName','Arial',...
    'LineStyle','--',...
    'EdgeColor',[1 1 0],...
    'LineWidth',2,...
    'BackgroundColor',[0.9 0.9 0.9],...
    'Color',[0.84 0.16 0]);