
function [out] = kh_z1 (z)
zsol = 0.16;
lsol = 0.065;
S = 0.32;
N = 35;
k0 = 100;

for jj = 0:1:N
    pha1 = mystep1 (z - (zsol + jj*S - lsol/2));
    pha2 = mystep1 (zsol + jj*S +lsol/2 - z);
    HE (jj+1, :) = pha1 + pha2 -1;
end

HE = sum(HE);
out = k0*HE';
end