function cmpr()

LP_Constants;

load('simu/MHD13.dat');
load('simu/DSMC13.dat');

mean(MHD13(:,5))
mean(DSMC13(:,9))

MT = MHD13(:,12)./((MHD13(:,5) + DSMC13(:,9)).*CO.Kb1).*CO.Kb2;
nT = DSMC13(:,8).*CO.Kb2;
meanT = (MT + nT)./2;
weightT = (MHD13(:,5).*MT+DSMC13(:,9).*nT)./(MHD13(:,5)+DSMC13(:,9));

TH = linspace(0,MHD13(360,1),360);

subplot(4,1,1), plot(TH,MT)
title('T_{MHD}  utraknat fran MHD-modell med antagandet T_{e} = T_{i} = T_{n}')
xlabel('Distance (m)')
ylabel('T (eV)')
hold on
subplot(4,1,2), plot(TH,nT)
title('T_{n}  rakt fran modell')
xlabel('Distance (m)')
ylabel('T (eV)')
subplot(4,1,3), plot(TH,meanT)
title('(T_{MHD}+T_{n}) / 2')
xlabel('Distance (m)')
ylabel('T (eV)')
subplot(4,1,4), plot(TH,weightT)
title('WeightT')
xlabel('Distance (m)')
ylabel('T (eV)')