import Mathlib

#check Real.sqrt_eq_zero_of_nonpos
#check div_zero

theorem sqrt_negative_is_zero : Real.sqrt (-3 : ℝ) = 0 := by
  simp [Real.sqrt_eq_zero_of_nonpos]

theorem gamma_relation_fails_at_two :
    (1 / Real.sqrt (1 - (2 : ℝ)^2))^2 * (1 - (2 : ℝ)^2) = 0 := by
  norm_num [Real.sqrt_eq_zero_of_nonpos]