liti10.m

来自「此文件包括蒙特卡罗方法的基本思想及应用」· M 代码 · 共 21 行

mm=10000
Randnum=unifrnd(0,1,1,mm);
xRandnum=zeros(1,mm);
for ii=1:mm
    if Randnum(1,ii)<=0.1
       xRandnum(1,ii)=10;
       else if  Randnum(1,ii)<=0.3
            xRandnum(1,ii)=20;
                else if  Randnum(1,ii)<=0.7
                         xRandnum(1,ii)=30;
                         else  if  Randnum(1,ii)<=0.95   
                               xRandnum(1,ii)=40;
                               else                          
                               xRandnum(1,ii)=50;
                             end
                     end
            end
      end     
 end      
xRandnum
cdfplot(xRandnum)