clear all;
close all;
clc;

% Zadej hodnoty
alpha = input('Zadej úhel: ');
v0 = input('Zadej počáteční rychlost: ');
%Krok času
dt = 0.1;

%Počáteční podmínky
vx(1) = v0*cos(deg2rad(alpha));
vy(1) = v0*sin(deg2rad(alpha));

px(1) = 0;
py(1) = 0;

%Vypočet šikmého vrhu
t = 0;
I = 1;

while py(I)>=0
    t = t + dt;
    I = I + 1;

    ax(I) = 0;
    ay(I) = -9.81;

    vx(I) = vx(I-1) + ax(I)*dt;
    vy(I) = vy(I-1) + ay(I)*dt; 

    px(I) = px(I-1) + vx(I)*dt;
    py(I) = py(I-1) + vy(I)*dt;
    
    plot(px,py);
    axis([0, 100, 0, 50]);
    pause(dt);
end