Starting of your definition in the question if $F_1(-c,0)$ and $F_2(c,0)$ are the foci ($c>0$) and $P\in\mathcal E$ with $P(x,y)$ and $O=(0,0)$ you have
$$\mathcal E: \quad |PF_1|+|PF_2|=2a$$
Hence:
$$(x+c)^2+y^2+(x-c)^2+y^2+2\sqrt{[(x+c)^2+y^2][(x-c)^2+y^2]}=4a^2$$
$$2x^2+2c^2+2y^2+2\sqrt{[(x+c)^2+y^2][(x-c)^2+y^2]} = 4a^2$$
$$x^2+c^2+y^2+\sqrt{[(x+c)^2+y^2][(x-c)^2+y^2]} = 2a^2$$
Isolate the root and square it again....
$$\sqrt{[(x+c)^2+y^2][(x-c)^2+y^2]} = 2a^2-x^2-y^2-c^2$$
$$(x^2+c^2+2xc+y^2) (x^2+c^2-2xc+y^2) = 4a^4+x^4+y^4+c^4-4a^2x^2-4a^2y^2-4a^2c^2+2x^2y^2+2x^2c^2+2y^2c^2$$
$$x^4+c^4+y^4+2x^2c^2+2x^2y^2+2c^2y^2-4x^2c^2 =4a^4+x^4+y^4+c^4-4a^2x^2-4a^2y^2-4a^2c^2+2x^2y^2+2x^2c^2+2y^2c^2$$
From which, simplifying,
$$4a^4-4a^2x^2-4a^2y^2-4a^2c^2+4x^2c^2 = 0$$
$$a^4-a^2x^2-a^2y^2-a^2c^2+x^2c^2 = 0$$
$$x^2(c^2-a^2)-a^2y^2 = a^2c^2 -a^4$$
Finally, remembering that $a^2-c^2 \gt 0$ we can put $b^2=a^2-c^2$ and have like this
$$-b^2x^2-a^2y^2=-b^2a^2$$ hence the proof
$$\begin{equation} \frac{x^2}{a^2}+\frac{y^2}{b^2} = 1 \end{equation}$$