suiji4.m

来自「有兴趣学系matlab的朋友」· M 代码 · 共 16 行

n=10000;
a=[];
for i=1:n
    m=floor(rand(1,6)*9+1);
    for l=1:6
        if floor(rand*2+1)==1
            m(l)=1/m(l);
        end
    end
    m=log(m);
    p=[(m(1)+m(4)-m(2))^2 (m(1)+m(5)-m(3))^2 (m(2)+m(6)-m(3))^2];
    s=sum(p);
    a=[a s];
end
b=sum(a)/n