thresh = 50.0;
probs = 0.5:0.01:1.0;
reps = 1000;
meancount = zeros(length(probs), 1);
stdcount = zeros(length(probs), 1);
for sp = 1:length(probs)
	p = probs(sp);
	positions = zeros(reps, 1);
	for k = 1:reps
		x = randomwalkthresh(p, thresh);
		positions(k) = length(x);
	end
	meancount(sp) = mean(positions);
	stdcount(sp) = std(positions);
end
semilogy(probs, meancount, 'displayname', 'mean');
hold on;
semilogy(probs, stdcount, 'displayname', 'std');
hold off;
xlabel('p');
xlim([0.5 1.0])
ylabel('number of steps');
legend('show');