clear;

// *** 計算の設定 ***
r = 1;			// 1ステップの間に移動する距離
tnum = 10000;		// 時間ステップ数
t = (0:1:tnum)';	// 時間のベクトル

// *** 位置の計算 ***
theta = 2 * %pi * rand(1,tnum);
phi   = 2 * %pi * rand(1,tnum);
xi = r .* sin(theta) .* cos(phi);
yi = r .* sin(theta) .* sin(phi);
zi = r .* cos(phi);

S = [zeros(3,1), [cumsum(xi,'c'); cumsum(yi,'c'); cumsum(zi,'c')]];

// *** グラフのプロット ***
param3d(S(1,:),S(2,:),S(3,:));
xlabel("x position");
ylabel("y position");
zlabel("z position");