%FY1001/TFY4145 Mekanisk fysikk hoesten 2012.
%Oving 2, oppgave 4c: Plotting av 
%bane z(x) for en serie med utgangsvinkler.
clear all;
%N = antall punkter langs x-aksen
N = 100;
%g = tyngdens akselerasjon (m/s^2)
g = 9.81;
%L = horisontal avstand mellom start- og landingsrampe (m)
L = 8.0;
%h = vertikal forskjell mellom kanten av start- og landingsrampe (m)
h = 4.0;
%v0 = hastighet (m/s)
v0 = 32/3.6;
%vi lager tabell med N x-verdier mellom x_min=0 og x_max=L
x_min = 0;
x_max = L;
x = linspace(x_min,x_max,N);
%Med kommandoen hold on; vil gjentatte plot-kommandoer tegne
%kurvene i samme figur. Det passer fint her. Fra mathworks.com:
%"hold: Retain current graph when adding new graphs"
%Fjern figuren fra skjermen hvis du vil kjoere programmet paa
%nytt med nye tallverdier.
hold on;
%Vi skal lage og plotte banen z(x) for ekvidistante vinkler theta.
for i = 1:4
   theta = i*20*pi/180;
   %lager tabell med z(x) for aktuell hastighet v0
   z = x*tan(theta)-g*x.*x/(2*v0*v0*cos(theta)*cos(theta));
   %NB: Merk x.*x (og ikke x*x) i linja over!
   %Vi plotter z(x) med linjetykkelse = 1.0
   plot(x,z,'linewidth',1.0);
end;
%axis([xmin xmax zmin zmax]) gir plottets utstrekning
axis([0 L 0 h]);
% Aksetekster:
xlabel('x (m)','fontsize',14);
ylabel('z (m)','fontsize',14);
% Tittel:
%title('Øving 2, oppgave 2: Sykkelhopp','fontsize',14);