Mercurial > hg > medcouple

changeset 41:2de6dcb8e9c2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
slowmedcouple: turn into a function
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Mon, 09 May 2016 19:15:49 -0400
parents 2493b21d7d8b
children 551327fe329b
files slowmedcouple.m
diffstat 1 files changed, 15 insertions(+), 14 deletions(-) [+]
line wrap: on
line diff
--- 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