Creates uniformly distributed random spheres which could be overlapping or not. Calculates following properties: pore connectivity, porosity, specific surface area, pore diameter, surface roughness (this one might not work properly). Plots 3D pore network, spheres, pores. Can add 3rd phase (water droplet), add extra space on top of porous media, record matrix to txt file.