%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% friksjon.m:
% Oving 3 Fysikk 2013:
% Plotting av eksperimentelle data for anvendt
% snorkraft med snor rundt plastsylinder.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Nullstiller alle variable
clear;      
% phi = kontaktvinkel (radianer)
% S = snordrag/g (gram)
phi = [1.57080,3.14159,4.71239,6.28319,7.85398,9.42478,10.99557,12.56637];
S = [240,300,440,600,800,1000,1100,1400];
% S0 = 185 gram = loddets masse
S0 = 185;
% Funksjonen length(tabell) returnerer antall elementer i tabell
imax = length(phi);
% Regner ut ln S/S0:
y = log(S./S0);
% Regner ut mu-verdier:
mu = y./phi;
musum = sum(mu);
% Middelverdi av mu:
mumiddel = musum/imax;
kvadratavvik = sum((mu-mumiddel).^2);
% Standardavvik og standardfeil i mu:
deltamu = sqrt(kvadratavvik/(imax-1));
deltamumiddel = deltamu/sqrt(imax);
% Skriver ut middelverdi, standardavvik og standardfeil
fprintf('mumiddel = %3.4f\n',mumiddel);
fprintf('deltamu = %3.4f\n',deltamu);
fprintf('deltamumiddel = %3.4f\n',deltamumiddel);
% Plotter 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)');
saveas(1,'matlabfigur.pdf');