Охват точек линией минимальной длинны в MATLAB

clear, clc
% случайные точки
N = 10;
x = randn(N,1);
y = randn(N,1);
plot(x,y,'.') % рисуем исходные точки
hold on 
p = x + 1i*y; % переводим координаты в комплексный вид, чтобы было проще записывать
% берем самую левую точку - это первая точка оболочки z
z(1) = p( x == min(x) ); 
k = 1; % счетчик точек
% алгоритм Джарвиса
while z(end) ~= z(1) | k == 1 % пока не вернемся к начальной точке (или только первый шаг)
    % находим разности всех точек и текущей: p-z(k) 
    % вычисляем аргумент этих разностей (angle)
    % функция unwrap убирает разрывы:
    f1 = unwrap( angle( p-z(k) ) ) ;
    % игнорируем точки, у которых аргумент получился +/-пи, 
    % то есть они находятся позади, чтобы не возвращалось:
    f1( abs(f1)-pi==0 ) = NaN; 
    % сохраняем в массив z точку с минимальным аргументом:
    z(k+1) = p( f1==min(f1) );
    p( p==z(k+1) ) = []; % убираем из исходного набора эту точку
    k = k + 1; % увеличиваем счетчик
end
plot(z,'d-r') % рисуем оболочку
hold off
% периметр - это сумма (sum) отрезков оболочки,
% длинна каждого равна модулю (abs) разности соседних (diff) точек
P = sum( abs( diff(z) ) );
title(['Периметр P = ',num2str(P)]) % подписываем