function k=univker(x,y)
global KERD

k = 0;
for ai = 1:length(KERD{3})
  bx = run_dfsa(KERD{3}(ai),x);
  by = run_dfsa(KERD{3}(ai),y);
  k = k + bx*by;
end

return