Some definitions
Let $(X,\mathbb X,\mu)$ be a measure space. A real-valued function is simple if it has only a finite number of values. A simple $\mathbb X$-measurable function $\varphi$ can be represented in the form $$\varphi=\sum_{j=1}^n a_n\chi_{E_j}$$ where $a_j\in\mathbb R$ and $\chi_{E_j}$ is the characteristic function of a set $E_j\in\mathbb X$. If we add the restriction that the $a_j$ be distinct and the $E_j$ form a partition of X, then the representation is unique and is called the standard representation of $\varphi$.
If $\varphi$ is a simple function in $M^+(X,\mathbb X)$ with the standard representation above, we define the integral of $\varphi$ with respect to $\mu$ to be the extended real number $$\int\varphi\,d\mu=\sum_{j=1}^n a_j\mu(E_j)$$
My question
If the simple function $\varphi\in M^+(X,\mathbb X)$ has the (not necessarily standard) representation $$\varphi=\sum_{k=1}^m b_k\chi_{F_k}$$ where $b_k\in\mathbb R$ and $F_k\in\mathbb X$, it can be shown that $$\int\varphi\,d\mu = \sum_{k=1}^m b_k\,\mu(F_k).$$ My problem is that I cannot find a clean yet rigorous step-by-step proof of that result.
My idea is to rewrite the function $\varphi$ as $\sum_{k=1}^n a_k\chi_{\phi^{-1}(a_k)}$ where $a_k$ is the sum of some $b_k$ terms and $\phi^{-1}(a_k)$ is the union of intersections of some $F_k$ terms. After some manipulations I note that I can put back together all the "pieces" and find $\sum_{k=1}^m b_k\,\mu(F_k)$. Unfortunately, that very last passage is left to the reader. Is there a way to make all the process explicit and yet clean and easy to follow?

