R=5;r=1;L=R+r; %lähtöarvot
s=0; %apumuuttuja
K=10^4; %arvausten lukumäärä
for k=1:K %aloitetaan silmukka ja
x=L*rand(1); %arvotaan piste (x,y,z)
y=L*rand(1);
z=r*rand(1);
if le(abs(sqrt(x^2+y^2)-R)^2 +z^2, r^2),
s=s+1; %lasketaan toruksen sisällä olevat pisteet
else
end
end
osuus=s/K; %lasketaan osuus
Vt=8*osuus*L^2*r %lasketaan tilavuus