t = 0:0.5:6.28;
seedx = 100;
seedy = 100;
scalefactor = 20;
P(:,1) = (seedx(1) + scalefactor*cos(t))';
P(:,2) = (seedy(1) + scalefactor*sin(t))';
plot(P(:,1),P(:,2));

%for i=1:3
    %[pointx, pointy] = ginput(1);    
    %q = [pointx pointy];
    q = P(13,:)
    str = BOPointInPolygon(P,q)
%end