# HG changeset patch
# User Jordi Gutiérrez Hermoso <jordigh@octave.org>
# Date 1421542386 18000
# Node ID 6339e1a06aa9f0a2bc7e91a5b6a46ca0226067bd
# Parent  8e35dcdb8dec4f0822eda671b51385573719be02
slowmedcouple: sort in reverse

diff --git a/slowmedcouple.m b/slowmedcouple.m
--- a/slowmedcouple.m
+++ b/slowmedcouple.m
@@ -8,7 +8,7 @@
 zminus = z(z <= 0)';
 zz = (zplus + zminus)./(zplus - zminus);
 zz(isnan(zz)) = 0;
-zzsort = sort(zz(:));
+zzsort = sort(zz(:), "descend");
 medc_idx = ceil(n2^2/2);
 mc = zzsort(medc_idx);
 printf("%.16g\n", mc);