各位大哥，求助！

rand('state', sum(100*clock));
ns= 2^16;
np=400;
ss = 1;
p = zeros(np,2);
for a = 1:ns
for c = 1:np
theta = rand*2*pi;
p(c,1)= p(c,1)+ss*cos(theta);
if p(c,1) > 100
p(c,1)= p(c,1)-2*ss*cos(theta);
end
if p(c,1) < -100
p(c,1)=p(c,1)-2*ss*cos(theta);
end
p(c,2)= p(c,2)+ss*sin(theta);
if p(c,2) > 100
p(c,2)= p(c,2)-2*ss*sin(theta);
end
if p(c,2) < -100
p(c,2)=p(c,2)-2*ss*sin(theta);
end
end
end

题目——假设在一个空旷封闭的房间中心滴上一滴香水，挥发的香水分子随机的与空气中的粒子发生碰撞，最终会扩散到整个房间。这里我们将通过计算机模拟 400 个分子在二维平面内的随机运动情况，在 MATLAB 中定义一个 400×2 的矩阵来存放这些分子的位置。其中每一行代表一个分子，第一列存放位置的 x 分量，类似的第二列存放其 y 分量。