There is a much easier way to solve the problem. Quick derivation: if something is moving past you at speed $u = \alpha~c$ then in a reference frame travelling at the speed $\beta ~c$ in that direction, the Lorentz boost puts its trajectory as:$$\gamma\begin{bmatrix}1&-\beta\\-\beta&1\end{bmatrix}\begin{bmatrix}c t\\ut\end{bmatrix} = \gamma \begin{bmatrix}c t - \beta \,ut\\ut - \beta\,ct\end{bmatrix} $$and substituting in you find a speed $c~\frac{\alpha - \beta}{1 - \alpha \beta}.$ But suppose we define the hyperbolic sine and cosine as $\sinh x = (e^x - e^{-x})/2,$ $\cosh x = (e^x + e^{-x})/2.$ Then we find that the products of these various things are$$\begin{align}
\sinh a\,\sinh b =& \frac 14 \left(e^{a + b} - e^{a - b} - e^{-a + b} + e^{-a - b} \right),\\
\sinh a\,\cosh b =& \frac 14 \left(e^{a + b} + e^{a - b} - e^{-a + b} - e^{-a - b} \right),\\
\cosh a\,\sinh b =& \frac 14 \left(e^{a + b} - e^{a - b} + e^{-a + b} - e^{-a - b} \right),\text{ and}\\
\cosh a\,\cosh b =& \frac 14 \left(e^{a + b} + e^{a - b} + e^{-a + b} + e^{-a - b} \right).
\end{align}$$This yields something very much like the trigonometric "angle sum rules" when you try to cancel out the $e^{\pm a \mp b}$ terms, $$\begin{align}\sinh(a + b) =& \sinh a\,\cosh b + \cosh a\,\sinh b \\
\cosh(a + b) =& \cosh a\,\cosh b + \sinh a\,\sinh b
\end{align}$$For $a - b$, simply replace both $+$ signs with $-$ signs on the right hand side, as $\sinh$ is odd and $\cosh$ is even.
Now rewrite each of the above expressions as $\alpha = \sinh a / \cosh a,\; \beta = \sinh b / \cosh b$. It turns out that this function $\tanh x$ is totally invertible, so you can always do this. Then multiply both top and bottom of the fraction by $\cosh a \, \cosh b$ to find: $$c \frac{\alpha - \beta}{1 - \alpha\beta} = c~\frac{\sinh a\,\cosh b - \cosh a\,\sinh b}{\cosh a\,\cosh b - \sinh a\,\sinh b}= c~\frac{\sinh(a - b)}{\cosh(a - b)} = c \tanh(a - b)$$So when you see something moving at a speed $u = c~\tanh a$ and you boost into a reference frame with speed $c~\tanh b$ relative to where you started, you find that the new velocity is just $u' = c~\tanh(b - a).$
These numbers $w = \tanh^{-1}(u/c)$ are usually called rapidities. In special relativity, rapidities are linear. So for example, if you want to know what "constant acceleration" looks like in the relativistic limit, it is $v = c~\tanh(a t / c).$
So, find the rapidity, divide it by 2, boost into the frame with that rapidity, and you should find that both of these particles are coming towards you with the same rapidity.