clear all
clc

% Vypocet pohybu hmotneho bodu v gravitacnom poli zeme
% s pociatocnou rychlostou bez uvazovania aodporu vzduchu

m = 1;          % [kg] hmotnost kamena
g = 9.80665;    % [m/s^2] granitacne zrychlenie
v0 = 6;         % [m/s] velkost pociatocnej rychlosti
theta = 60*pi/180;  % % [rad] uhol azumutu

v0x = v0*cos(theta);
v0y = v0*sin(theta);

t = 0:0.01:2;   % casovy vektor

% pohyb v smere osy X
vx = v0x;
x = v0x*t;
% pohyb v smere osy Y
vy = -g*t + v0y;
y = -0.5*g*t.^2 + v0y*t;

% RYCHLOSTI
figure(1)
subplot(2,1,1)
plot(t,vx)
xlabel('t [s]'), ylabel('v_x [m/s]'),title('Rychlost hmotneho bodu v smere X')
grid on
subplot(2,1,2)
plot(t,vy)
xlabel('t [s]'), ylabel('v_y [m/s]'),title('Rychlost hmotneho bodu v smere Y')
grid on

% POLOHY
figure(2)
subplot(2,1,1)
plot(t,x)
xlabel('t [s]'), ylabel('x [m]'),title('Poloha hmotneho bodu v smere X')
grid on
subplot(2,1,2)
plot(t,y)
xlabel('t [s]'), ylabel('y [m]'),title('Poloha hmotneho bodu v smere Y')
grid on

% TRAJEKTORIA
figure(3)
plot(x,y)
xlabel('x [m]'), ylabel('y [m]'), title('Trajektoria')
grid on