%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% randomwalk.m:
% Ekstraoppgave 3 FY1005/TFY4165 Termisk fysikk 2014:
% Simulering av virrevandring på 2D kvadratisk gitter
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Nullstiller variable og lukker figurer
clear all;
close all;
% Antall skritt
N = 1e6;
% Starter i origo
x = 0;
y = 0;
% Tabell for vandret vei
x1 = zeros(1,N);
y1 = zeros(1,N);
for i = 2:N
    ds = randi(4);
    if ds == 1
       x = x+1;
       x1(i) = x1(i-1)+1;
       y1(i) = y1(i-1);
    elseif ds == 2
       y = y+1;
       x1(i) = x1(i-1);
       y1(i) = y1(i-1)+1;
    elseif ds == 3
       x = x-1;
       x1(i) = x1(i-1)-1;
       y1(i) = y1(i-1);
    else
       y = y-1;
       x1(i) = x1(i-1);
       y1(i) = y1(i-1)-1;
    end;
end
% Regn ut vandret lengde
r = sqrt(x^2 + y^2)
% Plotter vandret vei
plot(x1,y1);
text(0,0,'START');
text(x1(N),y1(N),'END');