Analytical solution to linear coupled ODEs? I have this set of equations and I am trying to find $X(t)$ and $Y(t)$ analytically with initial values known such as $X(0)=X_0$ and $Y(0)=Y_0$. How should I approach to solve it?
$dX/dt=a \times X(t) + b \times Y(t)$
$dY/dt=c \times Y(t) + d \times X(t)$

Write your system of equations in matrix form:
$$\begin{pmatrix} \frac{dX}{dt} \\ \frac{dY}{dt} \end{pmatrix} =\begin{pmatrix} a & b \\ d & c \end{pmatrix} \begin{pmatrix} X \\ Y \end{pmatrix}$$
You can find the two eigenvalues $\lambda_1$ and $\lambda_2$ by letting $\det(A-\lambda I)=0$, and then evaluate corresponding eigenvectors $\vec{v_1}$ and $\vec{v_2}$ of the matrix $A=\begin{pmatrix} a & b \\ d & c \end{pmatrix}$ to give you the general solution to $X(t)$ and $Y(t)$.
Depending on your values of your eigenvalues $\lambda_1$ and $\lambda_2$ (which ultimately depend on the values of your constants $a,b,c,d$), there will be different general solutions.
There are 3 different solution cases:
Case 1: Real eigenvalues: $\lambda \in \mathbb{R}$
The solution will be of the form:
$$\begin{pmatrix} X(t) \\ Y(t) \end{pmatrix} =k_1 \vec{v_1} e^{\lambda_1 t}+k_2 \vec{v_2} e^{\lambda_2 t}$$
Case 2: Complex eigenvalues.
We denote a complex number by $\lambda=a+bi$.
Since $\lambda_{1,2}$ will both be complex, using $e^{i \theta}=\cos{\theta}+i\sin{\theta}$, we obtain the general solution:
$$\begin{pmatrix} X(t) \\ Y(t) \end{pmatrix}=k_1 \vec{v_1} e^{at} \cos{(bt)} + k_2 \vec{v_2} e^{at} \sin{(bt)}$$
Case 3: Real, repeated roots (i.e. $\lambda_1=\lambda_2$)
$$\begin{pmatrix} X(t) \\ Y(t) \end{pmatrix}=k_1 \vec{v_1} e^{\lambda t}+k_2 \vec{v_2} \cdot t e^{\lambda t}$$
Your initial values can be used to evaluate the constants $k_1$ and $k_2$.
Here is an Youtube video showing an example solve for Case 1.
