procedure kk(var v1,v2: vector; var u: ltuple);
var
  vt1,vt2: vector;
begin { kk }
  write('(kk'); flush(output);
  unpack(v1,vt1,u);
  if not uok(u,v1) then writeln('kk: error - u out of range');
  vscale(beta2,vt1,vt2);
  vsqr(tau,vt2,vt1);
  vconv(vdim,cov,vkappa,vt1,v2);
  write(')'); flush(output);
end { kk };

