%gaussiankernel.m
function [res]=gaussiankernel(x1,x2,sigma)
res=exp(-sum((x1-x2).^2)/(2*sigma^2));