clear; // *** 計算の設定 *** r = 1; // 1ステップの間に移動する距離 tnum = 10000; // 時間ステップ数 t = (0:1:tnum)'; // 時間のベクトル // *** 位置の計算 *** 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");