Points Inside polygone

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;