clear;

// *** 計算の設定 ***
u = 1;					// 1ステップの間の平均二乗変位平方根
av = 0;
tnum = 20;				// 時間ステップ数
t = (0:1:tnum)';			// 時間のベクトル

// *** 位置の計算 ***
// 移動距離が正規分布に従う場合
xigauss = grand(tnum,1,'nor',av,u);
Sgauss = [0; cumsum(xigauss)];
// 移動距離が一定の場合
xifix = 2 * u * (rand(tnum,1) >= 0.5) - u;
Sfix =  [0; cumsum(xifix)];

// *** グラフのプロット ***
plot(t,Sgauss,'-or');
plot(t,Sfix,'-ob');
xlabel("x position");
ylabel("y position");