Starting from the generalized form of my original formula
##\arcsin(\frac{ax}{2})+ \arcsin(\frac{bx}{2})+ \arcsin(\frac{cx}{2})+ \arcsin(\frac{dx}{2})=\pi. ##
we get;
##\sqrt{\left(1-\frac{a^2x^2}{4}\right) \left(1-\frac{b^2x^2}{4}\right)}\\ \quad+ \sqrt{\left(1-\frac{c^2x^2}{4}\right) \left(1-\frac{d^2x^2}{4}\right)} =\frac{(ab+cd)x^2}{4}##
Define. ##u=\frac{x^2}{4}##
For convenience, let
##F=\sqrt{(1-a^2u)(1-b^2u)} ##
and
##G=\sqrt{(1-c^2u)(1-d^2u)}##
Our equation is now simply
##F+G=(ab+cd)u##
Squaring gives
##F^2+2FG+G^2=(ab+cd)^2u^2##
Expand the individual terms:
## F^2=1-(a^2+b^2)u+a^2b^2u^2## and ##G^2=1-(c^2+d^2)u+c^2d^2u^2##
Substituting and rearranging,
##2FG=-2+(a^2+b^2+c^2+d^2)u+2abcd\,u^2##
Therefore,
## FG=-1+\frac{a^2+b^2+c^2+d^2}{2}u+abcd\,u^2##
Squaring the last equation gives
##F^2G^2= \left[-1+\frac{a^2+b^2+c^2+d^2}{2}u+abcd\,u^2\right]^2##
Substitute the expressions for ##F^2 ##and## G^2##:
##(1-a^2u)(1-b^2u)(1-c^2u)(1-d^2u)\\ =\left[-1+\frac{a^2+b^2+c^2+d^2}{2}u+abcd\,u^2\right]^2##
This has no square roots but is a bit of a challenge to expand and simplify!
1. Left-Hand Side (LHS)
##\text{LHS} = (1 - a^2u)(1 - b^2u)(1 - c^2u)(1 - d^2u)##
Expanded out term-by-term in powers of ##u##
##u^0## 1
##u^1##:## - (a^2 + b^2 + c^2 + d^2)##
##u^2##: ##+ (a^2b^2 + a^2c^2 + a^2d^2 + b^2c^2 + b^2d^2 + c^2d^2)##
##u^3##:## - (a^2b^2c^2 + a^2b^2d^2 + a^2c^2d^2 + b^2c^2d^2)##
##u^4##:## + a^2b^2c^2d^2 u^4##
2. Right-Hand Side (RHS)
##\text{RHS} = \left[-1 + \frac{a^2+b^2+c^2+d^2}{2}u + abcd,u^2\right]^2##
Expanded out term-by-term in powers of ##u##]:
##u^0##:## = 1 ##
##u^1##:##- (a^2 + b^2 + c^2 + d^2)##
##u^2##:##+ \left[\frac{(a^2+b^2+c^2+d^2)^2}{4} - 2abcd\right]##
##u^3##:##+ abcd(a^2 + b^2 + c^2 + d^2)##
##u^4##:##+ a^2b^2c^2d^2 ##
The zeroth, first and fourth powers of ##u## cancel out. Then ##u## will simply be (haha) the difference of the coefficients of the second power divided by the difference of the coefficients of the third power of ##u##.
##u =\Large \frac{2abcd + \frac{1}{2}\left(a^2b^2 + a^2c^2 + a^2d^2 + b^2c^2 + b^2d^2 + c^2d^2\right) - \frac{1}{4}\left(a^4 + b^4 + c^4 + d^4\right)}{a^2b^2c^2 + a^2b^2d^2 + a^2c^2d^2 + b^2c^2d^2 + abcd\left(a^2 + b^2 + c^2 + d^2\right)}##
The numerator factors as;
##N=\frac{-1}{4}(a+b+c−d)(a+b−c+d)(a−b+c+d)(a−b−c−d)##
or using the semiperimeter;
##4(s-a)(s-b)(s-c)(s-d)##
The denominator factors as;
##D=(ab+cd)(ac+bd)(ad+bc)##
The complete form is;
##u=\Large\frac{4(s-a)(s-b)(s-c)(s-d)}{(ab+cd)(ac+bd)(ad+bc)}##
Now, ##u=\frac{x^2}{4}=\frac{1}{4r^2}##
Thus ##r=\sqrt{\frac{1}{4u}}##
This gives us ##r=\Large\frac{1}{4}\sqrt{\frac{(ab+cd)(ac+bd)(ad+bc)}{(s-a)(s-b)(s-c)(s-d)}}##
I think I prefer
@Charles Link’s method!