Nous résolvons dans \( \mathbb{R} \times \mathbb{R} \), par la méthode de substitution, le système

\[
(E) 
\begin{cases} 
x + 3y = 9 \quad \text{(1)} \\ 
2x - y = 1 \quad \text{(2)} 
\end{cases}
\]

- Résolution par substitution :

1. **Étape 1** : Exprimer \( x \) en fonction de \( y \) à partir de l'équation (1) :
\[
x = 9 - 3y \quad \text{(3)}
\]

2. **Étape 2** : Substituer \( x \) dans l'équation (2) par l'expression trouvée en (3) :
\[
2(9 - 3y) - y = 1
\]
\[
18 - 6y - y = 1
\]
\[
18 - 7y = 1
\]

3. **Étape 3** : Résoudre pour \( y \) :
\[
-7y = 1 - 18
\]
\[
-7y = -17
\]
\[
y = \frac{17}{7}
\]

4. **Étape 4** : Substituer \( y = \frac{17}{7} \) dans l'équation (3) pour trouver \( x \) :
\[
x = 9 - 3\left(\frac{17}{7}\right)
\]
\[
x = 9 - \frac{51}{7}
\]
\[
x = \frac{63}{7} - \frac{51}{7}
\]
\[
x = \frac{12}{7}
\]
La solution du système est :

\[
\boxed{(x, y) = \left( \frac{12}{7}, \frac{17}{7} \right)}
\]