AC Power Calculator

import math

def ac_power(system, V, I, pf):
    """
    Calculate AC power components.

    system: 'single' or 'three'
    V: voltage (V) – line-to-line for three-phase
    I: current (A)
    pf: power factor (0 to 1)
    """

    if not (0 < pf <= 1):
        raise ValueError("Power factor must be between 0 and 1")

    phi = math.acos(pf)
    sin_phi = math.sin(phi)

    if system == 'three':
        S = math.sqrt(3) * V * I
    else:
        S = V * I

    P = S * pf
    Q = S * sin_phi

    return {
        "P_W": P,
        "Q_VAR": Q,
        "S_VA": S,
        "P_kW": P / 1000,
        "Q_kVAR": Q / 1000,
        "S_kVA": S / 1000,
        "pf": pf
    }


# Example
results = ac_power('three', 400, 50, 0.85)
for k, v in results.items():
    print(f"{k}: {v}")

function results = ac_power(system, V, I, pf)
    % AC power calculation
    % system: 'single' or 'three'
    % V = line-to-line voltage for three-phase

    if pf <= 0 || pf > 1
        error('Power factor must be between 0 and 1');
    end

    phi = acos(pf);
    sin_phi = sin(phi);

    if strcmp(system, 'three')
        S = sqrt(3) * V * I;
    else
        S = V * I;
    end

    results.P = S * pf;
    results.Q = S * sin_phi;
    results.S = S;
    results.P_kW = results.P / 1000;
    results.Q_kVAR = results.Q / 1000;
    results.S_kVA = results.S / 1000;
    results.pf = pf;
end

% Example
r = ac_power('three', 400, 50, 0.85);
disp(r);