harj33.m

% Lentävän esineen lentorata
xcoord = @(t) tan(t/10);
ycoord = @(t) exp(-t/3).*sin(t)-t.^2/20+3;

% Muodosta x ja y koordinaatit käyttäen funktioita
t = 0:0.3:10;
x = xcoord(t);
y = ycoord(t);

figure(1);
hold off
plot(x,y,'-+b')
xlabel('x');
ylabel('y');
hold on
plot([0 1.6],[0 0],'k-');

figure(2);
hold off
plot(t,y,'-+b')
xlabel('t');
ylabel('y');
hold on
plot([0 10],[0 0],'k-');

for i=1:length(t)
  if y(i) < 0
     break;
  end
end
i = i+1;
t0 = t(i) + (t(i+1)-t(i)) * (0-y(i)) / (y(i+1)-y(i))
x = x(i) + (x(i+1)-x(i)) * (t0-t(i)) / (t(i+1)-t(i))