function y=Kg(gamma,x1,x2)
    y = exp(-gamma*(x1'*x1 + x2'*x2 - 2*x1'*x2));