% matOhjelmistot HT 19 (palautettava)

%% oletukset
h = 1.5 % metriä
g = 9.8 %m/s^2
v = 4 % m/s
theta = 45 %astetta
%% Ajan diskretointi
t = linspace(0,1,1000);

%% Korkeus ja etäisyys

x = v*cos(theta*pi/180).*t;
y = h + v*sin(theta*pi/180).*t - 0.5*g*t.^2;

%% Etsi putoamiskohta
%disp('Lentoaika')
thit = find(y<=0, 1 ) ;
disp('Lentomatka')
x(thit)
disp('Lentoaika')
t(thit)

%% Visualisointi
% Pallon rata
plot(x,y,'b',[0,max(x)],[0,0],'k--');