suiji3.m

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

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