%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% eksp_friksjon.m:
% Oving 3 TFY4104 Fysikk 2013:
% Plotting av eksperimentelle data for anvendt
% snorkraft med snor rundt plastsylinder.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Nullstiller alle variable
clear;      
% Leser eksperimentelle data fra filen eksp_friksjon.txt som inneholder
% 2 kolonner: phi, S/g i enheter hhv radianer, gram.
% Det er 1 linje med overskrift, derfor 'headerlines',1
[phi,S] = textread('eksp_friksjon.txt','%f %f','headerlines',1);
% Funksjonen length(tabell) returnerer antall elementer i tabell
lengde = length(phi);
% Regner ut ln S/S0:
S0 = S(1);
y = log(S./S0);
% Regner ut mu-verdier:
for i = 2:lengde,
    mu(i) = y(i)/phi(i);
end
% Middelverdi av mu:
mumiddel = 0.0;
for i = 2:lengde,
    mumiddel = mumiddel + mu(i);
end
mumiddel = mumiddel/(lengde-1);
% Standardavvik og standardfeil i mu:
deltamu = 0.0;
for i = 2:lengde,
    deltamu = deltamu + (mu(i)-mumiddel)^2;
end
deltamu = sqrt(deltamu/(lengde-2));
deltamumiddel = deltamu/sqrt(lengde-1);

% Visualiserer resultatene:

figure(1);
hold on;
scatter(phi,y,'k','*');
plot(phi,phi*mumiddel,'b');
plot(phi,phi*(mumiddel+deltamu),'g');
plot(phi,phi*(mumiddel-deltamu),'r');
hleg1 = legend('$\ln(S/S0)$','$\overline{\mu}\phi$', ...
'$(\overline{\mu}+\Delta\mu)\phi$','$(\overline{\mu}-\Delta\mu)\phi$');
set(hleg1,'location','northwest');
set(hleg1,'interpreter','latex');
title('Øving 3, Oppgave 1');
xlabel('\phi (rad)');
