function [K] = gausskernel(Xt,X,param);

[n,m] = size(X);
[nt,m] = size(Xt);

K = 2*ones(nt,n)-2*Xt*X'; 
K = exp(-K/2/param);


