u

<pre>X = part1['X']
n = len(X)
a = np.zeros(n)
x = np.zeros(n)
for i in range(0,n):
    for j in range(25000,33333):
        da = 0.00001
        if (abs(x[i]-X[i])/X[i]<=0.01):
            a[i]=j*da
            break
            x[i] = (4*j*da-1)*((1-j*da)/(1-3*j*da))**0.5
def equations(vars):
    x, y, z = vars
    eq1 = (x/(1-x)) - (5/(8*np.pi*np.sin(z)**2)) * (1 + 4*(np.tan(z)))
    eq2 = (y/(1-y)) - (5/(8*np.pi * np.cos(z))) * (1 - (1/np.tan(z)))
    eq3 = np.tan(z) - ((1-x) / 3 * (1+y))
    return [eq1, eq2, eq3]

x, y, z =  fsolve(equations,(0.288,2,2))

print(x, y, z)