matlab如何做粒子模拟,求助，如何用matlab做蒙特卡罗模拟！！？？

我这儿正好有份程序，希望有所帮助。

代码：

一份蒙特卡洛程序

count=input('input the count:');%输入模拟粒子数

sigmaedata=[28370,13845,6908,2555,1223,2602,1925,905.3,479.7,164.5,74.24,23.86,66.60,36.62,22.29,9.978,5.298,1.668,0.7378,0.2361,0.1099,0.06211,0.03939,0.02030];

sigmacdata=[0.0220,0.0393,0.0568,0.0904,0.1209,0.1479,0.1722,0.2136,0.2480,0.3092,0.3486,0.3932,0.4153,0.4268,0.4319,0.4291,0.4215,0.3969,0.3691,0.3269,0.2944,0.2709,0.2512,0.2209];

Edata=[1,1.5,2,3,4,5,6,8,10,15,20,30,40,50,60,80,100,150,200,300,400,500,600,800];%截面数据

channel=zeros(1,ceil(662/5)+10);%多道数组

nget=0;%探测到的总计数

ntotal=0;%进入探测器的总计数

for ii=1:count%count个粒子循环

collidetime=0;%当前粒子碰撞次数

%粒子状态初始化

E0=622;

E=E0;

z=-2;

r=0;

theta=2*pi*rand(1);% 源抽样，z,r,theta坐标

miu=2*rand(