%prob1

T_1 = [1, 0.1, 0.01, 0.001];
n1=size(T_1);

for j=1:n1;
	alpha_1 = 1.0;
	T = T_1(j);
	s_1 = sqrt(T)*alpha_1;
	n=10/T;
	t_1 = (o:s_1:n*s_1);
	
	
	random_u = random('uniform',0,1,1,n+1);
	score = zeros(1,n+1);
	
	for k=1:n+1;
		if random_u(k) <=0.5;
			score(k) = -1;
		else
			score(k) = 1;
		end;
	end;


sum_1 = score(1,n+1);

for m=2:n+1;
	sum_1(m) = sum(score(1:m))*s;
end;

figure(j);
stars(t_1,sum_1);
xlabel('time')
ylabel('x_t')
title(['Total 'int2str(n)' Step random walk']);
saveas(figure(j),['Total 'int2str(n)' Step random walk.png'])
end;