// Exercice 3
m=2;
n=10000;
U=grand(n,m,'nor',0,1);
Y=sum(U.^2,2);
figure 
histplot(20,Y) // histogramme de Y
t1=[0:0.1:20];
plot2d(t1,dchisq(t1,m),2) // densité normale
YY=-gsort(-Y) // rangement par ordre croissant
k=length(YY)
figure 
plot2d(YY,[1:k]/k) // fonction de répartition empirique
t2=[0:0.1:20];
plot2d(t2,pchisq(t2,m),2)
