clear;

// *** 時間ステップ数 ***
tnum = 10000;		// 時間ステップ数
t = (0:1:tnum)';	// 時間のベクトル

// *** 位置の計算 ***
xi = 2 * (rand(tnum,1) >= 0.5) - 1;
S = [0;cumsum(xi)];

// *** グラフのプロット ***
plot(t,S,'-b');
xlabel("Time");
ylabel("Position");