function [res]=k55(x1,x2,c8)
x1=x1/sqrt(sum(x1.^2));
x2=x2/sqrt(sum(x2.^2));
res=1+(sum((x1-x2).^2)/c8);
res=1/res;
