The other answers are good, I just wanted to get the equation in a different way.
Starting with $a x^2 + b x y + c y^2 + d x + e y + f = 0$, right off the bat $f = 0$ since $(0,0)$ is on the curve and $d = 0$ since it's a tangency there and the tangent cone should be $e y$. Now $y = x + 1$ is along the major axis, so we get three more points $(-1,1), (6,8), (7,7) $ on the curve and the linear system:
$a - b + c + e = 0$
$49 a + 49 b +49 c +7e = 0$
$36 a + 48 b + 64 c + 8 e = 0$
We have the freedom to choose one parameter, let's use $e=-7$ to simplify the second equation:
$a - b + c - 7 = 0$
$a + b + c - 1 = 0$
$9 a + 12 b + 16 c - 14 = 0$
Now it's easily got that $a=2, b=-3, c=2$ and we have the equation $2x^2-3x y+2y^2-7y=0$.
Let's make that $4x^2-6x y+4y^2-14y=0$ and apply standard rotation techniques:
$\begin{pmatrix}x&y\end{pmatrix}\begin{pmatrix}4 &-3\\ -3 & 4\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}+\begin{pmatrix}0&-14\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}=0$, and we get eigenvalues 7 and 1 and $P=\frac1{\sqrt{2}}\begin{pmatrix}-1&1\\1&1\end{pmatrix}$, ${\bf x}^tP^t A P{\bf x}+KP{\bf x}=0$. That is $7x'^2+y'^2-7\sqrt{2}x'-7\sqrt{2}y'=0$ and completing squares we get $7(x'-\frac1{\sqrt{2}})^2+(y'-\frac{7}{\sqrt{2}})^2=28$ or $(\frac{x'-\frac1{\sqrt{2}}}{2})^2+(\frac{y'-\frac{7}{\sqrt{2}}}{2\sqrt{7}})^2=1$ Now we have the major axis $2\sqrt{7}$ and the minor axis $2$, so distance from focus to center is $\sqrt{(2\sqrt{7})^2-2^2}=\sqrt{24}$ and the foci are $(3+2\sqrt{3},4+2\sqrt{3})$ and $(3-2\sqrt{3},4-2\sqrt{3})$ from two 45-45-90 triangles with $\sqrt{24}$ as hypotenuse.