%% Teht. 19, Monte Carlo
close all
N=10;
x=rand(1,N);
y=rand(1,N);

nelio=[0 0;1 0;1 1;0 1;0 0]
plot(nelio(:,1),nelio(:,2))
kolmio=[0 0;1 0;1 1;0 0]
hold on
plot(kolmio(:,1),kolmio(:,2))
axis([-.2 1.2 -.2 1.2 ])

plot(x,y,'O')

Isisa=(y<=x)
Iulko=~Isisa
%%
xsisa=x(Isisa)
ysisa=y(Isisa)

plot(xsisa,ysisa,'*r')