# -------- Initialize ---------- isochrones = rand (73160, 29); lw(1, 1, 1, :) = isochrones (:, 29); j(1, 1, 1, :) = isochrones(:, 15); k(1, 1, 1, :) = isochrones(:, 17); ukidss = rand (2, 77); oj(1, 1, :, 1) = ukidss(:, 40); ok(1, 1, :, 1) = ukidss(:, 76); sj(1, 1, :, 1) = ukidss(:, 41); sk(1, 1, :, 1) = ukidss(:, 77); e = 0.1/2:0.1:19.75; KJ = 0.2833; KK = 0.11424; d = (0:0.1:12)'; m = convn (d, [1; 1]/2, 'valid'); md = 5*log10 (m*1000) - 5; # -------- The processing ---------- x = parcellfun (2, @exp_arg, {j, k}, {md, md}, {e, e}, {KJ, KK}, {sj, sk}, {oj, ok}, "UniformOutput", false);