
% Load All parameters

reactiondata;

% Set Critical parameters 

n0 = [30 10 0 0]; % initial population
S0 = 1e9*n0/NAV;       % initial concentration (nM)  
k = [0.6 0.1];    % rate constants
pars = {NAV,k,c,M,p,l,n0,S0,nu,tf,runs}; % parameters cell array

cells = stoch_gma(pars);
t_gilsp = cells{7};    
mean_a = cells{10}; 


figure
hplot = plot(t_gilsp,mean_a(:,1),'b-',t_gilsp,mean_a(:,2),'r-');
set(hplot,'LineWidth',1.5)
xlabel('time  (sec)','FontSize',8)
ylabel('mean of {\it\fontname{Times} a_\mu}','FontSize',8)

text(150,0.102,'{\it\fontname{Times} R_1}')
text(150,0.25,'{\it\fontname{Times} R_2}')

% Shifting the axis within the figure

posi = get(gca,'Position');
set(gca,'Position',posi.*[1.3 1.15 1 1],'FontSize',8);

% printing adjustments

set(gcf, 'PaperType', 'A4', 'PaperOrientation', 'portrait', 'PaperPositionMode', 'manual', 'PaperUnits', 'centimeters', 'PaperPosition', [2 2 8 6.3]);

textobj = findobj('type', 'text');
set(textobj, 'fontunits', 'points', 'fontsize', 8);

% Printing the figure into an .eps file:

print -depsc 'illu_a.eps'