System of non-linear differential equations with "guess". I'm currently reading about solving nonlinear differential equations, and I've stumpled upon an example that i can't seem to understand. Or rather a part of it.
Consider the system:
\begin{align}
x'&=x+y^2\\
y'&=-y
\end{align}
Solve explicitly, the second equation yields $y(t)=y_0e^{-t}$, inserting this in the first yields
\begin{align}
x'=x+y^2_0e^{-2t}.
\end{align}
The author of the book now guesses a particular solution to be of the form $ce^{-2t}$. And inserting this guess into the equation yields the particular solution:
$x(t)=-\frac{1}{3}y_0^2e^{-2t}$, he then writes that any function of the form
\begin{align}
x(t)=ce^t-\frac{1}{3}y_0^2e^{-2t}
\end{align}
And the general solution is then given by
\begin{align}
x(t)&=\bigg(x_0+\frac{1}{3}y_0^2\bigg)e^t-\frac{1}{3}y^2_0e^{-2t}\\
y(t)&=y_0e^{-t}
\end{align}.
My question is now, where does he insert the guess $ce^{-2t}$, is it $x=ce^{-2t}$, if so, i can't seems to get the particular solution he then presents? And how does he go from the particular solution to the general solution?
A last but less important question: Is there any good litterature or method about "guessing" the particular solution?
Source: Start of Chapter 8, in Differential Equations, Dynamical Systems & An Introduction to Chaos (2nd Ed) by Morris W. Hirsch, Stephen Smale and Robert L. Devaney.
 A: I believe the author tacitly uses the method of undetermined coefficients.
This is a method for finding particular solutions for (non-homogeneous) linear differential equations (DE) with constant coefficients. 
The method works if the non-homogeneous term itself is a solution of a linear DE with constant coefficients; i.e. the non-homogeneous term is a linear combination of (co)sines, exponentials and polynomials. Note that $x'=x + y_0^2 e^{-2t}$ is such a DE.
Basically the method starts with a good educated guess (an ansatz) for the particular solution. One then fills this "candidate particular solution" in to the DE, and then solves an algebraic equation in order to find the coefficients.
The link I added gives a brief description on how to make the good guess and some helpful examples. 

Addendum:
As an example, lets solve the DE $x' = x + y_0^2 e^{-2t}$. 
Step 1: Find the homogeneous solutions: $x_{h} = c e^t$, $c\in\mathbb{R}$.
Step 2: We want to find one particular solution $x_p$. Since the non-homogeneous term is $y_0^2 e^{-2t}$, we take a particular solution of the form $x_p = b e^{-2t}$ where $b$ is a constant we still need to determine. Plugging in this $x_p$ in the DE gives
$$ 
\begin{align*}
  -2 b e^{-2t} &=  b e^{-2t} + y_0^2 e^{-2t} \\
    -3 b e^{-2t} &= y_0^2 e^{-2t} \\
    b &= -\frac{1}{3}y_0^2.
\end{align*}
$$
Hence $x_p = -\frac{1}{3}y_0^2 e^{-2t}$ is a particular solution, so the general solution is $x = c e^t -\frac{1}{3}y_0^2 e^{-2t} $.
Step 3: If there are initial values, then in the last step you substitute these initial values in the general solution. It is not given explicitly in your question, but I guess that in your problem $x(0)=x_0$. Therefore 
$$
     x_0 = x(0) = c e^0 - \frac{1}{3}y_0^2 e^{-2\cdot 0} = c - \frac{1}{3}y_0^2.
$$
Therefore the solution is 
$$
x= \bigl(x_0 + \frac{1}{3}y_0^2\bigr) e^{t} -\frac{1}{3}y_0^2 e^{-2t}.
$$
