%TFY4106 Fysikk 2014.
%Øving 7, oppgave 2: Numerisk løsning av bevegelsen
%til en matematisk pendel. Eulermetoden ("Forward Euler").
%
clear all;
format long;
%Tilordne verdier for g, L og tidssteget dt:

%Lag "tidstabell" t med verdier mellom 0 og 2*perioden, med tidssteg dt.
%(For (A) harmonisk svingning kan perioden beregnes. For tilfellene (B)
%og (C), prøv deg fram.

%Nullstill tabeller for vinkel theta(t) og hastighet v(t):
theta = 0*t;
v = 0*t;
%Initialbetingelser for theta (rad) og v (m/s):
theta(1) = 0.0;
v(1) = ...;
%Beregn startverdi for mekanisk energi:
E(1) = ...;
%Løsning med Euler-integrasjon:
for n = 1:(length(t)-1)
    theta(n+1) = theta(n) + v(n)*dt/L;
    v(n+1) = v(n) - g*dt*sin(theta(n));
    %Mekanisk energi:
    E(n+1) = ...;
end;
%Plott theta(t); 'b' angir blå kurvefarge.
%Enheter: Grader for theta, sekunder for t.
figure(1);
plot(t,theta*180/pi,'b');
%Sett tekst på aksene med xlabel og ylabel:

%Lag tittel med title:

%I menyen til den aktuelle figuren, velg File - Save as og filtype pdf

%Plott E(t) i enheten J/kg:
figure(2);
plot(...);
%Sett tekst på aksene med xlabel og ylabel:

%Lag tittel med title:

