Random walk-kropp

Vi lager en kropp basert på stien etter en "agent" som beveger seg tiflfeldig i vG
NB: Har ingen test her på om agenten kommer uttafor vG, så her blir det lett feilmeldinger og man må nok kjøre flere ganger for å treffe (eller legge inn test for uttafor)
vG = zeros(50,50,50,'int8');
agentPos = int16([50/2 50/2 50/2]); % Utgangsposisjon, lar den være av type int siden vi snakker om indexposisjon i vG
for n = 1 : 200
vG(agentPos(1), agentPos(2), agentPos(3)) = 1; % Setter inn en voksel der agenten er
agentPos = agentPos + 2*int16(rand(1,3) > 0.5)-1; % Trekker ny random posisjon, sentrerer rundt 0, og caster til int
end
figure(1);clf;plotVg_safe(vG);
Legger på tykkelse i form av en 3x3 kube med voksler rundt agenten
vG = zeros(50,50,50,'int8');
agentPos = int16([50/2 50/2 50/2]);
for n = 1 : 200
vG(agentPos(1)-1:agentPos(1)+1, agentPos(2)-1:agentPos(2)+1, agentPos(3)-1:agentPos(3)+1) = 1; % Tykkelse med range-indexering
agentPos = agentPos + 2*int16(rand(1,3) > 0.5)-1;
end
figure(2);clf;plotVg_safe(vG);
Random kule-walk
radius = 20;
vG = zeros(250,250,250,'int8');
agentPos = int16([250/2 250/2 250/2]);
for n = 1 : 450
vG = vGsphere(vG,agentPos,1,radius);
agentPos = agentPos + 8*int16(rand(1,3) > 0.5)-4;
end
figure(3);clf;plotVg_safe(vG,'edgeOff');