% Grafik zavisnosti perioda oscilovanja od pocetnog ugla 

clear all

l = 0.7;    % duzina kanapa (m)
g = 9.81;   % gravitacija m/s^2

theta0=0:0.1:pi/2;
for i=1:length(theta0)
    T1(i)=integral(@(x) 1./sqrt(cos(x)-cos(theta0(i))),0,theta0(i)) *4 *sqrt(l/(2*g)); % za proizvoljno theta
    T2(i)=2*pi*sqrt(l/g); % za malo theta <<1
end

plot(theta0,T1,theta0,T2)
xlabel('theta')
ylabel('period T')