function res = jk(t,cj,k,hk)

s = 0;

for j = 1:k
    s = s + cj(j,k)*hk(j,t);
end

res = s;




end

