use numerically stable quadratic form in simple_efficiency

Author: aman-coder03

pvlib/transformer.py

@@ -114,7 +114,8 @@ def simple_efficiency(
     b = 1
     c = no_load_loss - input_power_normalized
 
-    output_power_normalized = (-b + (b**2 - 4*a*c)**0.5) / (2 * a)
+    disc = (b*b - 4*a*c)**0.5
+    output_power_normalized = 2*c / (-b - disc)
 
     output_power = output_power_normalized * transformer_rating
     return output_power