Untitled

for kk = 1:180 % određivanje domena (stepeni)
    teta(kk) = kk * pi / 180;
    TT = teta(kk);
    S(kk) = 0.5 * (1 - cos(TT)); % izraz za S
end

V = diff(S); % brzina
A = diff(V); % ubrzanje
J = diff(A); % udar


% Plot S, V, A, and J
figure; % Kreiranje novog prozora za dijagrame

subplot(4, 1, 1); % Kreiranje dijagrama za S
plot(teta,S,'k-','LineWidth',1.6); grid on;
xlabel('teta [rad]');
ylabel('S', 'FontWeight','bold');
title('Funkcija S(teta)');



subplot(4, 1, 2); % Kreiranje dijagrama za V
plot(teta(1:end-1),V,'b-','LineWidth',1.6); grid on; % Note: V has one less element than S
xlabel('teta [rad]');
ylabel('V','FontWeight','bold');
title({'';'Funkcija V(teta)'}); %{'';'Funkcija V(teta)'} - prazni apostrofi dodaju prazan red. 
                                %Naslov moze da bude u 2 reda
ax=gca;
ax.YAxis.Exponent=-2;

subplot(4, 1, 3); % Kreiranje dijagrama za A
plot(teta(1:end-2),A,'m-','LineWidth',1.6); grid on; % Note: A has two fewer elements than S
xlabel('teta [rad]');
ylabel('A','FontWeight','bold');
title('Funkcija A(teta)');



subplot(4, 1, 4); % Kreiranje dijagrama za J
plot(teta(1:end-3),J,'r-','LineWidth',1.6); grid on; % Note: J has three fewer elements than S
xlabel('teta [rad]');
ylabel('J','FontWeight','bold');
title('Funkcija J(teta)');


% Adjust subplot spacing
%set(gcf, 'subplotspacing', 1);