diff --git a/calculations/ECCC.jl b/calculations/ECCC.jl
index 1a9035d..8b4b00a 100644
--- a/calculations/ECCC.jl
+++ b/calculations/ECCC.jl
@@ -44,7 +44,7 @@ target_c = range(0.78, 0.45, 7) # original: range(0.75, 0.40, 8)
 target_E = [quick_pole_E(V_system(c)) for c in target_c]
 
 # Solve coefficients as a linear system
-M_left_element(c, i) = complex(c - c0)^(i/2)
+M_left_element(c, i) = alt_sqrt(c - c0)^i
 M_left = M_left_element.(training_c, (0:order)')
 a = zeros(ComplexF64, dim, order+1)
 b = zeros(ComplexF64, dim, order+1)
@@ -58,7 +58,7 @@ for i in 1:dim
 end
 
 # Pade approximant
-polynomial(a, c) = sum(i -> a[:, i+1] .* complex(c - c0)^(i/2), 0:order)
+polynomial(a, c) = sum(i -> a[:, i+1] .* alt_sqrt(c - c0)^i, 0:order)
 pade_approx(c) = polynomial(a, c) ./ polynomial(b, c)
 
 # Extrapolate