suiji5.m

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

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