# -------- 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);
