%TFY4104 Fysikk hoesten 2013.
%Oving 4, oppgave 1: Numerisk losning av bevegelsen
%til en matematisk pendel. Eulermetoden ("Forward Euler").
%
clear all;
format long;
%g = tyngdens akselerasjon (m/s^2)
g = 9.81;
%L = pendelens lengde (m)
L = 1.0;
%dt = tidssteget (s)
dt = 0.0001;
%t = tidsintervallet (s)
t = 0:dt:10;
%Definerer vinkel alpha(t) og hastighet v(t)
alpha = 0*t;
v = 0*t;
%Initialbetingelser for alpha (rad) og v (m/s)
alpha(1) = 0;
v(1) = 0.1;
%Losning med Euler-integrasjon.
for n = 1:(length(t)-1)
    alpha(n+1) = alpha(n) + v(n)*dt/L;
    v(n+1) = v(n) - g*dt*sin(alpha(n));
end;
%Plotter alpha(t)
grid on;
hold on;
plot(t,alpha*180/pi,'b');
xlabel('t (s)');
ylabel('vinkelutslag (grader)');
