%{ 
   This Matlab script plot was used to create a plot using LaTeX labelling.
   If the download does not work directly, you can probably cut-and-paste from
   here.
%}

clear all
close all

E = [ ...
   0.01   0.015   0.02   0.03   0.04   0.05   0.06   0.08 ...
   0.10   0.15    0.2    0.3    0.4    0.5    0.6    0.8  ...
   1.0    1.25    1.5    2.0    3.0    4.0    5.0    6.0    8.0    10. ];
mumu = [ ...
   5.21   1.6     0.778  0.371  0.267  0.225  0.205  0.185 ...
   0.171  0.151   0.137  0.119  0.106  0.0966 0.0894 0.0785 ...
   0.0707 0.0641  0.0575 0.0493 0.0396 0.0340 0.0303 0.0277 0.0243 0.0222];
mutr = [ ...
   4.79   1.28    0.512  0.149  0.0677 0.0418 0.0320 0.0262 ...
   0.0256 0.0277  0.0297 0.0319 0.0328 0.0330 0.0329 0.0321 ... 
   0.0311 0.02975 0.0284 0.0262 0.0229 0.0209 0.0195 0.0185 0.0170 0.0162];
muen = [ ...
   4.79   1.28    0.512  0.149  0.0677 0.0418 0.0320 0.0262 ...
   0.0256 0.0277  0.0297 0.0319 0.0328 0.0330 0.0329 0.0321 ... 
   0.0309 0.02955 0.0282 0.0260 0.0227 0.0206 0.0191 0.0180 0.0166 0.0157];

figure(1)
loglog(E,mumu,'k-',E,mutr,'k--',E,muen,'k-')
xlabel('$E_{\gamma}$ (MeV)','Interpreter','LaTex','FontSize',20)
ylabel('$\mu, \mu_{tr}, \mu_{en}$ for H$_2$O','Interpreter','LaTex','FontSize',20)
legendHandle = legend('$\mu$','$\mu_{tr}$','$\mu_{en}$');
set(legendHandle,'Interpreter','LaTex','FontSize',20)

figure(2)
loglog(E,1./mumu,'k-')
xlabel('$E_{\gamma}$ (MeV)','Interpreter','LaTex','FontSize',20)
ylabel('$1/\mu$ (cm)','Interpreter','LaTex','FontSize',20)
legendHandle = legend('mean free path $1/\mu$ in H$_2$O (cm)');
set(legendHandle,'Interpreter','LaTex','FontSize',16)

