import numpy as np
def f(w_a, w_b1, w_b2, w_c11, w_c12, w_c21, w_c22, z):
ret = (z**2/(1.e-6 - 2*w_a)**2 + (2*z**2)/(1.e-6 - 2*w_b1)**2 +
(2*z**2)/(1.e-6 + 2*w_a - 2*w_b1)**2 + (1.e-6 + w_a - w_b1)**(-2) +
z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_a + w_b1)**(-2) +
z**2/(1.e-6 - w_a + w_b1)**2 + (2*z**2)/(1.e-6 - 2*w_a + 2*w_b1)**2 +
(2*z**2)/(1.e-6 - 2*w_b2)**2 + (2*z**2)/(1.e-6 + 2*w_a - 2*w_b2)**2 +
z**2/(1.e-6 + 2*w_a - 2*w_b1 - 2*w_b2)**2 + (1.e-6 + w_a - w_b2)**(-2) +
z**2/(1.e-6 + w_a - w_b2)**2 + z**2/(1.e-6 + 2*w_a - w_b1 - w_b2)**2 +
(1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_a + w_b2)**(-2) +
z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 - w_b1 + w_b2)**(-2) +
z**2/(1.e-6 - 2*w_a + w_b1 + w_b2)**2 + (2*z**2)/(1.e-6 - 2*w_a + 2*w_b2)**2 +
((1.e-6 + w_a - w_b1)**(-2) + (1.e-6 + w_b1 - w_c11)**(-2))/
(1.e-6 + 2*(w_a - w_c11))**2 +
((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 + w_b1 - w_c11)**(-2))/
(1.e-6 + 2*(w_b2 - w_c11))**2 +
(z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**2 + z**2/(1.e-6 - 2*w_c11)**2 +
(2*z**2)/(1.e-6 + 2*w_b1 - 2*w_c11)**2 + (1.e-6 + w_b1 - w_c11)**(-2) +
z**2/(1.e-6 + w_b1 - w_c11)**2 + (1.e-6 - w_b1 + w_c11)**(-2) +
z**2/(1.e-6 - w_b1 + w_c11)**2 + (2*z**2)/(1.e-6 - 2*w_b1 + 2*w_c11)**2 +
((1.e-6 - w_a + w_b1)**(-2) + (1.e-6 - w_b1 + w_c11)**(-2))/
(1.e-6 + 2*(-w_a + w_c11))**2 +
((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_b1 + w_c11)**(-2))/
(1.e-6 + 2*(-w_b2 + w_c11))**2 +
(z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**2 +
((1.e-6 + w_a - w_b1)**(-2) + (1.e-6 + w_b1 - w_c12)**(-2))/
(1.e-6 + 2*(w_a - w_c12))**2 +
((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 + w_b1 - w_c12)**(-2))/
(1.e-6 + 2*(w_b2 - w_c12))**2 +
(z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**2 +
(z**2*((1.e-6 + w_a - w_b1)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**2 +
(z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**2 +
(z**4*(z**2/(1.e-6 - w_a + w_b2)**2 + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**2 +
z**2/(1.e-6 - 2*w_c12)**2 + (2*z**2)/(1.e-6 + 2*w_b1 - 2*w_c12)**2 +
z**2/(1.e-6 + 2*w_b1 - 2*w_c11 - 2*w_c12)**2 + (1.e-6 + w_b1 - w_c12)**(-2) +
z**2/(1.e-6 + w_b1 - w_c12)**2 + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2 +
(1.e-6 + w_c11 - w_c12)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2) +
z**2/(1.e-6 - w_b1 + w_c12)**2 + (1.e-6 - w_c11 + w_c12)**(-2) +
z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2 +
(2*z**2)/(1.e-6 - 2*w_b1 + 2*w_c12)**2 +
((1.e-6 - w_a + w_b1)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2))/
(1.e-6 + 2*(-w_a + w_c12))**2 +
((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2))/
(1.e-6 + 2*(-w_b2 + w_c12))**2 +
(z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**2 +
(z**2*((1.e-6 - w_a + w_b1)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**2 +
(z**4*(z**2/(1.e-6 + w_a - w_b2)**2 + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**2 +
(z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**2 +
((1.e-6 + w_a - w_b2)**(-2) + (1.e-6 + w_b2 - w_c21)**(-2))/
(1.e-6 + 2*(w_a - w_c21))**2 +
((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 + w_b2 - w_c21)**(-2))/
(1.e-6 + 2*(w_b1 - w_c21))**2 +
(z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**2 + z**2/(1.e-6 - 2*w_c21)**2 +
(2*z**2)/(1.e-6 + 2*w_b2 - 2*w_c21)**2 + (1.e-6 + w_b2 - w_c21)**(-2) +
z**2/(1.e-6 + w_b2 - w_c21)**2 + (1.e-6 - w_b2 + w_c21)**(-2) +
z**2/(1.e-6 - w_b2 + w_c21)**2 + (2*z**2)/(1.e-6 - 2*w_b2 + 2*w_c21)**2 +
((1.e-6 - w_a + w_b2)**(-2) + (1.e-6 - w_b2 + w_c21)**(-2))/
(1.e-6 + 2*(-w_a + w_c21))**2 +
((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 - w_b2 + w_c21)**(-2))/
(1.e-6 + 2*(-w_b1 + w_c21))**2 +
(z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**2 +
((1.e-6 + w_a - w_b2)**(-2) + (1.e-6 + w_b2 - w_c22)**(-2))/
(1.e-6 + 2*(w_a - w_c22))**2 +
((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 + w_b2 - w_c22)**(-2))/
(1.e-6 + 2*(w_b1 - w_c22))**2 +
(z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**2 +
(z**2*((1.e-6 + w_a - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**2 +
(z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**2 +
(z**4*(z**2/(1.e-6 - w_a + w_b1)**2 + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**2 +
z**2/(1.e-6 - 2*w_c22)**2 + (2*z**2)/(1.e-6 + 2*w_b2 - 2*w_c22)**2 +
z**2/(1.e-6 + 2*w_b2 - 2*w_c21 - 2*w_c22)**2 + (1.e-6 + w_b2 - w_c22)**(-2) +
z**2/(1.e-6 + w_b2 - w_c22)**2 + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2 +
(1.e-6 + w_c21 - w_c22)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2) +
z**2/(1.e-6 - w_b2 + w_c22)**2 + (1.e-6 - w_c21 + w_c22)**(-2) +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2 +
(2*z**2)/(1.e-6 - 2*w_b2 + 2*w_c22)**2 +
((1.e-6 - w_a + w_b2)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2))/
(1.e-6 + 2*(-w_a + w_c22))**2 +
((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2))/
(1.e-6 + 2*(-w_b1 + w_c22))**2 +
(z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**2 +
(z**4*(z**2/(1.e-6 + w_a - w_b1)**2 + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**2 +
(z**2*((1.e-6 - w_a + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**2 +
(z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**2)
return ret
def df(w_a, w_b1, w_b2, w_c11, w_c12, w_c21, w_c22, z):
ret = ((4*z**2)/(1.e-6 - 2*w_a)**3 - (8*z**2)/(1.e-6 + 2*w_a - 2*w_b1)**3 -
2/(1.e-6 + w_a - w_b1)**3 - (2*z**2)/(1.e-6 + w_a - w_b1)**3 +
2/(1.e-6 - w_a + w_b1)**3 + (2*z**2)/(1.e-6 - w_a + w_b1)**3 +
(8*z**2)/(1.e-6 - 2*w_a + 2*w_b1)**3 - (8*z**2)/(1.e-6 + 2*w_a - 2*w_b2)**3 -
(4*z**2)/(1.e-6 + 2*w_a - 2*w_b1 - 2*w_b2)**3 - 2/(1.e-6 + w_a - w_b2)**3 -
(2*z**2)/(1.e-6 + w_a - w_b2)**3 - (4*z**2)/(1.e-6 + 2*w_a - w_b1 - w_b2)**3 +
2/(1.e-6 - w_a + w_b2)**3 + (2*z**2)/(1.e-6 - w_a + w_b2)**3 +
(4*z**2)/(1.e-6 - 2*w_a + w_b1 + w_b2)**3 +
(8*z**2)/(1.e-6 - 2*w_a + 2*w_b2)**3 -
2/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_c11))**2) -
(4*((1.e-6 + w_a - w_b1)**(-2) + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_c11))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**3 +
(2*z**4)/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**2) +
2/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_c11))**2) +
(4*((1.e-6 - w_a + w_b1)**(-2) + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_c11))**3 -
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**3 -
(2*z**4)/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**2) -
2/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_c12))**2) -
(4*((1.e-6 + w_a - w_b1)**(-2) + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_c12))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**3 +
(2*z**4)/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**2) -
(2*z**2)/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**2) -
(4*z**2*((1.e-6 + w_a - w_b1)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**3 +
(4*z**4*(z**2/(1.e-6 - w_a + w_b2)**2 +
z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**3 +
(2*z**6)/
((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**2) +
2/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_c12))**2) +
(4*((1.e-6 - w_a + w_b1)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_c12))**3 -
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**3 -
(2*z**4)/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**2) +
(4*z**2*((1.e-6 - w_a + w_b1)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**3 +
(2*z**2)/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**2) -
(4*z**4*(z**2/(1.e-6 + w_a - w_b2)**2 +
z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**3 -
(2*z**6)/
((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**2) -
2/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - w_c21))**2) -
(4*((1.e-6 + w_a - w_b2)**(-2) + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_c21))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**3 +
(2*z**4)/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**2) +
2/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + w_c21))**2) +
(4*((1.e-6 - w_a + w_b2)**(-2) + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_c21))**3 -
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**3 -
(2*z**4)/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**2) -
2/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - w_c22))**2) -
(4*((1.e-6 + w_a - w_b2)**(-2) + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_c22))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**3 +
(2*z**4)/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**2) -
(2*z**2)/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**2) -
(4*z**2*((1.e-6 + w_a - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**3 +
(4*z**4*(z**2/(1.e-6 - w_a + w_b1)**2 +
z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**3 +
(2*z**6)/
((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**2) +
2/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + w_c22))**2) +
(4*((1.e-6 - w_a + w_b2)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_c22))**3 -
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**3 -
(2*z**4)/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**2) -
(4*z**4*(z**2/(1.e-6 + w_a - w_b1)**2 +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**3 -
(2*z**6)/
((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**2) +
(4*z**2*((1.e-6 - w_a + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**3 +
(2*z**2)/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**2),
(8*z**2)/(1.e-6 - 2*w_b1)**3 + (8*z**2)/(1.e-6 + 2*w_a - 2*w_b1)**3 +
2/(1.e-6 + w_a - w_b1)**3 + (2*z**2)/(1.e-6 + w_a - w_b1)**3 -
2/(1.e-6 - w_a + w_b1)**3 - (2*z**2)/(1.e-6 - w_a + w_b1)**3 -
(8*z**2)/(1.e-6 - 2*w_a + 2*w_b1)**3 +
(4*z**2)/(1.e-6 + 2*w_a - 2*w_b1 - 2*w_b2)**3 +
(2*z**2)/(1.e-6 + 2*w_a - w_b1 - w_b2)**3 - 2/(1.e-6 + w_b1 - w_b2)**3 +
2/(1.e-6 - w_b1 + w_b2)**3 - (2*z**2)/(1.e-6 - 2*w_a + w_b1 + w_b2)**3 +
(2/(1.e-6 + w_a - w_b1)**3 - 2/(1.e-6 + w_b1 - w_c11)**3)/
(1.e-6 + 2*(w_a - w_c11))**2 +
(2/(1.e-6 - w_b1 + w_b2)**3 - 2/(1.e-6 + w_b1 - w_c11)**3)/
(1.e-6 + 2*(w_b2 - w_c11))**2 -
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**3 -
(8*z**2)/(1.e-6 + 2*w_b1 - 2*w_c11)**3 - 2/(1.e-6 + w_b1 - w_c11)**3 -
(2*z**2)/(1.e-6 + w_b1 - w_c11)**3 -
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**2*(1.e-6 + w_b1 - w_c11)**3) +
2/(1.e-6 - w_b1 + w_c11)**3 + (2*z**2)/(1.e-6 - w_b1 + w_c11)**3 +
(8*z**2)/(1.e-6 - 2*w_b1 + 2*w_c11)**3 +
(-2/(1.e-6 - w_a + w_b1)**3 + 2/(1.e-6 - w_b1 + w_c11)**3)/
(1.e-6 + 2*(-w_a + w_c11))**2 +
(-2/(1.e-6 + w_b1 - w_b2)**3 + 2/(1.e-6 - w_b1 + w_c11)**3)/
(1.e-6 + 2*(-w_b2 + w_c11))**2 +
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**3 +
(2*z**2)/((1.e-6 - w_b1 + w_c11)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**2) +
(2/(1.e-6 + w_a - w_b1)**3 - 2/(1.e-6 + w_b1 - w_c12)**3)/
(1.e-6 + 2*(w_a - w_c12))**2 +
(2/(1.e-6 - w_b1 + w_b2)**3 - 2/(1.e-6 + w_b1 - w_c12)**3)/
(1.e-6 + 2*(w_b2 - w_c12))**2 -
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**3 +
(z**2*(2/(1.e-6 + w_a - w_b1)**3 - (4*z**2)/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3))/
(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**2 -
(4*z**2*((1.e-6 + w_a - w_b1)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**3 -
(4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**3 +
(z**2*(2/(1.e-6 - w_b1 + w_b2)**3 -
(4*z**2)/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**2 -
(8*z**4*(z**2/(1.e-6 - w_a + w_b2)**2 +
z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**3 -
(8*z**2)/(1.e-6 + 2*w_b1 - 2*w_c12)**3 -
(4*z**2)/(1.e-6 + 2*w_b1 - 2*w_c11 - 2*w_c12)**3 - 2/(1.e-6 + w_b1 - w_c12)**3 -
(2*z**2)/(1.e-6 + w_b1 - w_c12)**3 -
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**2*(1.e-6 + w_b1 - w_c12)**3) -
(4*z**2)/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3 -
(4*z**6)/
((1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**2*
(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3) + 2/(1.e-6 - w_b1 + w_c12)**3 +
(2*z**2)/(1.e-6 - w_b1 + w_c12)**3 +
(4*z**2)/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**3 +
(8*z**2)/(1.e-6 - 2*w_b1 + 2*w_c12)**3 +
(-2/(1.e-6 - w_a + w_b1)**3 + 2/(1.e-6 - w_b1 + w_c12)**3)/
(1.e-6 + 2*(-w_a + w_c12))**2 +
(-2/(1.e-6 + w_b1 - w_b2)**3 + 2/(1.e-6 - w_b1 + w_c12)**3)/
(1.e-6 + 2*(-w_b2 + w_c12))**2 +
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**3 +
(2*z**2)/((1.e-6 - w_b1 + w_c12)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**2) +
(4*z**2*((1.e-6 - w_a + w_b1)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**3 +
(z**2*(-2/(1.e-6 - w_a + w_b1)**3 +
(4*z**2)/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**3))/
(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**2 +
(8*z**4*(z**2/(1.e-6 + w_a - w_b2)**2 +
z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**3 +
(4*z**6)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**2) +
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**3 +
(z**2*(-2/(1.e-6 + w_b1 - w_b2)**3 +
(4*z**2)/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**3))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**2 -
2/((1.e-6 + w_b1 - w_b2)**3*(1.e-6 + 2*(w_b1 - w_c21))**2) -
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(w_b1 - w_c21))**3 -
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**3 -
(2*z**4)/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**2) +
2/((1.e-6 - w_b1 + w_b2)**3*(1.e-6 + 2*(-w_b1 + w_c21))**2) +
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(-w_b1 + w_c21))**3 +
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**3 +
(2*z**4)/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**2) -
2/((1.e-6 + w_b1 - w_b2)**3*(1.e-6 + 2*(w_b1 - w_c22))**2) -
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(w_b1 - w_c22))**3 -
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**3 -
(2*z**4)/((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**2) -
(2*z**2)/((1.e-6 + w_b1 - w_b2)**3*(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**2) -
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**3 -
(4*z**4*(z**2/(1.e-6 - w_a + w_b1)**2 +
z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**3 -
(2*z**6)/
((1.e-6 - w_a + w_b1)**3*(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**2) +
2/((1.e-6 - w_b1 + w_b2)**3*(1.e-6 + 2*(-w_b1 + w_c22))**2) +
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(-w_b1 + w_c22))**3 +
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**3 +
(2*z**4)/((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**2) +
(4*z**4*(z**2/(1.e-6 + w_a - w_b1)**2 +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**3 +
(2*z**6)/
((1.e-6 + w_a - w_b1)**3*(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**2) +
(4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**3 +
(2*z**2)/((1.e-6 - w_b1 + w_b2)**3*(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**2),
(8*z**2)/(1.e-6 - 2*w_b2)**3 + (8*z**2)/(1.e-6 + 2*w_a - 2*w_b2)**3 +
(4*z**2)/(1.e-6 + 2*w_a - 2*w_b1 - 2*w_b2)**3 + 2/(1.e-6 + w_a - w_b2)**3 +
(2*z**2)/(1.e-6 + w_a - w_b2)**3 + (2*z**2)/(1.e-6 + 2*w_a - w_b1 - w_b2)**3 +
2/(1.e-6 + w_b1 - w_b2)**3 - 2/(1.e-6 - w_a + w_b2)**3 -
(2*z**2)/(1.e-6 - w_a + w_b2)**3 - 2/(1.e-6 - w_b1 + w_b2)**3 -
(2*z**2)/(1.e-6 - 2*w_a + w_b1 + w_b2)**3 -
(8*z**2)/(1.e-6 - 2*w_a + 2*w_b2)**3 -
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(w_b2 - w_c11))**3 -
2/((1.e-6 - w_b1 + w_b2)**3*(1.e-6 + 2*(w_b2 - w_c11))**2) -
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**3 -
(2*z**4)/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**2) +
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(-w_b2 + w_c11))**3 +
2/((1.e-6 + w_b1 - w_b2)**3*(1.e-6 + 2*(-w_b2 + w_c11))**2) +
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**3 +
(2*z**4)/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**2) -
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(w_b2 - w_c12))**3 -
2/((1.e-6 - w_b1 + w_b2)**3*(1.e-6 + 2*(w_b2 - w_c12))**2) -
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**3 -
(2*z**4)/((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**2) -
(4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**3 -
(2*z**2)/((1.e-6 - w_b1 + w_b2)**3*(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**2) -
(4*z**4*(z**2/(1.e-6 - w_a + w_b2)**2 +
z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**3 -
(2*z**6)/
((1.e-6 - w_a + w_b2)**3*(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**2) +
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(-w_b2 + w_c12))**3 +
2/((1.e-6 + w_b1 - w_b2)**3*(1.e-6 + 2*(-w_b2 + w_c12))**2) +
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**3 +
(2*z**4)/((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**2) +
(4*z**4*(z**2/(1.e-6 + w_a - w_b2)**2 +
z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**3 +
(2*z**6)/
((1.e-6 + w_a - w_b2)**3*(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**2) +
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**3 +
(2*z**2)/
((1.e-6 + w_b1 - w_b2)**3*(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**2) +
(2/(1.e-6 + w_a - w_b2)**3 - 2/(1.e-6 + w_b2 - w_c21)**3)/
(1.e-6 + 2*(w_a - w_c21))**2 +
(2/(1.e-6 + w_b1 - w_b2)**3 - 2/(1.e-6 + w_b2 - w_c21)**3)/
(1.e-6 + 2*(w_b1 - w_c21))**2 -
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**3 -
(8*z**2)/(1.e-6 + 2*w_b2 - 2*w_c21)**3 - 2/(1.e-6 + w_b2 - w_c21)**3 -
(2*z**2)/(1.e-6 + w_b2 - w_c21)**3 -
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**2*(1.e-6 + w_b2 - w_c21)**3) +
2/(1.e-6 - w_b2 + w_c21)**3 + (2*z**2)/(1.e-6 - w_b2 + w_c21)**3 +
(8*z**2)/(1.e-6 - 2*w_b2 + 2*w_c21)**3 +
(-2/(1.e-6 - w_a + w_b2)**3 + 2/(1.e-6 - w_b2 + w_c21)**3)/
(1.e-6 + 2*(-w_a + w_c21))**2 +
(-2/(1.e-6 - w_b1 + w_b2)**3 + 2/(1.e-6 - w_b2 + w_c21)**3)/
(1.e-6 + 2*(-w_b1 + w_c21))**2 +
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**3 +
(2*z**2)/((1.e-6 - w_b2 + w_c21)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**2) +
(2/(1.e-6 + w_a - w_b2)**3 - 2/(1.e-6 + w_b2 - w_c22)**3)/
(1.e-6 + 2*(w_a - w_c22))**2 +
(2/(1.e-6 + w_b1 - w_b2)**3 - 2/(1.e-6 + w_b2 - w_c22)**3)/
(1.e-6 + 2*(w_b1 - w_c22))**2 -
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**3 +
(z**2*(2/(1.e-6 + w_a - w_b2)**3 - (4*z**2)/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3))/
(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**2 +
(z**2*(2/(1.e-6 + w_b1 - w_b2)**3 -
(4*z**2)/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**2 -
(4*z**2*((1.e-6 + w_a - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**3 -
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**3 -
(8*z**4*(z**2/(1.e-6 - w_a + w_b1)**2 +
z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**3 -
(8*z**2)/(1.e-6 + 2*w_b2 - 2*w_c22)**3 -
(4*z**2)/(1.e-6 + 2*w_b2 - 2*w_c21 - 2*w_c22)**3 - 2/(1.e-6 + w_b2 - w_c22)**3 -
(2*z**2)/(1.e-6 + w_b2 - w_c22)**3 -
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**2*(1.e-6 + w_b2 - w_c22)**3) -
(4*z**2)/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3 -
(4*z**6)/
((1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**2*
(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3) + 2/(1.e-6 - w_b2 + w_c22)**3 +
(2*z**2)/(1.e-6 - w_b2 + w_c22)**3 +
(4*z**2)/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**3 +
(8*z**2)/(1.e-6 - 2*w_b2 + 2*w_c22)**3 +
(-2/(1.e-6 - w_a + w_b2)**3 + 2/(1.e-6 - w_b2 + w_c22)**3)/
(1.e-6 + 2*(-w_a + w_c22))**2 +
(-2/(1.e-6 - w_b1 + w_b2)**3 + 2/(1.e-6 - w_b2 + w_c22)**3)/
(1.e-6 + 2*(-w_b1 + w_c22))**2 +
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**3 +
(2*z**2)/((1.e-6 - w_b2 + w_c22)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**2) +
(8*z**4*(z**2/(1.e-6 + w_a - w_b1)**2 +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**3 +
(4*z**6)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**2) +
(4*z**2*((1.e-6 - w_a + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**3 +
(z**2*(-2/(1.e-6 - w_a + w_b2)**3 +
(4*z**2)/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**3))/
(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**2 +
(4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**3 +
(z**2*(-2/(1.e-6 - w_b1 + w_b2)**3 +
(4*z**2)/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**3))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**2,
(4*((1.e-6 + w_a - w_b1)**(-2) + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_c11))**3 +
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(w_b2 - w_c11))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**3 + (4*z**2)/(1.e-6 - 2*w_c11)**3 +
(8*z**2)/(1.e-6 + 2*w_b1 - 2*w_c11)**3 + 2/(1.e-6 + w_b1 - w_c11)**3 +
(2*z**2)/(1.e-6 + w_b1 - w_c11)**3 +
2/((1.e-6 + 2*(w_a - w_c11))**2*(1.e-6 + w_b1 - w_c11)**3) +
2/((1.e-6 + 2*(w_b2 - w_c11))**2*(1.e-6 + w_b1 - w_c11)**3) +
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c11))**2*(1.e-6 + w_b1 - w_c11)**3) -
2/(1.e-6 - w_b1 + w_c11)**3 - (2*z**2)/(1.e-6 - w_b1 + w_c11)**3 -
(8*z**2)/(1.e-6 - 2*w_b1 + 2*w_c11)**3 -
2/((1.e-6 - w_b1 + w_c11)**3*(1.e-6 + 2*(-w_a + w_c11))**2) -
(4*((1.e-6 - w_a + w_b1)**(-2) + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(-w_a + w_c11))**3 -
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(-w_b2 + w_c11))**3 -
2/((1.e-6 - w_b1 + w_c11)**3*(1.e-6 + 2*(-w_b2 + w_c11))**2) -
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c11)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**3 -
(2*z**2)/((1.e-6 - w_b1 + w_c11)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c11))**2) +
(4*z**2*((1.e-6 + w_a - w_b1)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**3 +
(4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**3 +
(4*z**4*(z**2/(1.e-6 - w_a + w_b2)**2 +
z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**3 +
(4*z**2)/(1.e-6 + 2*w_b1 - 2*w_c11 - 2*w_c12)**3 +
(2*z**2)/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3 +
(2*z**4)/
((1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**2*(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3)\
+ (2*z**4)/
((1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**2*(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3)\
+ (2*z**6)/
((1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**2*
(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3) - 2/(1.e-6 + w_c11 - w_c12)**3 +
2/(1.e-6 - w_c11 + w_c12)**3 - (2*z**2)/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**3 -
(4*z**2*((1.e-6 - w_a + w_b1)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**3 -
(2*z**4)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**2)\
- (4*z**4*(z**2/(1.e-6 + w_a - w_b2)**2 +
z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**3 -
(2*z**6)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**2) -
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**3 -
(2*z**4)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**2)
,(4*((1.e-6 + w_a - w_b1)**(-2) + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_c12))**3 +
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(w_b2 - w_c12))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b2)**2 + (1.e-6 + w_b1 - w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**3 +
(4*z**2*((1.e-6 + w_a - w_b1)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**3 +
(4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) + z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**3 +
(4*z**4*(z**2/(1.e-6 - w_a + w_b2)**2 +
z**2/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**2))/
(1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**3 +
(4*z**2)/(1.e-6 - 2*w_c12)**3 + (8*z**2)/(1.e-6 + 2*w_b1 - 2*w_c12)**3 +
(4*z**2)/(1.e-6 + 2*w_b1 - 2*w_c11 - 2*w_c12)**3 + 2/(1.e-6 + w_b1 - w_c12)**3 +
(2*z**2)/(1.e-6 + w_b1 - w_c12)**3 +
2/((1.e-6 + 2*(w_a - w_c12))**2*(1.e-6 + w_b1 - w_c12)**3) +
2/((1.e-6 + 2*(w_b2 - w_c12))**2*(1.e-6 + w_b1 - w_c12)**3) +
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c12))**2*(1.e-6 + w_b1 - w_c12)**3) +
(2*z**2)/(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3 +
(2*z**4)/
((1.e-6 + 2*(w_a + w_b1 - w_c11 - w_c12))**2*(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3)\
+ (2*z**4)/
((1.e-6 + 2*(w_b1 + w_b2 - w_c11 - w_c12))**2*(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3)\
+ (2*z**6)/
((1.e-6 + 2*(-w_a + 2*w_b1 + w_b2 - w_c11 - w_c12))**2*
(1.e-6 + 2*w_b1 - w_c11 - w_c12)**3) + 2/(1.e-6 + w_c11 - w_c12)**3 -
2/(1.e-6 - w_b1 + w_c12)**3 - (2*z**2)/(1.e-6 - w_b1 + w_c12)**3 -
2/(1.e-6 - w_c11 + w_c12)**3 - (2*z**2)/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**3 -
(8*z**2)/(1.e-6 - 2*w_b1 + 2*w_c12)**3 -
2/((1.e-6 - w_b1 + w_c12)**3*(1.e-6 + 2*(-w_a + w_c12))**2) -
(4*((1.e-6 - w_a + w_b1)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(-w_a + w_c12))**3 -
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(-w_b2 + w_c12))**3 -
2/((1.e-6 - w_b1 + w_c12)**3*(1.e-6 + 2*(-w_b2 + w_c12))**2) -
(4*z**2*(z**2/(1.e-6 + w_a - w_b2)**2 + (1.e-6 - w_b1 + w_c12)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**3 -
(2*z**2)/((1.e-6 - w_b1 + w_c12)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c12))**2) -
(4*z**2*((1.e-6 - w_a + w_b1)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**3 -
(2*z**4)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*(1.e-6 + 2*(-w_a - w_b1 + w_c11 + w_c12))**2)\
- (4*z**4*(z**2/(1.e-6 + w_a - w_b2)**2 +
z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**3 -
(2*z**6)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*
(1.e-6 + 2*(w_a - 2*w_b1 - w_b2 + w_c11 + w_c12))**2) -
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 - 2*w_b1 + w_c11 + w_c12)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**3 -
(2*z**4)/
((1.e-6 - 2*w_b1 + w_c11 + w_c12)**3*(1.e-6 + 2*(-w_b1 - w_b2 + w_c11 + w_c12))**2)
,(4*((1.e-6 + w_a - w_b2)**(-2) + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_c21))**3 +
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(w_b1 - w_c21))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**3 + (4*z**2)/(1.e-6 - 2*w_c21)**3 +
(8*z**2)/(1.e-6 + 2*w_b2 - 2*w_c21)**3 + 2/(1.e-6 + w_b2 - w_c21)**3 +
(2*z**2)/(1.e-6 + w_b2 - w_c21)**3 +
2/((1.e-6 + 2*(w_a - w_c21))**2*(1.e-6 + w_b2 - w_c21)**3) +
2/((1.e-6 + 2*(w_b1 - w_c21))**2*(1.e-6 + w_b2 - w_c21)**3) +
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c21))**2*(1.e-6 + w_b2 - w_c21)**3) -
2/(1.e-6 - w_b2 + w_c21)**3 - (2*z**2)/(1.e-6 - w_b2 + w_c21)**3 -
(8*z**2)/(1.e-6 - 2*w_b2 + 2*w_c21)**3 -
2/((1.e-6 - w_b2 + w_c21)**3*(1.e-6 + 2*(-w_a + w_c21))**2) -
2/((1.e-6 - w_b2 + w_c21)**3*(1.e-6 + 2*(-w_b1 + w_c21))**2) -
(4*((1.e-6 - w_a + w_b2)**(-2) + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(-w_a + w_c21))**3 -
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(-w_b1 + w_c21))**3 -
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c21)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**3 -
(2*z**2)/((1.e-6 - w_b2 + w_c21)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c21))**2) +
(4*z**2*((1.e-6 + w_a - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**3 +
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**3 +
(4*z**4*(z**2/(1.e-6 - w_a + w_b1)**2 +
z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**3 +
(4*z**2)/(1.e-6 + 2*w_b2 - 2*w_c21 - 2*w_c22)**3 +
(2*z**2)/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3 +
(2*z**4)/
((1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**2*(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3)\
+ (2*z**4)/
((1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**2*(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3)\
+ (2*z**6)/
((1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**2*
(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3) - 2/(1.e-6 + w_c21 - w_c22)**3 +
2/(1.e-6 - w_c21 + w_c22)**3 - (2*z**2)/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**3 -
(4*z**4*(z**2/(1.e-6 + w_a - w_b1)**2 +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**3 -
(2*z**6)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**2) -
(4*z**2*((1.e-6 - w_a + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**3 -
(2*z**4)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**2)\
- (4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**3 -
(2*z**4)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**2)
,(4*((1.e-6 + w_a - w_b2)**(-2) + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_c22))**3 +
(4*((1.e-6 + w_b1 - w_b2)**(-2) + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(w_b1 - w_c22))**3 +
(4*z**2*(z**2/(1.e-6 - w_a + w_b1)**2 + (1.e-6 + w_b2 - w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**3 +
(4*z**2*((1.e-6 + w_a - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**3 +
(4*z**2*((1.e-6 + w_b1 - w_b2)**(-2) + z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**3 +
(4*z**4*(z**2/(1.e-6 - w_a + w_b1)**2 +
z**2/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**2))/
(1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**3 +
(4*z**2)/(1.e-6 - 2*w_c22)**3 + (8*z**2)/(1.e-6 + 2*w_b2 - 2*w_c22)**3 +
(4*z**2)/(1.e-6 + 2*w_b2 - 2*w_c21 - 2*w_c22)**3 + 2/(1.e-6 + w_b2 - w_c22)**3 +
(2*z**2)/(1.e-6 + w_b2 - w_c22)**3 +
2/((1.e-6 + 2*(w_a - w_c22))**2*(1.e-6 + w_b2 - w_c22)**3) +
2/((1.e-6 + 2*(w_b1 - w_c22))**2*(1.e-6 + w_b2 - w_c22)**3) +
(2*z**2)/((1.e-6 + 2*(-w_a + w_b1 + w_b2 - w_c22))**2*(1.e-6 + w_b2 - w_c22)**3) +
(2*z**2)/(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3 +
(2*z**4)/
((1.e-6 + 2*(w_a + w_b2 - w_c21 - w_c22))**2*(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3)\
+ (2*z**4)/
((1.e-6 + 2*(w_b1 + w_b2 - w_c21 - w_c22))**2*(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3)\
+ (2*z**6)/
((1.e-6 + 2*(-w_a + w_b1 + 2*w_b2 - w_c21 - w_c22))**2*
(1.e-6 + 2*w_b2 - w_c21 - w_c22)**3) + 2/(1.e-6 + w_c21 - w_c22)**3 -
2/(1.e-6 - w_b2 + w_c22)**3 - (2*z**2)/(1.e-6 - w_b2 + w_c22)**3 -
2/(1.e-6 - w_c21 + w_c22)**3 - (2*z**2)/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**3 -
(8*z**2)/(1.e-6 - 2*w_b2 + 2*w_c22)**3 -
2/((1.e-6 - w_b2 + w_c22)**3*(1.e-6 + 2*(-w_a + w_c22))**2) -
2/((1.e-6 - w_b2 + w_c22)**3*(1.e-6 + 2*(-w_b1 + w_c22))**2) -
(4*((1.e-6 - w_a + w_b2)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(-w_a + w_c22))**3 -
(4*((1.e-6 - w_b1 + w_b2)**(-2) + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(-w_b1 + w_c22))**3 -
(4*z**2*(z**2/(1.e-6 + w_a - w_b1)**2 + (1.e-6 - w_b2 + w_c22)**(-2)))/
(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**3 -
(2*z**2)/((1.e-6 - w_b2 + w_c22)**3*(1.e-6 + 2*(w_a - w_b1 - w_b2 + w_c22))**2) -
(4*z**4*(z**2/(1.e-6 + w_a - w_b1)**2 +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**3 -
(2*z**6)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*
(1.e-6 + 2*(w_a - w_b1 - 2*w_b2 + w_c21 + w_c22))**2) -
(4*z**2*((1.e-6 - w_a + w_b2)**(-2) + z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**3 -
(2*z**4)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*(1.e-6 + 2*(-w_a - w_b2 + w_c21 + w_c22))**2)\
- (4*z**2*((1.e-6 - w_b1 + w_b2)**(-2) +
z**2/(1.e-6 - 2*w_b2 + w_c21 + w_c22)**2))/
(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**3 -
(2*z**4)/
((1.e-6 - 2*w_b2 + w_c21 + w_c22)**3*(1.e-6 + 2*(-w_b1 - w_b2 + w_c21 + w_c22))**2))
return np.array(ret)