const I1_ef = i1 * Math.sqrt(1 + thd1 ** 2);
    const I2_ef = i2 * Math.sqrt(1 + thd2 ** 2);
    const I3_ef = i3 * Math.sqrt(1 + thd3 ** 2);

    const phiA = Math.acos(cosPhiA); // Ángulo de Ia 0°
    const phiB = Math.acos(cosPhiB) + (2 * Math.PI / 3); // (120° desfasado)
    const phiC = Math.acos(cosPhiC) - (2 * Math.PI / 3); // (-120° desfasado)

    // Convertir corrientes a componentes reales e imaginarias
    const Ia_real = I1_ef * Math.cos(phiA);  const Ia_imag = I1_ef * Math.sin(phiA);
    const Ib_real = I2_ef * Math.cos(phiB);  const Ib_imag = I2_ef * Math.sin(phiB);
    const Ic_real = I3_ef * Math.cos(phiC);  const Ic_imag = I3_ef * Math.sin(phiC);

    // Sumar las componentes reales e imaginarias de las corrientes de fase
    const In_real = -(Ia_real + Ib_real + Ic_real); 
    const In_imag = -(Ia_imag + Ib_imag + Ic_imag);

    // Calcular la magnitud de la corriente del neutro
    const In_magnitud = Math.sqrt(In_real ** 2 + In_imag ** 2);