General solution of an inhomogeneous linear equation

Statement

Let $I$ be an open interval, let $a_0, \dots, a_{n-1}, b : I \to \mathbb{R}$ be continuous, and let $L y = y^{(n)} + a_{n-1}(x) y^{(n-1)} + \dots + a_0(x) y$. Let $y_p$ be a solution of $L y = b$ on $I$, and let $(y_1, \dots, y_n)$ be a fundamental system of $L y = 0$ on $I$. Then a function $y : I \to \mathbb{R}$ is a solution of $L y = b$ on $I$ if and only if $$y = y_p + c_1 y_1 + \dots + c_n y_n$$ for some real numbers $c_1, \dots, c_n$.