function [distmatrix,xy] = gentspmat(npts)
% function [distmatrix,xy] = gentspmat(npts)

load usborder

X=[];
Y=[];
% Form the US border in imaginary coords for the INSIDE routine
w=xx+i*yy;
n=0;
while n<npts,
    a=rand*1.4+i*rand;
    if inside(a,w),
        X=[X; real(a)];
        Y=[Y; imag(a)];
        n=n+1;
    end;
end;
xy=[X Y];
distmatrix = zeros(npts);
for count1=1:npts,
    for count2=1:count1,
        x1 = xy(count1,1);
        y1 = xy(count1,2);
        x2 = xy(count2,1);
        y2 = xy(count2,2);
        distmatrix(count1,count2)=sqrt((x1-x2)^2+(y1-y2)^2);
        distmatrix(count2,count1)=distmatrix(count1,count2);
    end;
end;