128*Qm*d**5*np.log(d)**2*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 128*Qm*d**5*np.log(d)*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 32*Qm*d**5*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 128*Qm*d**3*np.log(d)*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 64*Qm*d**3*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 32*Qm*d*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 16*U_cilia*d**7*np.log(d)*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 8*U_cilia*d**7*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 64*U_cilia*d**5*np.log(d)**2*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 32*U_cilia*d**5*np.log(d)*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 8*U_cilia*d**5*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 48*U_cilia*d**3*np.log(d)*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) + 8*U_cilia*d**3*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 8*U_cilia*d*dddt/(16*d**8*np.log(d)**2 - 24*d**8*np.log(d) + 9*d**8 + 32*d**6*np.log(d) - 24*d**6 - 8*d**4*np.log(d) + 22*d**4 - 8*d**2 + 1) - 16*Qm*d*np.log(d)*dddt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 4*U_cilia*d**3*dddt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 8*U_cilia*d*np.log(d)*dddt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 4*U_cilia*d*dddt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - d**4*dU_cdt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 8*d**2*np.log(d)*dQmdt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 4*d**2*np.log(d)*dU_cdt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + 4*d**2*dQmdt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) - 4*dQmdt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1) + dU_cdt/(4*d**4*np.log(d) - 3*d**4 + 4*d**2 - 1)