load ../test_data/honshu.mat;
figure;
plot(honshu(:,1),honshu(:,2)) ;
hold on;
load ../test_data/seto_grid.mat;
inds = mexinside(xgrid,ygrid,honshu(:,1),honshu(:,2));
ind2 = find(inds);
plot(xgrid(ind2),ygrid(ind2),'r*');
ind3 = find(1-inds);
plot(xgrid(ind3),ygrid(ind3),'g*');
plot(honshu(:,1),honshu(:,2));      


titlestr = sprintf ( 'Test of Mexinside:  Southern coast of Honshu, Japan' );
title ( titlestr );

text ( 3.5e5, 3.65e6, 0, 'Red dots should be inside coast (north)' );
text ( 3.5e5, 3.62e6, 0, 'Green dots should be outside coast (south)' );



