Triangulación y convexidad

Este es un programa en MATLAB, clickea en el axes y luego de varias veces haz clic derecho. Verás la triangulación por Delaunay y la capsula convexa. Si desean ver un poco de teoría sobre convexidad y cápsula convexa, revisen : Functional Analisys with Applications , Zeidler.

clear all
clc
figure(1)
xlim([0 1])
ylim([0 1])
hold on
i=0;
X=[];
Y=[];
while 1
    [x,y,buttom]=ginput(1);
    if buttom==3
        break
    end
    plot(x,y,'ok')
    i=i+1;
    X(i)=x;
    Y(i)=y;  
end

% si usas Octave solo cambia esto por lo último

%T=delaunay(X',Y');

%k = convhull(X',Y');

%triplot(T,X,Y)

%plot(X'(k),Y'(k), '-rs')

dt=DelaunayTri(X',Y');
k = convexHull(dt);
triplot(dt)
plot(dt.X(k,1),dt.X(k,2), 'rs')