This paper consists of two independent, but related parts. In the first part we show how to use symbolic computation to derive explicit expressions for the first few moments of the number of implicants that a random Boolean function has, or equivalently the number of fixed-dimensional subcubes contained in a random subset of the $n$-dimensional cube. These explicit expressions suggest, but do not prove, that these random variables are always asymptotically normal. The second part presents a full, human-generated proof, of this asymptotic normality, first proved by Urszula Konieczna.