How to find the orthogonal trajectories of the family of all the circles through the points $(1,1)$ and $(-1,-1)$? I'm trying to find the orthogonal trajectories of the family of circles through the points $(1,1)$ and $(-1, -1)$. Now such a family can be given by an equation of the form $$ x^2 + y^2 + 2g(x-y) - 2 = 0, $$ where $g$ is a parameter. 
Now upon differentiation with respect to x, we obtain $$ 2x + 2y y^\prime + 2g (1 - y^\prime ) = 0, $$ where $y^\prime$ denotes the derivative of $y$ with respect to $x$. Upon dividing out by $2$, we arrive at $$ x + y y^\prime + g(1 - y^\prime ) = 0, $$ from which we get $$ g = \frac{x + y y^\prime}{y^\prime - 1}. $$ Putting this value of $g$ into the equation of the family of circles, we get $$ x^2 + y^2 +2 \frac{x + y y^\prime}{y^\prime - 1} ( x - y ) - 2 = 0, $$ so $$ (x^2 + y^2 -2 ) (y^\prime - 1 ) + 2 (x + y y^\prime ) (x - y)  = 0$$ or $$ (x^2 + y^2 - 2 + 2xy - 2y^2 ) y^\prime + (2x^2 - 2xy - x^2 - y^2 + 2) = 0  $$ or $$ ( x^2 + 2xy - y^2 - 2) y^\prime +  (x^2 - 2xy  - y^2 + 2) = 0, $$  from which we get $$ y^\prime = - \frac{ x^2 - 2xy  - y^2 + 2}{ x^2 + 2xy - y^2 - 2}. $$ Now for the orthogonal trajectories, we get $$ y^\prime = \frac{x^2 + 2xy - y^2 - 2}{x^2 - 2xy  - y^2 + 2}. $$ How to solve this differential equation? 
 A: Let $u = y+x$, $v = y-x$, we have
$$\begin{align}
u' = y'+1 = & \frac{2x^2-2y^2}{x^2 - 2xy - y^2+2} = \frac{-2uv}{x^2-2xy-y^2+2}\\
v' = y'-1 = & \frac{4xy-4}{x^2 -2xy- y^2+22} = \frac{u^2 - v^2-4}{x^2-2xy-y^2+2}\\
\end{align}$$
From this, we get
$$u^2 \left( \frac{v^2}{u} + u +\frac{4}{u} \right)' = u^2\left(\frac{v^2}{u}\right)' + (u^2 - 4)u' = 2uvv' + (u^2 - v^2 - 4)u' = 0$$
and hence for some integration constant $K$, one has
$$\frac{v^2}{u} + u +\frac{4}{u} = 2K
\;\;\iff\;\; v^2 + u^2 + 4 = 2K u
\;\;\iff\;\; x^2 + y^2 + 2 = K(x+y)
$$
i.e the orthogonal trajectories is another family of circles.
Update
There is a pure geometric argument why the orthogonal trajectories are circles.
Let $p = (-1,-1), q = (1,1)$. Let $\mathscr{C}$ be a circle centered at $p$ with radius
$|pq| = 2\sqrt{2}$. Consider the inversion $\mu$ with respect to circle $\mathscr{C}$, we have:
$$\mu(p) = \infty\quad\text{ and }\quad \mu(q) = q.$$
This means under $\mu$, the family of circles passing through $p$ and $q$ get mapped to the family of straight lines passing through $q$. The orthogonal trajectories of this family of straight lines is the family of circles centered at $q$. 
It is known that inversion is conformal. i.e. preserve angles and hence map orthogonal trajectories to orthogonal trajectories. This means the orthogonal trajectories of the family of circles passing through $p$ and $q$ are the inverse image of $\mu$ of the family of circles centered at $q$.
It is also known that inversion $\mu$ map circles not passing through $q$ to circles.
As a result, the family of orthogonal trajectories we seek are circles themselves.
A: Assuming your calculations are correct up to the last part, I begin with:
$$ y^\prime = \frac{x^2 + 2xy - y^2 - 2}{x^2 - 2xy  - y^2 + 2}. $$
Observe that:
$$ y^\prime = \frac{1 + 2y/x - (y/x)^2 - 2/x^2}{1 - 2y/x  - (y/x)^2 + 2/x^2}. $$
Thus, let $v = y/x$ for which $y=vx$ and $y' = xv'+v$ and the given problem changes to:
$$ xv'+v = \frac{1 + 2v - v^2 - 2/x^2}{1 - 2v  - v^2 + 2/x^2}. $$
I suspect this can be solved. But, I'll stop here for now. This is a standard substitution for ODEs which can be written as a function of $y/x$ which reflects a certain symmetry, namely that $x$ and $y$ can be scaled by the same factor.
