randomMatrixDom[rnm_]:=Table[RandomReal[{-1,1}]&[i],{i,rnm},{j,rnm}]+(rnm+1)*IdentityMatrix[rnm];
test[mat_]:=Module[{absmat,absmatBezDiagonaly,diagonala,soucty},
absmat=mat/.a_:>Abs[a];
diagonala=Diagonal[absmat];
absmatBezDiagonaly=absmat-DiagonalMatrix[diagonala];
soucty=Total/@absmatBezDiagonaly;
Max[(soucty/diagonala)]<1

]
test[randomMatrixDom[3]]

Timing[Integrate[((Sin[x])^(10^3)),x]][[1]]

Range[3]

Thread[{a,b,c}=={d,e,f}]

randomMatrixDom[5].((Subscript[x, #])&/@Range[5])

dejRovnice[rnm_]:=Thread[randomMatrixDom[rnm].((Subscript[x, #])&/@Range[rnm])==Range[rnm]];
Timing[Solve[dejRovnice[50]]][[1]]