# HG changeset patch
# User Jordi Gutiérrez Hermoso <jordigh@octave.org>
# Date 1421418729 18000
# Node ID 3958c6e774429b6c0699a8aebb2fe880e2ed30c2
# Parent  549846d50ddb0f92ef0c198747d7d3caf3965260
Add slow Octave medcouple implementation

diff --git a/slowmedcouple.m b/slowmedcouple.m
new file mode 100644
--- /dev/null
+++ b/slowmedcouple.m
@@ -0,0 +1,14 @@
+z = load("gamma2");
+n = length(z);
+n2 = ceil(n/2);
+z = sort(z, "descend");
+zmed = median(z);
+z -= zmed;
+zplus = z(z >= 0);
+zminus = z(z <= 0)';
+zz = (zplus + zminus)./(zplus - zminus);
+zz(isnan(zz)) = 0;
+zzsort = sort(zz(:));
+medc_idx = ceil(n2^2/2);
+mc = zzsort(medc_idx);
+printf("%.16g\n", mc);