十分感谢，但可能是我转化问题时不够确切，转化后的问题是可以通过排序的方法解决。但实际问题是a、b两个变量分别是不同的经纬度坐标，需要计算a中每一个坐标到b中每一个坐标的距离，然后取最近的三个距离之和，作为变量c的值，这就无法使用排序，然后sas的双重循环我不知道怎么表达。不过还是谢谢你回答我的问题。

如果你是坐标值，那么两个variable表示一个位置，你需要的是一个笛卡尔积，用sql实现，参见如下程序

data test;
input a1 a2  b1 b2;
cards;
3 1 7 6
4 2 3 7
1 3 4 8
2 4 6 9
6 5 5 5
;
data test1;set test;run;
proc sql;create table test2 as select a.a1,a.a2, b.b1,b.b2 from test as a cross join test1 as b;quit;
data test3;set test2;c0=sqrt((a1-b1)**2+ (a2-b2)**2); proc sort;by a1 a2  c0;
data test4;
  retain c;
  set test3;
  by a1 a2 c0;
  if first.a2 then do; ord=0;c=c0;end;
  else do; ord+1; c=c+c0; end;
  if ord=3;
  drop c0 ord;
run;