clear;
clc;

global m a b g r;

m = 1;
a = 0.15;
b = 0.01;
g = 9.81;
r = -0.5;

V0 = 10;
alfa0 = (pi/180)*45;

X0 = [0 0 V0*sin(alfa0) V0*cos(alfa0)]';
T = 0:1e-3:1.5;
[T1,Y1] = ode45('rzutfun', T, X0);
plot(T1,Y1);
pause;

% Wizualizacja
close all;
axis([-0.5 9 -1 9]);
set(gca,'box','on');
grid on;

% Wspolrzedne kuli
X = Y1(:,1);
Y = Y1(:,2);

% Rysujemy okrag, ktory bedziemy przesuwac
ii = 0:2*pi/20:(2*pi)+2*pi/80;
xx = 0.15*sin(ii);
yy = 0.15*cos(ii);
id = line(xx+X(1),yy+Y(1),'erasemode', 'xor');

for i=1:10:length(T1)
    % Przesuwamy kule
    set(id, 'xdata', xx+X(i),'ydata', yy + Y(i));

    % Wyswietlamy i czekamy
    drawnow;
    pause(0.025);
end;