%%
%% TP3 Exercice 3
%% Approximation numérique du volume d'un convexe
%%

N=input('Entrer la valeur du paramètre N: ');
X=2*rand(N,1)-1;
Y=2*rand(N,1)-1;
Z=2*rand(N,1)-1;

S=(X.^2 + Y.^2 + Z.^2 < 1);
I=8*sum(S)/N;


plot3(X.*S,Y.*S,Z.*S,'d');


