Yes. It can be done but as you mentioned it is non-elegant.
Perform step (1) to (4).
Then by power of a point, do $(3).(4) = (1).X_b$
After simplification, you should get the following quadratic equation in X(b) , if my calculation did not make any mistake.
$X_{b}^2 (L^2 – Y^2) + 2Y^2XX_{b} –(X^2 + Y^2 – L^2)Y^2 = 0$
Edit: Let $K = X^2 + Y^2 – L^2$.
Then, $[√K – √(X_b^2 + Y^2)][√(X_b^2 + Y^2)] = [X.X_b – X_b^2]$
$√K√(X_b^2 + Y^2) – (X_b^2 + Y^2) = [X.X_b – X_b^2]$
$√K√(X_b^2 + Y^2) = [X.X_b – X_b^2] + (X_b^2 + Y^2)$
$√K√(X_b^2 + Y^2) = [X.X_b + Y^2]$
$K(X_b^2 + Y^2) = X^2.X_b^2 + 2XY^2X_b+ Y^4$
$(X^2 – K)X_b^2+ 2XY^2X_b + Y^2(Y^2 – K)= 0$
Putting K back to get
$(L^2 – Y^2)X_b^2 + 2XY^2X_b + Y^2(L^2 – X^2)= 0$