function [res]=k33(x1,x2,c4,c5)
x1=x1/sqrt(sum(x1.^2));
x2=x2/sqrt(sum(x2.^2));
res=exp(-sum((x1-x2).^2)/(c4^2))^(c5);