We combute hi th number to elemnts in the i-th equidistant subdivision of the interval [min L, max L] into n parts
i1 : M=(randomChainComplex({20,20},{20},ZeroMean=>true)).dd_1;
40 40
o1 : Matrix ZZ <--- ZZ
|
i2 : (svds,U,Vt)=SVD(M**RR_53); |
i3 : (entries matrix {svds})_0/log
o3 = {6.37106, 6.31472, 6.27245, 6.10348, 6.02102, 5.98252, 5.92934, 5.83927, 5.72509, 5.63923, 5.51957, 5.51441, 5.45378,
----------------------------------------------------------------------------------------------------------------------------
5.3237, 5.14787, 5.11063, 4.80679, 4.71988, 4.56427, 3.9834, -30.0523, -30.1114, -30.197, -30.4056, -30.4379, -30.4944,
----------------------------------------------------------------------------------------------------------------------------
-30.622, -30.7326, -30.7829, -30.8541, -31.0555, -31.2725, -31.5534, -31.7353, -31.8644, -31.9208, -32.1599, -32.8495,
----------------------------------------------------------------------------------------------------------------------------
-33.7088, -35.1847}
o3 : List
|
i4 : maximalEntry M o4 = 138 o4 : RR (of precision 53) |
i5 : histogram(svds/log,10)
o5 = {10, 10, 0, 0, 0, 0, 0, 0, 0, 20}
o5 : List
|
i6 : histogram(svds_{0..19}/log,10)
o6 = {1, 0, 1, 2, 2, 1, 4, 2, 4, 3}
o6 : List
|
i7 : histogram(svds_{20..39}/log,10)
o7 = {1, 0, 1, 0, 1, 1, 3, 2, 5, 6}
o7 : List
|