Subversion Repositories pentevo

Rev

Rev 1251 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed

  1. #from sympy import *
  2. import math
  3.  
  4.  
  5. r1 = 1.2e3
  6. r2 = 47e3
  7. r3 = 68e3
  8. c3 = 10e-12
  9.  
  10.  
  11. def R4C1C2(w1,w2,q,r1,r2,r3,c3):
  12.  
  13.     A=1.0/w1 + 1.0/(w2*q)
  14.     B=1.0/(q*w1*w2) + 1.0/(w2*w2)
  15.     C=1.0/(w1*w2*w2)
  16.  
  17.     r4=(-(-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))/(3*(math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)) - (math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)/3 - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))/(3*c3*r1*(r2 + r3)*(r1 + r2 + r3))
  18.  
  19.     c1=-c3*(r1 + r2 + r3)*((-(-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))/(3*(math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)) - (math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)/3 - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))/(3*c3*r1*(r2 + r3)*(r1 + r2 + r3)))/(r1*r2) + (A - c3*r3)*(r1 + r2)/(r1*r2)
  20.  
  21.     c2=C*(r1 + r2)/(c3*r3*(-c3*(r1 + r2 + r3)*((-(-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))/(3*(math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)) - (math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)/3 - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))/(3*c3*r1*(r2 + r3)*(r1 + r2 + r3)))**2 + (A - c3*r3)*(r1 + r2)*((-(-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))/(3*(math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)) - (math.sqrt(-4*((-3*B*r1*r2*(r1 + r2)*(r1 + r2 + r3) + 6*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) - 3*r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**2/(c3**2*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**2))**3 + ((-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + 2*(c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**2)/2 + (-27*B*r1*r2*(A - c3*r3)*(r1 + r2)**2/c3 + 27*C*r1*r2*(r1 + r2)**2/c3 + 27*r1*r2*r3*(A - c3*r3)**2*(r1 + r2)**2)/(2*c3**2*r1*(r2 + r3)*(r1 + r2 + r3)**2) - (9*c3*r1*r2*r3*(r1 + r2 + r3) - 18*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))*(B*r1*r2*(r1 + r2)*(r1 + r2 + r3) - 2*c3*r1*r2*r3*(A - c3*r3)*(r1 + r2)*(r1 + r2 + r3) + r1*(A - c3*r3)**2*(r1 + r2)**2*(r2 + r3))/(2*c3**3*r1**2*(r2 + r3)**2*(r1 + r2 + r3)**3) + (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))**3/(c3**3*r1**3*(r2 + r3)**3*(r1 + r2 + r3)**3))**(1/3)/3 - (c3*r1*r2*r3*(r1 + r2 + r3) - 2*r1*(A - c3*r3)*(r1 + r2)*(r2 + r3))/(3*c3*r1*(r2 + r3)*(r1 + r2 + r3)))))
  22.  
  23.     return (r4,c1,c2)
  24.  
  25. f1=30000
  26. f2=f1
  27. q=1.0
  28.  
  29. w1=f1*2*math.pi
  30. w2=f2*2*math.pi
  31.  
  32. r4,c1,c2=R4C1C2(w1,w2,q,r1,r2,r3,c3)
  33.  
  34. print('\nfrom r1={},r2={},r3={},c3={}; f1={},f2={},q={}'.format(r1,r2,r3,c3,f1,f2,q))
  35. print('r4={},c1={},c2={}'.format(r4,c1,c2))
  36.  
  37.