%## Autor: Peter Bokes  <peter.bokes zo stuba.sk>
%## podporené z projektu KEGA 006STU-4/2022
%## 
%## Skript pre výpočet náboja na nabotých vodivých valčekoch
%## Dokumentácia: kf-lin.elf.stuba.sk:~bokes/public_html/DPPFI/F2-Coulomb-law-experiment.pdf
%##
function Coulomb
%##
%## manuál pre octave: https://docs.octave.org/interpreter/index.html#SEC_Contents
%## anoc nedokaze pracovat s pause() -- treba odkomentovat
%##
% hmotnosť jedného valčeka v kg
m = 0.042e-3
% gravitačné zýrchlenie v m.s^-2
g = 9.81                                                  
% dĺžka zavesu v m
l = 0.4                                                   
% konštanta vo vzorci
k = sqrt(4.0*pi*8.85e-12*m*g)

% function q = qfun(d,l,k)
%    q = k .* sqrt( d.^3 ./ sqrt( 4.0.*l.^2 - d.^2));
% end

% nemarané hodnoty vzdialenosti valčekov
d_vals = [0.065, 0.04, 0.02];      #  14.9.2022
q_vals = qfun(d_vals,l,k);

% formátovaný výpis pre namerané hodnoty 
for d = d_vals
    q = qfun(d,l,k);
    fprintf('d = %f m, q = %e C\n', d, q);
end

% graf
d_vals = linspace(0.0, 0.1, 20);
q_vals = qfun(d_vals,l,k);

plot(d_vals,q_vals);
xlabel("d (m)");
ylabel("q (C)");
grid on

% pause();
%##
%##  zadefunovanie používanej funkcie
%## 
function q = qfun(d,l,k)
    q = k .* sqrt( d.^3 ./ sqrt( 4.0.*l.^2 - d.^2));
end

end
