Add files via upload

private/myhistc.m

@@ -0,0 +1,20 @@
+function col = myhistc(x, t)
+% col = myhistc(x, t)
+% as [~, col] = histc(x,t) but with specific col value for x == max(t)
+% that does not take the last index, meaning the last bin contains also the
+% right bound
+
+[~, col] = histc(x,t); %#ok
+n = length(t);
+b = col == n;
+if any(b)
+    tmax = t(end);
+    for colmax = n:-1:1
+        if t(colmax) ~= tmax
+            break
+        end
+    end
+    col(b) = colmax;
+end
+
+end % myhistc

private/mynormest.m

@@ -0,0 +1,10 @@
+function nA = mynormest(A)
+% nA = mynormest(A)
+% Estimate 2-norm of matrix A
+
+if issparse(A)
+    nA = normest(A);
+else
+    nA = norm(A);
+end
+end