@ -43,6 +43,8 @@ import Base.LinAlg.det
function det(M::Array{Mod,2})
if size(M,1) ≠ size(M,2)
d = Mod(0,M[1,1].mod)
elseif size(M,1) == 1
d = M[1,1]
elseif size(M,1) == 2
d = M[1,1]*M[2,2] - M[1,2]*M[2,1]
else