TB066

 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}"
     )