@@ -318,16 +318,15 @@ set MS12 := (S1 * S2)%R.
318318set C := ((((y * s1 ^+ 2)%R + (m1 * s2 ^+ 2)%R)%R - m2 * s1 ^+ 2) / DS12)%R.
319319under eq_integral do rewrite expRD EFinM.
320320rewrite ge0_integralZr//=; last first.
321- apply/measurable_EFinP.
322- apply: measurableT_comp => //.
323- apply: measurable_funM => //.
324- apply: measurableT_comp => //.
325- apply: (@measurableT_comp _ _ _ _ _ _ (fun t : R => t ^+ 2)%R) => //.
326- exact: measurable_funD.
321+ apply/measurable_EFinP.
322+ apply: measurableT_comp => //.
323+ apply: measurable_funM => //.
324+ apply: measurableT_comp => //.
325+ apply: (@measurableT_comp _ _ _ _ _ _ (fun t : R => t ^+ 2)%R) => //.
326+ exact: measurable_funD.
327327rewrite /normal_peak /normal_fun.
328328rewrite [in RHS]EFinM.
329- rewrite [in RHS]sqr_sqrtr//; last first.
330- by rewrite addr_ge0// sqr_ge0.
329+ rewrite [in RHS]sqr_sqrtr//; last by rewrite addr_ge0// sqr_ge0.
331330rewrite muleA; congr *%E; last by rewrite -mulNr.
332331(* gauss integral *)
333332have MS12DS12_gt0 : (0 < MS12 / DS12)%R.
@@ -340,14 +339,11 @@ transitivity (((Num.sqrt S1 * Num.sqrt S2 * pi *+ 2)^-1)%:E
340339 congr *%E.
341340 apply: eq_integral => x _.
342341 rewrite -EFinM; congr EFin.
343- rewrite normal_pdfE; last first.
344- apply: lt0r_neq0.
345- by rewrite sqrtr_gt0.
342+ rewrite normal_pdfE; last by rewrite lt0r_neq0// sqrtr_gt0.
346343 rewrite mulrA mulVf// ?mul1r//.
347344 rewrite lt0r_neq0// invr_gt0 sqrtr_gt0 pmulrn_lgt0// mulr_gt0// ?pi_gt0//.
348- rewrite exprn_even_gt0//=.
349- by rewrite lt0r_neq0// sqrtr_gt0.
350- rewrite ge0_integralZl//; last 3 first.
345+ by rewrite exprn_even_gt0//= lt0r_neq0// sqrtr_gt0.
346+ rewrite ge0_integralZl//=; last 3 first.
351347 apply/measurable_EFinP.
352348 exact: measurable_normal_pdf.
353349 move=> x _.
0 commit comments