% FY1001/TFY4145 Mekanisk fysikk hosten 2015.
% Oving 7, oppg. 2: Matematisk pendel

% Definere parametre
% Velg tilfelle = 'A','B', eller 'C'
tilfelle='B';
% For tilfelle A:
dt=1e-2;	% Størrelsen på tidssteg
T=2.0;		% Periode
v0=1.0;		% Starthastighet

% For tilfelle B:
if (tilfelle=='B')
dt=1e-2;	% Størrelsen på tidssteg
T=8.7531;	% Periode
v0=6.2641;	% Starthastighet
end

% For tilfelle C:
if (tilfelle=='C')
dt=1e-2;	% Størrelsen på tidssteg
T=1.684;	% Periode
v0=6.5;		% Starthastighet
end

theta0=0.0;	% Startvinkel
L=1.0;		% Pendelens lengde
N=floor(2.0*T/dt);% Antall tidssteg som tilsvarer to perioder

% Kalle funksjonen mat_pendel som er definert i filen mat_pendel.m
mat_pendel(L,v0,theta0,dt,N,tilfelle);
