H.J. Sommer

Excellent function, a real time saver

Raghuram, the correct answer is given if you reorder your vertices so they are given clockwise around the outside of the polygon, rather than just specified randomly, e.g,

% Your vertices
% x = [1.0000 0.5000 0.8333 0.5694];
% y = [1.0000 0.1667 0.5000 0.5694] ;

% reordered moving clockwise around the polygon from (0.5, 0.1667)

xy = [0.5, 0.1667;
      0.8333, 0.5;
      1, 1;
      0.5694, 0.5694; ];
      
x = xy(:,1);
y = xy(:,2);

for i = 1:numel(x)-1
    line(x(i:i+1), y(i:i+1), [0,0], 'Color', 'b');
    hold on
end

line([x(end), x(1)], [y(end), y(1)], [0,0], 'Color', 'b');

[ geom, iner, cpmo ] = polygeom( x, y )

plot(geom(2), geom(3), '+r')

hold off

This prog is wonderful, just what I needed. If you have time please improve it to handle the third dimension!!