function result = Interpolate2L(a, b, r);
factor_ar = 1/( 1+ (a-r)^2 )^.5;
factor_rb = 1/( 1+ (r-b)^2 )^.5;
result = (factor_ar * a + factor_rb * b ) / (factor_ar + factor_rb);

