function factor_ab = Interpolate2Lfactor(a,b);
factor_ab = 1/ (( 1+ (a-b)^2 )^.5 );