theorem
Papers.AnsariRockel2024.elliptical_rho_conditional_comparison
{r : ℝ}
(hr : r ∈ Set.Icc (-1) 1)
(a b c : ℝ)
(hab : 0 < a ^ 2 + b ^ 2)
(hac : 0 < a ^ 2 + c ^ 2)
:
(((ProbabilityTheory.gaussianReal 0 1).prod (ProbabilityTheory.gaussianReal 0 1)).prod
((ProbabilityTheory.gaussianReal 0 1).prod (ProbabilityTheory.gaussianReal 0 1))).real
{p : (ℝ × ℝ) × ℝ × ℝ | b * p.1.2 ≤ a * p.1.1 ∧ c * p.2.2 ≤ a * (r * p.1.1 + √(1 - r ^ 2) * p.2.1)} = 1 / 4 + Real.arcsin (r * a ^ 2 / (√(a ^ 2 + b ^ 2) * √(a ^ 2 + c ^ 2))) / (2 * Real.pi)