MATLAB 蒙特卡洛

有多少个6阶0-1矩阵是非奇异矩阵？

蒙特卡洛算法估算的MATLAB代码：

y=0;
n=100000;
for i=1:n
    a=round(rand(6));
    if(det(a)~=0)
        y=y+1;
    end
end
answer=y/n*4096*4096*4096;