clear;

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

// *** 位置の計算 ***
// 移動する距離
r = abs(grand(1,tnum,'nor',av,u));
// 移動する方向
theta = 2 * %pi * rand(1,tnum);
// 移動量
xi = r .* cos(theta);
yi = r .* sin(theta);
// 最終的な座標
S = [zeros(2,1), [cumsum(xi,'c'); cumsum(yi,'c')]];

// *** グラフのプロット ***
plot(S(1,:),S(2,:),'-b');
xlabel("x position");
ylabel("y position");