%pallonheitto

%% alkuehdot
h0 = 1.5;
g = 9.81;
v0 = 4;
theta = pi/4;

%% ajan diskretointi

t = linspace(0,1,1000);

%% paikan approksimointi

x = v0*cos(theta)*t;
y = h0 + v0*sin(theta)*t - 0.5*g*t.^2;
%% arvioidaan maahanosumishetkeä.
ind = find(y<0,1,'first');
lentomatka = x(ind)
lentoaika = t(ind)

%% piirretään tulokset

plot(x,y);
hold on 
plot([0 x(end)], [0 0],'k--')
plot(x(ind),y(ind),'ro')
xlabel('matka')
ylabel('korkeus')
hold off