% Program name: monte_diffusion.m
%
for k_trace = 0: 2
N = 200
imax = 100;
jmax = 100;
nframes = 20; % number of frames for movie
i = 50;
j = 50;
rand('state', k_trace) % this sets the generator to an initial state
for k = 1: 4000
% left = 0; right = 0;
% top = 0; bottom = 0;
for kk = 1: 2
rr = rand;
% if((rr >= 0 ) & ( rr < 0.25) & (right == 1))
if((rr >= 0 ) & ( rr < 0.25) )
i = i + 1;
j = j;
elseif((rr >= 0.25 ) & ( rr < 0.50) )
i = i;
j = j + 1;
elseif((rr >= 0.50 ) & ( rr < 0.75))
i = i - 1;
j = j;
elseif((rr >= 0.75 ) & ( rr <= 1.00))
i = i;
j = j - 1;
end
end
if i >= imax
break
end
if i <= 1
break
end
if j >= jmax
break
end
if j <= 1
break
end
axis([0 imax 0 jmax]);
if k_trace == 0
plot(i, j, 'r:.');
elseif k_trace == 1
plot(i, j, 'b:.');
elseif k_trace == 2
docsity.com
