-
Notifications
You must be signed in to change notification settings - Fork 0
/
TKD_residual.py
17 lines (16 loc) · 87.8 KB
/
TKD_residual.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
# -*- coding: utf-8 -*-
"""
"""
from __future__ import division
import numpy as np
from numpy import sqrt
one6 = 1./6.
one3 = 1./3.
two3 = 2./3.
four3 = 4./3.
def Equations(mu, fo, d, alpha, q1y, q2x, q3z, lam1, lam2, lam3, pressure):
R=np.zeros(3)
R[0] = 0.0222300324977063*alpha*mu*(fo - 1)**2*(-4*q1y - four3)*(q1y**2 + two3*q1y - q2x**2 - two3*q2x)**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**3) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) - 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) + 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) - 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q1y + one3)**2 + (q2x + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q1y - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) + 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q1y + one3)**2 + (q2x + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q1y + one3)*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*(q1y + one3)*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q1y + one3)*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q1y + one3)*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q1y + one3)*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q1y + one3)*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*(q1y + one3)*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(4*q1y + four3)*(q1y**2 + two3*q1y - q2x**2 - two3*q2x)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-two3*lam2 + 2.0*q1y + 0.5*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-one3*lam2 + q1y + 0.25*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + one3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-two3*lam2 + 2.0*q1y + 0.5*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-one3*lam2 + q1y + 0.25*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + one3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-four3*lam2 + 4*q1y + 1.0*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 1.0*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + four3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-two3*lam2 + 2.0*q1y + 0.5*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))*(-two3*lam2 + 2.0*q1y + 0.5*(q1y + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0889201299908253*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(q3z + one3) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2*(2.66666666666667*lam2 - 4*q1y - four3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**3) + 0.0889201299908253*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(q3z + one3) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))*(four3*lam2 - 2.0*q1y - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - 0.666666666666666)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - 0.666666666666666)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam2 - 2*q1y - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam2 + q1y + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(q3z + one3) - 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))*(four3*lam2 - 2.0*q1y - 1.66533453693773e-16*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - two3 - 1.66533453693773e-16*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam2 - 2*q1y - two3)*((0.5*q3z + one6)*(q3z + one3) - 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(q3z + one3) + 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))*(four3*lam2 - 2.0*q1y + 1.66533453693773e-16*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - two3 + 1.66533453693773e-16*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam2 + q1y + one3)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam2 - 2*q1y - 0.666666666666666)*((0.5*q3z + one6)*(q3z + one3) + 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam2 - 2*q1y - two3)*((0.5*q3z + one6)*(q3z + one3) - 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam2 - 2*q1y - two3)*((0.5*q3z + one6)*(q3z + one3) + 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) - 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))*(2.0*q1y - 1.66533453693773e-16*(q1y + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) - 1.66533453693773e-16*(q1y + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) + 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))*(2.0*q1y + 1.66533453693773e-16*(q1y + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 1.66533453693773e-16*(q1y + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(8.99999999999999*q1y + 0.942809041582064*(-q1y - one3)/((q1y + one3)**2 + (q2x + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(8.99999999999999*q1y + 0.942809041582064*(-q1y - one3)/((q1y + one3)**2 + (q2x + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q1y + 0.942809041582063*(-q1y - one3)/((-q1y - one3)**2 + (-q2x - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q1y + 0.942809041582063*(-q1y - one3)/((-q1y - one3)**2 + (q2x + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q1y + 0.942809041582063*(-q1y - one3)/((q1y + one3)**2 + (-q2x - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam2 + 9.0*q1y + 3.0 + 0.942809041582063*(two3*lam2 - q1y - one3)/((-q3z - one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam2 + 9.0*q1y + 3.0 + 0.942809041582063*(two3*lam2 - q1y - one3)/((-q3z - one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam2 + 9.0*q1y + 3.0 + 0.942809041582063*(two3*lam2 - q1y - one3)/((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam2 + 8.99999999999999*q1y + 3.0 + 0.942809041582064*(two3*lam2 - q1y - one3)/((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam2 + 8.99999999999999*q1y + 3.0 + 0.942809041582064*(two3*lam2 - q1y - one3)/((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q1y + one3)**2 + (-q2x - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q1y - one3)**2 + (q2x + one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q1y - one3)**2 + (-q2x - one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q1y + one3)**2 + (q2x + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q1y + one3)**2 + (q2x + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam2 + q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam2 + q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam2 + q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam2 + q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q3z - one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam2 + q1y/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q3z - one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/4))
R[1] = 0.0222300324977063*alpha*mu*(fo - 1)**2*(-4*q2x - four3)*(q1y**2 + two3*q1y - q2x**2 - two3*q2x)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-4*q2x - four3)*(q1y**2 + two3*q1y - q2x**2 - two3*q2x)**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**3) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) + 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) - 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) - 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q1y + one3)**2 + (q2x + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q2x - two3)*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) + 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**2*((q1y + one3)**2 + (q2x + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q2x + one3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*(q2x + one3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q2x + one3)*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q2x + one3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q2x + one3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q2x + one3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*(q2x + one3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)**(3/2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + 2.0*q2x + 0.5*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-one3*lam1 + q2x + 0.25*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + one3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + 2.0*q2x + 0.5*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-one3*lam1 + q2x + 0.25*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + one3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-four3*lam1 + 4*q2x + 1.0*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 1.0*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + four3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + 2.0*q2x + 0.5*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + 2.0*q2x + 0.5*(q2x + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - 0.666666666666666)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954126*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954126*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - 0.666666666666666)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam1 - 2*q2x - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam1 + q2x + one3)/((0.75*d - 1)**2*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0889201299908253*alpha*mu*(fo - 1)**2*((-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2*(2.66666666666667*lam1 - 4*q2x - four3)/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**3) + 0.0889201299908253*alpha*mu*(fo - 1)**2*((-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))*(-four3*lam1 + 2.0*q2x + two3)/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))*(-four3*lam1 + 2.0*q2x - 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3 - 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam1 - 2*q2x - 0.666666666666666)*(-8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam1 - 2*q2x - two3)*(-8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))*(-four3*lam1 + 2.0*q2x + 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3 + 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam1 + q2x + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam1 - 2*q2x - two3)*(8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam1 - 2*q2x - two3)*(8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) - 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))*(-2.0*q2x - 1.66533453693773e-16*(q2x + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) - 1.66533453693773e-16*(q2x + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q1y + one6)*(q1y + one3) - (0.5*q2x + one6)*(q2x + one3) + 8.32667268468867e-17*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q1y + one3)**2 + (q2x + one3)**2))*(-2.0*q2x + 1.66533453693773e-16*(q2x + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) + 1.66533453693773e-16*(q2x + one3)*sqrt((q1y + one3)**2 + (q2x + one3)**2)/sqrt((q1y + one3)**2 + (q2x + one3)**2) - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q1y + one3)**2 + (q2x + one3)**2)) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(8.99999999999999*q2x + 0.942809041582064*(-q2x - one3)/((q1y + one3)**2 + (q2x + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(8.99999999999999*q2x + 0.942809041582064*(-q2x - one3)/((q1y + one3)**2 + (q2x + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q2x + 0.942809041582063*(-q2x - one3)/((-q1y - one3)**2 + (-q2x - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q2x + 0.942809041582063*(-q2x - one3)/((-q1y - one3)**2 + (q2x + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q2x + 0.942809041582063*(-q2x - one3)/((q1y + one3)**2 + (-q2x - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam1 + 9.0*q2x + 3.0 + 0.942809041582063*(two3*lam1 - q2x - one3)/((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam1 + 9.0*q2x + 3.0 + 0.942809041582063*(two3*lam1 - q2x - one3)/((-two3*lam1 + q2x + one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam1 + 9.0*q2x + 3.0 + 0.942809041582063*(two3*lam1 - q2x - one3)/((two3*lam1 - q2x - one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam1 + 8.99999999999999*q2x + 3.0 + 0.942809041582064*(two3*lam1 - q2x - one3)/((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam1 + 8.99999999999999*q2x + 3.0 + 0.942809041582064*(two3*lam1 - q2x - one3)/((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q1y + one3)**2 + (-q2x - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q1y - one3)**2 + (q2x + one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q1y - one3)**2 + (-q2x - one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q1y + one3)**2 + (q2x + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q1y + one3)**2 + (q2x + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam1 + q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam1 + q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((two3*lam1 - q2x - one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam1 + q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam1 + q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-two3*lam1 + q2x + one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam1 + q2x/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/4))
R[2] = 0.0889201299908253*alpha*mu*(fo - 1)**2*(-4*q3z - four3)*((0.5*q3z + one6)*(q3z + one3) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**3) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(q3z + one3) - 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(q3z + one3) + 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))**2/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(q3z + one3) - 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-2*q3z - two3)*((0.5*q3z + one6)*(q3z + one3) + 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))**2/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q3z + one3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q3z + one3)*((0.5*q1y + one6)*(-two3*lam2 + q1y + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*(q3z + one3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q3z + one3)*((q1y + one3)*(-two3*lam2 + q1y + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q3z + one3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0444600649954126*alpha*mu*(fo - 1)**2*(q3z + one3)*((0.5*q1y + one6)*(-one3*lam2 + 0.5*q1y + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*(q3z + one3)*((q1y + one3)*(-one3*lam2 + 0.5*q1y + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2))/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2)) + 0.0889201299908253*alpha*mu*(fo - 1)**2*(2.0*q3z + two3)*((0.5*q3z + one6)*(q3z + one3) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - 0.666666666666666)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954126*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-two3*lam1 + q2x + one3) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954126*alpha*mu*(fo - 1)**2*((0.5*q2x + one6)*(-one3*lam1 + 0.5*q2x + one6) + 0.125*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q2x + one3)*(-two3*lam1 + q2x + one3) + 0.5*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - 0.666666666666666)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) + 0.25*sqrt((q1y + one3)**2 + (q2x + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + q3z + one3)/((0.75*d - 1)**2*sqrt((q1y + one3)**2 + (q2x + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + 2.0*q3z + 0.5*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-one3*lam3 + q3z + 0.25*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-two3*lam3 + q3z + one3) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + 2.0*q3z + 0.5*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.177840259981651*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(-one3*lam3 + 0.5*q3z + one6) + 0.125*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-one3*lam3 + q3z + 0.25*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + one3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0222300324977063*alpha*mu*(fo - 1)**2*((q3z + one3)*(-two3*lam3 + q3z + one3) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-four3*lam3 + 4*q3z + 1.0*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 1.0*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + four3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + 2.0*q3z + 0.5*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))**2*(four3*lam3 - 2*q3z - two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((q3z + one3)*(-one3*lam3 + 0.5*q3z + one6) + 0.25*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))*(-two3*lam3 + 2.0*q3z + 0.5*(q3z + one3)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 0.5*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0889201299908253*alpha*mu*(fo - 1)**2*((-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2*(2.66666666666667*lam3 - 4*q3z - four3)/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**3) + 0.0889201299908253*alpha*mu*(fo - 1)**2*((-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))*(four3*lam3 - 2.0*q3z - two3)/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(-8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))*(four3*lam3 - 2.0*q3z - 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) - two3 - 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam3 - 2*q3z - two3)*(-8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam3 - 2*q3z - two3)*(-8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))*(four3*lam3 - 2.0*q3z + 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) - two3 + 1.66533453693773e-16*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*(-two3*lam3 + q3z + one3)/sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2))/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam3 - 2*q3z - two3)*(8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2) + 0.0444600649954127*alpha*mu*(fo - 1)**2*(four3*lam3 - 2*q3z - 0.666666666666666)*(8.32667268468867e-17*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)*sqrt((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2) + (-two3*lam1 + q2x + one3)*(-one3*lam1 + 0.5*q2x + one6) - (-two3*lam3 + q3z + one3)*(-one3*lam3 + 0.5*q3z + one6))**2/((0.75*d - 1)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**2*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(q3z + one3) - 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))*(2.0*q3z - 1.66533453693773e-16*(q3z + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - 1.66533453693773e-16*(q3z + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.0444600649954127*alpha*mu*(fo - 1)**2*((0.5*q3z + one6)*(q3z + one3) + 8.32667268468867e-17*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) - (-two3*lam2 + q1y + one3)*(-one3*lam2 + 0.5*q1y + one6))*(2.0*q3z + 1.66533453693773e-16*(q3z + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + 1.66533453693773e-16*(q3z + one3)*sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)/sqrt((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2) + two3)/((0.75*d - 1)**2*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(8.99999999999999*q3z + 0.942809041582064*(-q3z - one3)/((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(8.99999999999999*q3z + 0.942809041582064*(-q3z - one3)/((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q3z + 0.942809041582063*(-q3z - one3)/((-q3z - one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q3z + 0.942809041582063*(-q3z - one3)/((-q3z - one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(9.0*q3z + 0.942809041582063*(-q3z - one3)/((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/2) + 3.0)/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam3 + 9.0*q3z + 3.0 + 0.942809041582063*(two3*lam3 - q3z - one3)/((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam3 + 9.0*q3z + 3.0 + 0.942809041582063*(two3*lam3 - q3z - one3)/((-two3*lam1 + q2x + one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.235702260395516*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam3 + 9.0*q3z + 3.0 + 0.942809041582063*(two3*lam3 - q3z - one3)/((two3*lam1 - q2x - one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam3 + 8.99999999999999*q3z + 3.0 + 0.942809041582064*(two3*lam3 - q3z - one3)/((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) + 0.117851130197758*mu*(1.18518518518519 - 1.18518518518519*fo)*(-6.0*lam3 + 8.99999999999999*q3z + 3.0 + 0.942809041582064*(two3*lam3 - q3z - one3)/((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/2))/(5.65685424949238 - 4.24264068711928*d) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q3z + one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q3z - one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-q3z - one3)**2 + (-two3*lam2 + q1y + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((q3z + one3)**2 + (two3*lam2 - q1y - one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam3 + q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam3 + q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((two3*lam1 - q2x - one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/4)) - 0.686589047969039*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam3 + q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((two3*lam1 - q2x - one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam3 + q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-two3*lam1 + q2x + one3)**2 + (two3*lam3 - q3z - one3)**2)**(3/4)) - 1.37317809593808*pressure*(1.18518518518519 - 1.18518518518519*fo)*(-one3*lam3 + q3z/2 + one6)/((5.65685424949238 - 4.24264068711928*d)*((-two3*lam1 + q2x + one3)**2 + (-two3*lam3 + q3z + one3)**2)**(3/4))
return R