System of ODE with non constant coefficients I have this system
$$\begin{bmatrix}\frac{dy_1(x)}{dx} \\ \frac{dy_2(x)}{dx} \end{bmatrix}=\begin{bmatrix}
     \frac{1}{x} & \frac{1}{x} \\
     \frac{4}{x} & \frac{1}{x} 
  \end{bmatrix}\begin{bmatrix}y_1(x) \\ y_2(x) \end{bmatrix}$$
I can solve system of ODE with constant coefficients and write the general solution but this is not the case. Does this system have a solution you can write with elementary functions? If this is the case how can I do? Is there a systematic procedure? Could you explain me?
Any help will be appreciated.

Let $A=\begin{bmatrix}1 & 1 \\ 4 & 1\end{bmatrix}$, $y=\begin{bmatrix}y_1 \\ y_2\end{bmatrix}$ and write the equation as
$$
           \frac{dy}{dx}-\frac{1}{x}Ay = 0.
$$
The characteristic polynomial of $A$ is $(\lambda-1)^2-4=(\lambda-3)(\lambda+1)$. So $A$ diagonalizable. It has a basis of eigenvectors
$$
        X_1=\begin{bmatrix}1 \\ 2\end{bmatrix}, \; X_2 =\begin{bmatrix}-1 \\ 2\end{bmatrix}
$$
$AX_1 = 3X_1$ and $AX_2 = -X_2$.Then $y(x)=a(x)X_1+b(x)X_2$ can be substituted into the original equation to obtain
$$
         a'(x)X_1+b'(x)X_2-\frac{3}{x}a(x)X_1+\frac{1}{x}b(x)X_2=0.
$$
This gives component equations
$$
         a'(x)-\frac{3}{x}a(x)=0,\;\; b'(x)+\frac{1}{x}b(x)=0 \\
           \frac{a'}{a}=\frac{3}{x},\;\; \frac{b'}{b}=-\frac{1}{x} \\
           a = Ax^3,\;\; b=B\frac{1}{x}.
$$
The general solution is
$$
      y(x) = Ax^3\begin{bmatrix}1 \\ 2\end{bmatrix}+B\frac{1}{x}\begin{bmatrix}-1 \\ 2\end{bmatrix},
$$
where $A$ and $B$ are arbitrary constants. So $y_1(x)=Ax^3-B/x$ and $y_2(x)=2Ax^3+2B/x$ are your solutions.
