llin=load('distaztsts.dat')
a(:,1)=dcalc([0 0],llin);
a(:,2)=dcalc(llin,[0 0]);
c=a(:,1)-a(:,2)

a(:,3)=dcalc([0 180],llin);
a(:,4)=dcalc(llin,[0 180]);
c=a(:,3)-a(:,4)

llin=load('distazmoretsts.dat')
a(:,5)=dcalc(llin(end,:),llin(1:end-1,:));
a(:,6)=dcalc(llin(1:end-1,:),llin(end,:));
c=a(:,5)-a(:,6)

function a=dcalc(p1,p2)
[d a]=distance(p1,p2);
d
a
end
