r**2*((2*d**2 - 2)*U_cilia/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 4*Qm/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)) + ((-4*d**4 + 4*d**2)*U_cilia/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 8*Qm*d**2/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1))*np.log(r) + (4*d**4*np.log(d) - 3*d**4 + 2*d**2 + 1)*U_cilia/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 4*Qm/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)
