蒙特卡洛模拟法

function [out1 out2 out3 out4]=Msc（N）
%N开炮次数
%out1射中概率
%out2平均每次击中的火炮个数
%out3击中敌人一门火炮的概率
%out4击中敌人两门火炮的概率
%K1是击中零门火炮的次数
%K2是击中一门火炮的次数
%K3是击中两门火炮的次数
N=50;
k1=0;
k2=0;
k3=0;
out1=（k2+k3）/N;
out2=（0*k1+k2+2*k3）/20;
out3=k2/N;
out4=k3/N;
for i=1:50
    x0=randperm（2）-1;%随机产生1和2两个数
    y0=x0（1）;
    if y0==1
        fprintf（‘第%d次：指示正确||‘i）;
        x1=randperm（6）;