// Exercice 4

X1=grand(1,1000,'nor',0,1);
X2=grand(1,1000,'nor',0,1);
M=[2 -1;1 -1];
N=1/M;
Z=N*[X1;X2];
Y1=Z(1,:);
Y2=Z(2,:);
l1=Y1>=0.4;
l2=Y2<=3;
l=bool2s(l1).*bool2s(l2);
p=sum(l)/1000;
