I'm studying the spectral theorem as appears in Reed and Simon's Functional Analysis.
Assume we have constructed the continuous functional calculus for a self adjoint bounded operator $A$ on a hilbert space $H$.
The construction of the Borel functional calculus is through a function $\Phi : \mathbb{B}(\mathbb{R}) \to L(H)$ from the bounded Borel functions to bounded linear operators on $H$, which extends the continuous functional calculus.
Please see this link for the explicit construction (page 4 of the pdf):
http://www.math.mcgill.ca/jakobson/courses/ma667/mendelsontomberg-spectral.pdf
I want to show $AB = BA \implies B\Phi(f) = \Phi(f)B$.
The claim is clearly true for polynomials (for $f(x) = x, \Phi(f) = A$), and hence for continuous functions, being the uniform limit of polynomials and as $\Phi$ is continuous under the supremum norm.
My attempt:
Step 1 is to prove $\forall x \in H (x, B\Phi(f) x) = (x, \Phi(f) B x)$.
Step 2 is to notice that the above implies $\forall x,y \in H$ $(x, B\Phi(f) y) = (x, \Phi(f) B y)$ by the polarization identity.
Step 3 the above of course implies that $\Phi(f) B = B\Phi(f)$.
So we only need a proof for 1:
Set $x \in H$.
Define $\phi_k = B^*x +i^kx$ for $k \in \{0,1,2,3\}$, $\psi_k = x +i^kBx$ for $k \in \{0,1,2,3\}$.
Define $\mu = \sum_{k = 0}^{3} |\mu_{\phi_k}| + |\mu_{\psi_k}|$. Where $\mu_{\phi_k}$, $\mu_{\psi_k}$ are the spectral measures associated to the corresponding vectors. As each is a regular Borel measure it follows that $\mu$ is a regular Borel measure.
Choose $\{f_n\} \subset C(\sigma(A))$ s.t $\int f_n d\mu \to \int f d\mu$. It follows that $\int f_n \mu_{\psi_k} \to \int f \mu_{\psi_k}$, and $\int f_n \mu_{\phi_k} \to \int f \mu_{\phi_k}$ for all $k \in \{0,1,2,3\}$.
So we have that (see the Mcgill link (page 4) for the second and sixth equalities):
$(x, B\Phi(f)x) = (B^*x, \Phi(f)x) \overset{definition}{=} \frac{1}{4} \sum_{k = 0}^{3} i^k \int f d\mu_{\phi_k} = lim_n \frac{1}{4} \sum_{k = 0}^{3} i^k \int f_n d\mu_{\phi_k} = lim_n(B^*x, \Phi(f_n)x) \overset{f_n \in C(\sigma(A))}{=} lim_n(x, \Phi(f_n)Bx) = lim_n \frac{1}{4} \sum_{k = 0}^{3} i^k \int f_n d\mu_{\psi_k} = \frac{1}{4} \sum_{k = 0}^{3} i^k \int f d\mu_{\psi_k} = (x, \Phi(f)Bx)$.
Btw, also posted this question on MSE, but the methods suggested draw from backgrounds I don't have yet, and I'd like to get feedback for this answer.