function [result] = closestPoints(x,y)



TRI = delaunay(x,y);


A = zeros(length(TRI),1);
ind = zeros(length(TRI),1);

for i = 1 :length(TRI)

    tempX = abs(x(TRI(i,1))-x(TRI(i,2)));
    tempY = abs(y(TRI(i,1))-y(TRI(i,2)));

    tempX1 = abs(x(TRI(i,1))-x(TRI(i,3)));
    tempY1 = abs(y(TRI(i,1))-y(TRI(i,3)));

    tempX2 = abs(x(TRI(i,2))-x(TRI(i,3)));
    tempY2 = abs(y(TRI(i,2))-y(TRI(i,3)));


    dist = norm([tempX,tempY]);
    dist1 = norm([tempX1,tempY1]);
    dist2 = norm([tempX2,tempY2]);

    [A(i),ind(i)] = min([dist,dist1,dist2]);



end

[~,ind1] = min(A);

switch ind(ind1)
    
    case 1
        
        result = [[x(TRI(ind1,1));x(TRI(ind1,2))],[y(TRI(ind1,1));y(TRI(ind1,2))]];
        
    case 2
        
        result = [[x(TRI(ind1,1));x(TRI(ind1,3))],[y(TRI(ind1,1));y(TRI(ind1,3))]];
    case 3
        
        result = [[x(TRI(ind1,2));x(TRI(ind1,3))],[y(TRI(ind1,2));y(TRI(ind1,3))]];
end
    

   plot(x,y,'.')
   hold on
   plot(result(:,1),result(:,2),'.r')
   
end