Let $\left(\xi_n\right)_{n\geqslant1}$ be a sequence of independent identically distributed random variables with distribution $\mathcal{N}(0,1)$. Let $X_0$ be a random variable in $\mathbb{R}$ with distribution $\mu$, independent from the sequence $\left(\xi_n\right)_{n\geqslant1}$. Let us define the sequence $\left(X_n\right)_{n\geqslant1}$ by $$X_n=\alpha X_{n-1}+\xi_n$$ where $\alpha\in(0,1)$.
Show that $\left(X_n\right)_{n\geqslant0}$ is a Markov chain and compute its transition kernel. Find $\sigma^2>0$ such that $\mathcal{N}\left(0,\sigma^2\right)$ is a stationnary distribution for this chain.
My work
I have computed $X_n=\alpha^nX_0+\sum_{j=1}^n \alpha^{n-j}\xi_j$, however when writing $$\mathbb{P}\left(X_0\in C_0,\dots,X_n\in C_n\right)$$ for $C_0,\dots,C_n$ borel sets of $\mathbb{R}$, I do not know how to use independence to make it look like $$\int_{\substack{x_0\in C_0\\\dots\\x_n\in C_n}}\mu(dx_0)P(x_0,dx_1)\dots P(x_{n-1},dx_n)$$