import itertools
from math import sqrt

# Capacitores
E6_C = [1.0, 1.5, 2.2, 3.3, 4.7, 6.8]

# resistencias de 5% de error
E24_R = [
    1.0, 1.1, 1.2, 1.3, 1.5, 1.6, 1.8, 2.0, 2.2, 2.4, 2.7, 3.0,
    3.3, 3.6, 3.9, 4.3, 4.7, 5.1, 5.6, 6.2, 6.8, 7.5, 8.2, 9.1
]

E24_R_custom = [
    1.0, 1.1, 1.2, 1.3, 1.5, 1.6, 1.8, 2.0, 2.2, 2.4, 2.7, 3.0,
    3.3, 3.6, 3.9, 4.3, 4.7, 5.1, 5.6, 6.2, 6.8, 7.5, 8.2
]

# resistencias de 2% de error
E48_R = [
    1.00, 1.05, 1.10, 1.15, 1.21, 1.27, 1.33, 1.40, 1.47, 1.54, 1.62, 1.69,
    1.78, 1.87, 1.96, 2.05, 2.15, 2.26, 2.37, 2.49, 2.61, 2.74, 2.87, 3.01,
    3.16, 3.32, 3.48, 3.65, 3.83, 4.02, 4.22, 4.42, 4.64, 4.87, 5.11, 5.36,
    5.62, 5.90, 6.19, 6.49, 6.81, 7.15, 7.50, 7.87, 8.25, 8.66, 9.09, 9.53
]

# resistencias de 1% de error
E96_R = [
    1.00, 1.02, 1.05, 1.07, 1.10, 1.13, 1.15, 1.18, 1.21, 1.24, 1.27, 1.30,
    1.33, 1.37, 1.40, 1.43, 1.47, 1.50, 1.54, 1.58, 1.62, 1.65, 1.69, 1.74,
    1.78, 1.82, 1.87, 1.91, 1.96, 2.00, 2.05, 2.10, 2.15, 2.21, 2.26, 2.32,
    2.37, 2.43, 2.49, 2.55, 2.61, 2.67, 2.74, 2.80, 2.87, 2.94, 3.01, 3.09,
    3.16, 3.24, 3.32, 3.40, 3.48, 3.57, 3.65, 3.74, 3.83, 3.92, 4.02, 4.12,
    4.22, 4.32, 4.42, 4.53, 4.64, 4.75, 4.87, 4.99, 5.11, 5.23, 5.36, 5.49,
    5.62, 5.76, 5.90, 6.04, 6.19, 6.34, 6.49, 6.65, 6.81, 6.98, 7.15, 7.32,
    7.50, 7.68, 7.87, 8.06, 8.25, 8.45, 8.66, 8.87, 9.09, 9.31, 9.53, 9.76
]


OMEGA0_OBJ = 4000.0
Q_OBJ = 2.25

SERIE_R = E24_R_custom
SERIE_C = E6_C

DECADAS_R = [1e3, 1e4, 1e5]
DECADAS_C = [1e-9, 10e-9, 100e-9, 1e-6]

resistencias = [v*d for v, d in itertools.product(SERIE_R, DECADAS_R)]
capacitores  = [v*d for v, d in itertools.product(SERIE_C, DECADAS_C)]

def calc_params(R1, R2, C1, C2):
    omega0 = sqrt(1/(R1*R2*C1*C2))
    Q = omega0 / ((R1 + R2) / (C1*R1*R2))

    return omega0, Q

def cost(omega, Q):
    err_w = (omega - OMEGA0_OBJ) / OMEGA0_OBJ
    err_q = (Q - Q_OBJ) / Q_OBJ

    return err_w**2 + err_q**2

resultados = []

for R1, R2, C1, C2 in itertools.product(
        resistencias,
        resistencias,
        capacitores,
        capacitores):

    omega, Q = calc_params(R1, R2, C1, C2)

    E = cost(omega,Q)

    resultados.append((
        E,
        omega,
        Q,
        R1,
        R2,
        C1,
        C2
    ))

resultados.sort(key=lambda x: x[0])

print(f"  R1 [Ω]  R2 [Ω]    C1 [F]    C2 [F]      w0  dw0 [%]       Q  dQ [%]       cost")
for E, omega, Q, R1, R2, C1, C2 in resultados[:20]:

    errw = 100 * (omega - OMEGA0_OBJ)/OMEGA0_OBJ
    errq = 100 * (Q - Q_OBJ)/Q_OBJ

    print(
        f"{R1:8.0f}  "
        f"{R2:6.0f}  "
        f"{C1:8.2e}  "
        f"{C2:8.2e}  "
        f"{omega:6.2f}  "
        f"{errw:+7.2f}  "
        f"{Q:6.3f}  "
        f"{errq:+6.2f}  "
        f"{E:.3e}"
    )