%% PallonHeitto.m, MIT problem set 1 probl. 11
% 
% Maaritellaan vakioita:
h=1.5;   % m       alkukorkeus
g=9.81;  % m/s^2
v=4;     % m/s     alkunopeus
Theta=45 % astetta lahtokulma

t=linspace(0,1,1000);

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

%% Pallon putoamisaika ja paikka:

               % help find
ind=find(y<0,1)% Etsit??n eka indeksi,jossa (y(ind)<0
xetaisyys=x(ind)
disp(['Pallo putoaa maahan etaisyydella ', num2str(xetaisyys),' m'])

%% Piirretaan pallon rata.
close all
plot(x,y)
xlabel('Matka')
ylabel('Korkeus')
title('Pallon lentorata')

hold on
plot([0 max(x)],[0 0],'k --')
grid on


