sqrt((((((b + a)) * (a - (- b))) * (a + (- b))) * ( (b - a)))) / 4.0

(FPCore (b a)
 :precision binary64
 (/ (sqrt (* (* (* (+ b a) (- a (- b))) (+ a (- b))) (- b a))) 4.0))

double code(double b, double a) {
	return sqrt(((((b + a) * (a - -b)) * (a + -b)) * (b - a))) / 4.0;
}

real(8) function code(b, a)
    real(8), intent (in) :: b
    real(8), intent (in) :: a
    code = sqrt(((((b + a) * (a - -b)) * (a + -b)) * (b - a))) / 4.0d0
end function

public static double code(double b, double a) {
	return Math.sqrt(((((b + a) * (a - -b)) * (a + -b)) * (b - a))) / 4.0;
}

def code(b, a):
	return math.sqrt(((((b + a) * (a - -b)) * (a + -b)) * (b - a))) / 4.0

function code(b, a)
	return Float64(sqrt(Float64(Float64(Float64(Float64(b + a) * Float64(a - Float64(-b))) * Float64(a + Float64(-b))) * Float64(b - a))) / 4.0)
end

function tmp = code(b, a)
	tmp = sqrt(((((b + a) * (a - -b)) * (a + -b)) * (b - a))) / 4.0;
end

code[b_, a_] := N[(N[Sqrt[N[(N[(N[(N[(b + a), $MachinePrecision] * N[(a - (-b)), $MachinePrecision]), $MachinePrecision] * N[(a + (-b)), $MachinePrecision]), $MachinePrecision] * N[(b - a), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / 4.0), $MachinePrecision]

\begin{array}{l}

\\
\frac{\sqrt{\left(\left(\left(b + a\right) \cdot \left(a - \left(-b\right)\right)\right) \cdot \left(a + \left(-b\right)\right)\right) \cdot \left(b - a\right)}}{4}
\end{array}