Points Inside polygone

Finds those points (xi,yi) who are inside the polygon, defined by the vectors (x,y)

Find those points (xi,yi) who are inside the polygone defined by the vectors (x,y).

Syntax: [xin,yin,Ind]=Inside(X,Y,xi,yi)

where Ind is the index of (xi,yi) that are inside the polygone

Example:

%--- create polygone (x,y)
phi=linspace(0,2*pi,50);
X=(1+0.2*sin(6*phi)).*cos(phi);
Y=(1+0.2*sin(6*phi)).*sin(phi);
%--- random points (xi,yi)
Nr=20;
xi=randn(1,Nr);yi=randn(1,Nr);
[xin,yin,Ind]=Inside(X,Y,xi,yi);
%--- plot
plot(X,Y,X,Y,'.',xi,yi,'ro',xin,yin,'b.');axis equal;