К выводу квадратного уравнения из (37) и его решению (38)

kill(all)$
assume(l[12]>0,l[13]>0,l[23]>0,m[2]>0,m[3]>0);

eq1: sqrt(l[13]^2-m[2]^2*x^2/m[3]^2) - sqrt(l[12]^2-x^2) = sqrt(l[23]^2-(m[2]+m[3])^2*x^2/m[3]^2);

eq2: m[3]^2*ratsimp(eq1^2);

eq3: lhs(eq2) - rhs(eq2) - first(lhs(eq2)) = -first(lhs(eq2));

eq4: factor(eq3^2/m[3]^2);

eq5: rat(lhs(eq4) - rhs(eq4)) = 0;

a: expand(first(lhs(eq5))/x^2);

b: expand((lhs(eq5) - a*x^2)/m[3]^2);

tex(b);

combine((rhs(second(solve(eq4, x)))*2/m[3])^2);