Solve system of equations $ y + 3x = 4x^3$ and $x + 3y = 4y^3$ Problem: In the set of real numbers find all $x,y$ that satisfy
$ y + 3x = 4x^3,x + 3y = 4y^3$.
What I don't understand, is how do I know that a condition satisfies both equations(I will show later exactly).
The solution goes like this: By addition and substraction of the equations we get the following 2 equations:
$x + y = (x + y)(x^2 − xy + y^2),x − y = 2(x − y)(x^2 + xy + y^2)$.
My question is, how do I know which conditions have to be true? According to the solution either $x+y=0$ or $x-y=0$, which I understand. My problem comes when they mention that both $(x^2 + xy + y^2)=1/2$ and $(x^2 - xy + y^2)=1$ have to be true at the same time. Why both? How do systems of equations work and how do we know that there won't be another solution if we multiply or divide the equations? Any tips for these kinds of problems are greatly appreciated.
Question 2: Why can't we also try to find a solution which satisfies for example $x-y=0$ and $(x^2 − xy + y^2)=1$

Here is a natural way to solve the system. Continue with, 
$$x + y = (x + y)(x^2 − xy + y^2),\>\>\>\>\>
x − y = 2(x − y)(x^2 + xy + y^2)$$
or,
$$(x + y)(x^2 − xy + y^2-1)=(x − y)(x^2 + xy + y^2-\frac12)=0$$
So, there are four cases to be considered,
Case 1: Substitute $x-y=0$ into one of the original equations, say, $y + 3x = 4x^3$, which leads to,
$$4x(x^2-1)=0\implies x= 0, \pm 1$$
Case 2:  Similarly, $x+y=0$ leads to,
$$2x(2x^2-1)=0\implies x = 0,\pm \frac 1{\sqrt2}$$
Case 3: Substitute $x^2 − xy + y^2-1=0$ into $y+3x=4x^3$, leading to
$$16x^6-28x^4+13x^2-1=0,\>\> \text{or}$$
$$(x^2-1)(16x^4-12x^2+1)=0\implies x=\pm 1,\pm\frac{\sqrt5\pm 1}{4}$$
Case 4: Similarly $x^2 + xy + y^2-\frac12=0$ leads to
$$(2x^2-1)(16x^4-12x^2+1)=0\implies x=\pm \frac1{\sqrt2},\pm\frac{\sqrt5\pm 1}{4}$$
Thus, all the real solutions are:
$(0,0), (1,1), (-1,-1),
(\frac 1{\sqrt2},-\frac 1{\sqrt2}), (-\frac 1{\sqrt2},\frac 1{\sqrt2}), $
$(\frac {\sqrt5 + 1}{4},\frac {-\sqrt5 + 1}{4}), (\frac {\sqrt5 - 1}{4},\frac {-\sqrt5 - 1}{4}), 
(\frac {-\sqrt5 + 1}{4},\frac {\sqrt5 + 1}{4}), (\frac {-\sqrt5 - 1}{4},\frac {\sqrt5 - 1}{4})$
Note that the system of the two equations are implicitly of ninth order. So, it is expected to have nine pairs of solutions.
