%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% maxwell.m:
% Oving 5 FY1005/TFY4165 Termisk fysikk 2013:
% Maxwells hastighetsfordeling for skiver paa luftputebord.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Nullstiller variable og lukker figurer
clear all;
close all;
% Leser eksperimentelle data fra 49 filer mass8_1.txt, mass8_2.txt osv 
% Hver fil inneholder 3 kolonner: 
%      tid (i sekunder), horisontal posisjon, vertikal posisjon
% Posisjonsverdiene multipliseres med faktoren 0.22 for
% at de skal angi verdier i cm.
% Det er 2 linjer med overskrift, derfor 'headerlines',2
% Heltallet j holder styr paa totalt antall hastighetsmaalinger
j=0;
% Tidssteget er 1/15 sekund
dt = 1.0/15.0;
for ifil = 1:49
   filnavn = ['mass8_',int2str(ifil),'.txt'];
   [t,x,y] = textread(filnavn,'%f %f %f','headerlines',2);
   x = x*0.22;
   y = y*0.22;
   n = length(t);
   for i = 1:n-1, % For i=n kan det ikke beregnes forflytning.
      j = j+1;
      dx(j) = x(i+1) - x(i);
      dy(j) = y(i+1) - y(i);
   end
end
% Beregner forflytningene (absoluttverdi)
ds = ..........;
% Beregner fart v
v  = ..........;
%Sorterer v med tanke paa aa lage histogram
v = sort(v);
% binantall = antall histogram-intervaller for v
binantall = ..........;
% beregner intervallbredde
binbredde = ...........;
% nullstiller tabell for v_fordeling
fordeling = zeros(1,binantall);
k = 1;
ibin = 1;
while (ibin <= binantall)
    while (k <= j && v(k) <= binbredde*ibin)
        fordeling(ibin) = fordeling(ibin)+1;
        k=k+1;
    end
    vbin(ibin) = binbredde*ibin;
    ibin=ibin+1;
end
% Plotting
figure;
bar(vbin,fordeling);
xlabel('v (cm/s)');
ylabel('Antall fartsmålinger');
title('Fartsfordeling');
