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