Analytic solution of coupled nonlinear differential equation

I would like to solve the following differential equation system,

$$\frac{dx}{dt} = C_{1}\cdot x(t) + C_{2}\cdot x(t)^{2} +C_{3}\cdot x(t)\cdot y(t)\\ \frac{dy}{dt} = C_{4}\cdot y(t) + C_{5}\cdot y(t)^{2} +C_{6}\cdot x(t)\cdot y(t)$$

where $$C_{1},C_{2},C_{3},C_{4},C_{5},C_{6}$$ are constants and the initial conditions are given, $$x(t=0) = xt_{0}\qquad y(t=0) = yt_{0}$$ and the values of the variables in equilibrium when $$\frac{dx}{dt} = 0 \qquad \frac{dy}{dt} = 0$$ are also given: $$x_{0} = 0 \qquad y_{0} = 0$$

Could someone give me any advise or method to help me find the analytic solution of x(t), y(t).

• That is a quite interesting problem. I'm afraid that there is no analytical solution but for some special cases (perhaps for $C_2=C_3=C_5=C_6=0$). Since the RHS of your equations is not Lipschitz continuous for all $t \in R$, the Picard-Lindelöf theorem can not ensure uniqueness (nor existence) for solutions for infinite time. Considering non uniqueness, but existence, maybe it means that there is some stationary solution for infinite $t$ in the form of an attractor to which more then one solution converges. – rafa11111 Apr 17 '18 at 13:55
• @rafa11111 The RHS is locally Lipschitz in the variables $x,y$. Existence and uniqueness of a local solution is guaranteed. – Julián Aguirre Apr 17 '18 at 14:06
• $(0,0)$ is an equilibrium point. Its character (sink, source, saddle,...) will depend on the constants $C_i$, and can be found by linearization. Also, depending on the values of the constants, there may be more equilibria. – Julián Aguirre Apr 17 '18 at 14:22