# HG changeset patch
# User Jordi Gutiérrez Hermoso <jordigh@octave.org>
# Date 1462835749 14400
# Node ID 2de6dcb8e9c2e6d5b4000a63c6eee18362a94d31
# Parent  2493b21d7d8b498bf1ebbfeba3e9a90cbe5c0dc5
slowmedcouple: turn into a function

diff --git a/slowmedcouple.m b/slowmedcouple.m
--- a/slowmedcouple.m
+++ b/slowmedcouple.m
@@ -1,14 +1,15 @@
-z = load("rand53");
-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(:), "descend");
-medc_idx = ceil(n2^2/2);
-mc = zzsort(medc_idx);
-printf("%.16g\n", mc);
+function mc = slowmedcouple(z)
+  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(:), "descend");
+  medc_idx = ceil(n2^2/2);
+  mc = zzsort(medc_idx);
+  printf("%.16g\n", mc);
+endfunction