%FY1001/TFY4145 Mekanisk fysikk hoesten 2014.
%Oving 1, oppgave 3c: Plotting av funksjonen
%sin(theta)*cos(theta), som angir lengden av et
%skraatt kast som funksjon av utgangsvinkelen theta.
%Semikolon til slutt paa linja forhindrer utskrift til
%kommandovinduet.
clear all;
%N = antall punkter langs theta-aksen
N = 100;
%tabell med N theta-verdier mellom theta_min og theta_max
theta_min = 0;    %tilsvarer horisontalt kast
theta_max = pi/2; %tilsvarer vertikalt kast
theta = linspace(theta_min,theta_max,N);
%deretter lager vi en tabell med verdier for lengden L,
%eller rettere sagt, den dimensjonsloese stoerrelsen
%L*g/v_0^2;
%legg merke til at det skrives .* og ikke bare * for
%multiplikasjon. Punktumet foran * resulterer i at
%L blir en tabell med N verdier
L = sin(theta).*cos(theta);
%Et alternativ til denne ene linja er aa bruke en for-loekke:
%for i = 1:N
%    L(i) = sin(theta(i))*cos(theta(i));
%end;
%Med for ... end involveres kun ett element
%i theta-tabellen hver gang, og da er ikke
%punktumet foran * noedvendig. (Men det hadde
%heller ikke gjort noen skade :-))
%Vi plotter L(theta) med linjetykkelse = 2
plot(theta,L,'linewidth',2);
% Aksetekster:
xlabel('Utgangsvinkel (rad)','fontsize',14);
ylabel('Lg/v_0^2 (dimensjonsløs lengde)','fontsize',14);
% Tittel:
title('Øving 1, oppgave 3c','fontsize',14);
