$($c$)$ Using Monte Carlo simulation, write an algorithm

to calculate that part of the volume of the sphere

$x^2+y^2+z^{2}1$ that lies in the first octant,

$x>0,y>0,z>0.$