function mat_pendel(L,v0,theta0,dt,N,tilfelle)
g=9.81;
dtL=dt/L;
% Initialisere tabeller for hastighet, vinkel, energi og tid
v=zeros(N);
v(1)=v0;
theta=zeros(N);
theta(1)=theta0;
E=zeros(N);
E(1)=0.5*v(1)^2+2*g*L*(1-cos(theta(1)));
tid=zeros(N);
tid(1)=0;
% Starter iterasjonen i 1'te tidssteg
t=1;
% Øvre grense for iterasjonen
t_max=N;
while (t < t_max)
	v(t+1)=v(t)-g*sin(theta(t))*dt;
	theta(t+1)=theta(t)+v(t)*dtL;
	E(t+1)=0.5*v(t+1)^2+g*L*(1-cos(theta(t+1)));
	tid(t+1)=t*dt;
	t=t+1;
end
% vinkler større enn 2*pi flyttes til intervallet (-2*pi,2*pi)
theta=mod(theta,2*pi);
% Plotte hastigheten
figure;
plot(tid,v);
% Lage tittel
title(sprintf('Hastighet, tilfelle %s med v_0=%.4f',tilfelle,v0));
% Merke aksene
xlabel('t (s)');
ylabel('v (m/s)');
% Plotte vinkelen
figure;
plot(tid,theta*180/pi);
title(sprintf('Vinkel, tilfelle %s med v_0=%.4f',tilfelle,v0));
xlabel('t (s)');
ylabel('\theta (grader)');
% Plotte energien
figure;
plot(tid,E);
title(sprintf('Energi, tilfelle %s med v_0=%.4f',tilfelle,v0));
xlabel('t (s)');
ylabel('E/m (m^2/s^2)');
