Hallo!
Ich kann nicht die beiden unteren Matrizen mit Mathematica 4.1.2
(Lynux, SUSE 8.0) invertieren. (No memory available. Mathematica kernel
has shut down.) Ob es um Mathematica oder um den Computer Leistungen
geht, weiß ich nicht. Kann jemand antworten?
Grüße!
Ioana
stress = {{2*(1 - alpha15)*(eta1^3/9 - (2*eta1*eta2)/3 + eta3) -
> 2*alpha17*(eta10 + (eta1^2*eta8)/9 - (2*eta1*eta9)/3) -
> ((-eta1^2/3 + eta2)*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*(eta1^4/6 - (8*eta1^2*eta2)/9 + eta2^2/6 + eta1*eta3) -
> 2*alpha17*((2*eta1*eta10)/3 + ((9*detd + eta1*eta2)*eta8)/9 -
> ((3*eta1^2 - eta2)*eta9)/6) - ((-(eta1*eta2)/3 + eta3)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*((eta1^2*eta4)/9 - (eta2*eta4)/3 - (eta1*eta5)/3 +
> eta6) - ((-(eta1*eta4)/3 + eta5)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*(-(eta1*eta11)/3 +
> (eta1*eta4*eta8)/9 - (eta4*eta9)/3 + xi1),
> -2*alpha17*((-2*eta1*eta12)/3 + eta13) - 2*(1 + alpha16)*
> (((eta1^2 - eta2)*eta4)/2 - 2*eta1*eta5 + 3*eta6),
> -2*alpha17*(-dwd2k + (eta1*eta13)/3 - (eta12*(eta1^2 - eta2))/2) -
> 2*(1 + alpha16)*(((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/2 -
> ((3*eta1^2 - eta2)*eta5)/2 + 2*eta1*eta6),
> 2*(1 - alpha15)*(-((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/6 +
> (4*eta1^2*eta5)/9 - (eta2*eta5)/3 - (eta1*eta6)/3) -
> ((-((eta1^2 - eta2)*eta4)/2 + (5*eta1*eta5)/3 - 2*eta6)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta11*(eta1^2 - 3*eta2))/6 - (eta10*eta4)/2 -
> (((6*eta1^2 - 9*eta2)*eta4 - 14*eta1*eta5 + 18*eta6)*eta8)/18 +
> ((eta1*eta4 - eta5)*eta9)/3 - (eta1*xi1)/3 + xi2),
> -2*(1 + alpha16)*(-((eta1^2 - eta2)*eta4^2)/2 + 2*eta1*eta4*eta5 -
> eta5^2 - 2*eta4*eta6) - 4*(1 - alpha15)*eta7 -
> 2*alpha17*((4*eta1*eta14)/3 - eta13*eta4 + eta12*(eta1*eta4 - eta5) -
> 2*xi4), 2*(1 + alpha16)*eta12 + 2*(1 - alpha15)*
> (eta10 + (eta1^2*eta8)/9 - (eta2*eta8)/3 - (eta1*eta9)/3) -
> 2*alpha17*((-2*eta1*eta8^2)/9 + (2*eta8*eta9)/3) -
> ((-(eta1*eta8)/3 + eta9)*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 + alpha16)*eta13 + 2*(1 - alpha15)*((4*eta1*eta10)/3 +
> 2*detd*eta8 - ((7*eta1^2 - 3*eta2)*eta9)/9) -
> 2*alpha17*(eta10*eta8 - (4*eta1*eta8*eta9)/9 + eta9^2/3) -
> ((2*eta10 - (2*eta1*eta9)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*((-2*eta1*eta12)/3 + 2*eta13) - 2*alpha17*eta12*eta8 -
> (2*eta12*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(eta1*eta11 + (-(eta1*eta4)/2 + eta5)*eta8 +
> (eta4*eta9)/2 - 3*xi1)},
> {2*(1 - alpha15)*(eta1^4/6 - (8*eta1^2*eta2)/9 + eta2^2/6 + eta1*eta3) -
> 2*alpha17*((2*eta1*eta10)/3 + ((9*detd + eta1*eta2)*eta8)/9 -
> ((3*eta1^2 - eta2)*eta9)/6) - ((-(eta1*eta2)/3 + eta3)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*((eta1*eta2^2)/9 + (eta1^5 - 5*eta1^3*eta2 +
> 5*eta1^2*eta3 + eta2*eta3)/6) - 2*alpha17*
> ((eta10*(3*eta1^2 - eta2))/6 + ((9*detd*eta1 + eta2^2)*eta8)/9 -
> ((eta1^3 - eta3)*eta9)/3) - ((eta1^4/6 - eta1^2*eta2 + eta2^2/6 +
> (4*eta1*eta3)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*((eta1*(3*eta1^2 - 7*eta2)*eta4)/18 -
> ((3*eta1^2 - eta2)*eta5)/6 + eta1*eta6) -
> ((-(eta2*eta4)/3 + eta6)*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*(-(eta11*eta2)/3 - (eta10*eta4)/3 + (eta2*eta4*eta8)/9 +
> xi2), -2*alpha17*(dwd2k + eta1*eta13 - (eta12*(3*eta1^2 + eta2))/
> 6) - 2*(1 + alpha16)*(((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/2 -
> ((3*eta1^2 - eta2)*eta5)/2 + 2*eta1*eta6),
> -2*alpha17*((eta13*(3*eta1^2 - eta2))/6 - (eta12*(eta1^3 - eta3))/3) -
> 2*(1 + alpha16)*(((4*detd*eta1 + (eta1^2 - eta2)^2)*eta4)/4 -
> (eta1^3 - eta3)*eta5 + ((3*eta1^2 - eta2)*eta6)/2),
> 2*(1 - alpha15)*(-((eta1^4 - 2*eta1^2*eta2 - eta2^2 + 2*eta1*eta3)*
> eta4)/12 + (eta1*(3*eta1^2 - eta2)*eta5)/18 - (eta2*eta6)/3) +
> 2*(1 + alpha16)*eta7 - (((-eta1^3/3 + (eta1*eta2)/2 - eta3/6)*eta4 +
> eta1^2*eta5 - (eta2*eta5)/3 - eta1*eta6)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*((eta11*(eta1^3 - 2*eta1*eta2 - eta3))/
> 6 + (eta10*(-3*eta1*eta4 + eta5))/6 -
> ((9*eta1^3*eta4 - 15*eta1*eta2*eta4 - 18*eta1^2*eta5 + 8*eta2*eta5 +
> 18*eta1*eta6)*eta8)/36 - (((-3*eta1^2 - eta2)*eta4 + 6*eta6)*
> eta9)/12 - ((3*eta1^2 - eta2)*xi1)/6 + eta1*xi2),
> -2*(1 + alpha16)*(-((2*eta1^3 - 3*eta1*eta2 + eta3)*eta4^2)/6 +
> eta1^2*eta4*eta5 - eta1*eta4*eta6 - eta5*eta6) -
> 4*(1 - alpha15)*eta1*eta7 + (2*eta7*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*((2*eta14*(3*eta1^2 - eta2))/3 +
> eta1^2*eta12*eta4 - eta13*(eta1*eta4 + eta5) - 2*eta1*xi4 - 2*xi5),
> 2*(1 + alpha16)*eta13 - 2*alpha17*((2*eta10*eta8)/3 -
> (2*eta2*eta8^2)/9) + 2*(1 - alpha15)*(eta1*eta10 +
> (eta1*(3*eta1^2 - 7*eta2)*eta8)/18 - ((3*eta1^2 - eta2)*eta9)/6) -
> ((eta10 - (eta2*eta8)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*((eta10*(3*eta1^2 + eta2))/3 + 2*detd*eta1*eta8 -
> (2*eta1*(3*eta1^2 - eta2)*eta9)/9) -
> 2*alpha17*(detd*eta8^2 - ((9*eta1^2 - eta2)*eta8*eta9)/18 +
> (eta10*(3*eta1*eta8 + eta9))/3) -
> (((-2*eta2*eta9)/3 + 2*(eta1*eta10 + ((eta1^3 - 3*eta1*eta2 + 2*eta3)*
> eta8)/6 - ((eta1^2 - eta2)*eta9)/2))*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*(1 + alpha16)*(-(eta1*eta13) +
> (eta12*(eta1^2 - eta2))/2 + xi3),
> 2*(1 - alpha15)*(2*eta1*eta13 - (eta12*(3*eta1^2 - eta2))/3) -
> 2*alpha17*eta13*eta8 - (2*eta13*(epsilon + alpha14*k + (2*eta1*k)/3))/
> k - 2*(1 + alpha16)*(eta11*eta2 + (eta10*eta4)/2 +
> (-(eta2*eta4)/2 + eta6)*eta8 - 3*xi2)},
> {2*(1 - alpha15)*((eta1^2*eta4)/9 - (eta2*eta4)/3 - (eta1*eta5)/3 +
> eta6) - ((-(eta1*eta4)/3 + eta5)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*(-(eta1*eta11)/3 +
> (eta1*eta4*eta8)/9 - (eta4*eta9)/3 + xi1),
> 2*(1 - alpha15)*((eta1*(3*eta1^2 - 7*eta2)*eta4)/18 -
> ((3*eta1^2 - eta2)*eta5)/6 + eta1*eta6) -
> ((-(eta2*eta4)/3 + eta6)*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*(-(eta11*eta2)/3 - (eta10*eta4)/3 + (eta2*eta4*eta8)/9 +
> xi2), 2*(1 - alpha15)*((eta1*eta4^2)/9 - (eta4*eta5)/6) -
> 2*alpha17*(-(eta11*eta4)/6 + (eta4^2*eta8)/9) -
> (eta4^2*(epsilon + alpha14*k + (2*eta1*k)/3))/(6*k),
> -2*alpha17*(-eta14 - (eta12*eta4)/6), -2*(1 - alpha15)*eta7 -
> 2*alpha17*(-(eta13*eta4)/6 - xi4),
> 2*(1 - alpha15)*(-((eta1^2 - eta2)*eta4^2)/12 + (4*eta1*eta4*eta5)/9 -
> eta5^2/3 - (eta4*eta6)/3) - (eta4*eta5*(epsilon + alpha14*k +
> (2*eta1*k)/3))/(6*k) - 2*alpha17*((eta11*(eta1*eta4 - 2*eta5))/6 +
> (eta4*((-3*eta1*eta4)/2 + 5*eta5)*eta8)/18 + (eta4^2*eta9)/12 -
> (eta4*xi1)/3), (-2*alpha17*eta14*eta4)/3,
> -2*alpha17*((2*eta11*eta8)/3 - (2*eta4*eta8^2)/9) -
> ((eta11 - (eta4*eta8)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/k +
> 2*(1 - alpha15)*((eta1*eta4*eta8)/9 - (eta5*eta8)/3 - (eta4*eta9)/3 +
> xi1), -(((epsilon + alpha14*k + (2*eta1*k)/3)*((-2*eta4*eta9)/3 +
> 2*xi1))/k) - 2*alpha17*((eta11*eta9)/3 - (4*eta4*eta8*eta9)/9 +
> eta8*xi1) + 2*(1 - alpha15)*(-(eta11*(eta1^2 - eta2))/2 +
> (eta10*eta4)/3 + (((eta1^2 - eta2)*eta4)/2 - eta1*eta5 + eta6)*
> eta8 + ((-7*eta1*eta4 + 3*eta5)*eta9)/9 + 2*eta1*xi1 - xi2),
> 2*(1 - alpha15)*(-eta14 - (eta12*eta4)/6)},
> {-2*alpha17*((-2*eta1*eta12)/3 + eta13) - 2*(1 + alpha16)*
> (-((eta1^2 - eta2)*eta4)/2 + 2*eta1*eta5 - 3*eta6),
> -2*alpha17*(dwd2k + eta1*eta13 - (eta12*(3*eta1^2 + eta2))/6) -
> 2*(1 + alpha16)*(-((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/2 +
> ((3*eta1^2 - eta2)*eta5)/2 - 2*eta1*eta6),
> -2*alpha17*(-eta14 - (eta12*eta4)/6),
> 2*(1 - alpha15)*(-((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/2 +
> ((3*eta1^2 - eta2)*eta5)/2 - 2*eta1*eta6) -
> ((-((eta1^2 - eta2)*eta4) + 4*eta1*eta5 - 6*eta6)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta11*(3*eta1^2 - 5*eta2))/2 - (5*eta10*eta4)/2 -
> ((3*eta1^2*eta4 - 4*eta2*eta4 - 6*eta1*eta5 + 8*eta6)*eta8)/2 -
> (-2*eta1*eta4 + eta5)*eta9 - 4*eta1*xi1 + 6*xi2),
> 2*(1 - alpha15)*(-((4*detd*eta1 + (eta1^2 - eta2)^2)*eta4)/4 +
> (eta1^3 - eta3)*eta5 - ((3*eta1^2 - eta2)*eta6)/2) -
> 6*(1 + alpha16)*eta7 - ((-((eta1^3 - 2*eta1*eta2 + eta3)*eta4) +
> (3*eta1^2 - eta2)*eta5 - 4*eta1*eta6)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*
> ((eta11*(2*eta1^3 - 3*eta1*eta2 - eta3))/2 -
> (eta10*(3*eta1*eta4 + eta5))/2 + ((-(detd*eta4) - 2*eta1^3*eta4 +
> 3*eta1*eta2*eta4 + (4*eta1^2 - eta2)*eta5 - 5*eta1*eta6)*eta8)/
> 2 + ((5*eta1^2*eta4 - eta2*eta4 - 2*eta6)*eta9)/4 -
> (3*eta1^2 - eta2)*xi1 + 4*eta1*xi2),
> -2*(1 + alpha16)*(((eta1^2 - eta2)*eta4^2)/2 - 2*eta1*eta4*eta5 +
> eta5^2 + 2*eta4*eta6) + 2*(1 - alpha15)*eta7 -
> 2*alpha17*(-(eta1*eta14) + (eta13*eta4)/2 +
> (eta12*(-3*eta1*eta4 + 2*eta5))/6 + xi4),
> 2*(1 - alpha15)*(((2*eta1^3 - 3*eta1*eta2 + eta3)*eta4^2)/6 -
> eta1^2*eta4*eta5 + (eta1*eta4 + eta5)*eta6) -
> ((-(eta4*(-((eta1^2 - eta2)*eta4)/2 + 2*eta1*eta5 - 2*eta6)) +
> 2*(((eta1^2 - eta2)*eta4^2)/4 + eta5*(-(eta1*eta4) + eta5) +
> eta4*eta6))*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta10*eta4^2)/2 - eta11*(eta1*eta5 - 2*eta6) +
> (((eta1^2 - eta2)*eta4^2)/2 - (3*eta1*eta4*eta5)/2 + eta5^2 +
> eta4*eta6)*eta8 + ((-(eta1*eta4^2) + eta4*eta5)*eta9)/2 +
> (eta1*eta4 - eta5)*xi1 - eta4*xi2), 2*(1 - alpha15)*eta13 -
> (8*alpha17*eta12*eta8)/3 - (2*eta12*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*(1 + alpha16)*(eta1*eta11 +
> (-(eta1*eta4)/2 + eta5)*eta8 + (eta4*eta9)/2 - 3*xi1),
> -2*alpha17*(eta13*eta8 + (2*eta12*eta9)/3) -
> (2*eta13*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*((3*eta11*(eta1^2 - eta2))/2 - 2*eta10*eta4 -
> 3*(((eta1^2 - eta2)*eta4)/2 - eta1*eta5 + eta6)*eta8 +
> (2*eta1*eta4 - eta5)*eta9 - 4*eta1*xi1 + 3*xi2) +
> 2*(1 - alpha15)*(eta1*eta13 - (eta12*(eta1^2 - eta2))/2 - xi3),
> -2*(1 + alpha16)*(3*eta14 - (eta12*eta4)/2) -
> ((epsilon + alpha14*k + (2*eta1*k)/3)*(-2*eta1*eta11 -
> (-(eta1*eta4) + 2*eta5)*eta8 - eta4*eta9 + 6*xi1))/k -
> 2*alpha17*(-(eta1*eta11*eta8) - (-(eta1*eta4)/2 + eta5)*eta8^2 -
> (eta4*eta8*eta9)/2 + 3*eta8*xi1) + 2*(1 - alpha15)*
> (-(eta11*eta2) - (eta10*eta4)/2 - (-(eta2*eta4)/2 + eta6)*eta8 +
> 3*xi2)}, {-2*alpha17*(-dwd2k + (eta1*eta13)/3 -
> (eta12*(eta1^2 - eta2))/2) - 2*(1 + alpha16)*
> (-((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/2 + ((3*eta1^2 - eta2)*eta5)/
> 2 - 2*eta1*eta6), -2*alpha17*((eta13*(3*eta1^2 - eta2))/6 -
> (eta12*(eta1^3 - eta3))/3) - 2*(1 + alpha16)*
> (-((4*detd*eta1 + (eta1^2 - eta2)^2)*eta4)/4 + (eta1^3 - eta3)*eta5 -
> ((3*eta1^2 - eta2)*eta6)/2), -2*(1 - alpha15)*eta7 -
> 2*alpha17*(-(eta13*eta4)/6 - xi4),
> 2*(1 - alpha15)*(-((4*detd*eta1 + (eta1^2 - eta2)^2)*eta4)/4 +
> (eta1^3 - eta3)*eta5 - ((3*eta1^2 - eta2)*eta6)/2) +
> 6*(1 + alpha16)*eta7 - ((-((eta1^3 - 2*eta1*eta2 + eta3)*eta4) +
> (3*eta1^2 - eta2)*eta5 - 4*eta1*eta6)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*
> ((eta11*(2*eta1^3 - 3*eta1*eta2 - eta3))/2 -
> (eta10*(3*eta1*eta4 + eta5))/2 + ((-(detd*eta4) - 2*eta1^3*eta4 +
> 3*eta1*eta2*eta4 + (4*eta1^2 - eta2)*eta5 - 5*eta1*eta6)*eta8)/
> 2 + ((5*eta1^2*eta4 - eta2*eta4 - 2*eta6)*eta9)/4 -
> (3*eta1^2 - eta2)*xi1 + 4*eta1*xi2),
> 2*(1 - alpha15)*(((-3*eta1^5 + 5*eta1^3*eta2 + 3*eta1^2*eta3 -
> 5*eta2*eta3)*eta4)/12 + ((5*eta1^4 + 18*eta1^2*eta2 - 3*eta2^2 -
> 32*eta1*eta3)*eta5)/12 - (eta1^3 - eta3)*eta6) -
> ((-(((5*eta1^4)/6 - 2*eta1^2*eta2 + eta2^2/2 + (2*eta1*eta3)/3)*
> eta4) + 2*(eta1^3 - eta3)*eta5 + (-3*eta1^2 + eta2)*eta6)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta11*(7*eta1^4 - 6*eta1^2*eta2 + 3*eta2^2 -
> 16*eta1*eta3))/12 - (eta10*(5*eta1^2*eta4 - 3*eta2*eta4 +
> 4*eta6))/4 + (((-9*eta1^4 + 15*eta1^2*eta2 - 6*eta2^2 +
> 6*eta1*eta3)*eta4 + 16*(eta1^3 - eta3)*eta5 -
> 12*(2*eta1^2 - eta2)*eta6)*eta8)/12 +
> (5*(eta1^3 - eta3)*eta4*eta9)/6 - 2*(eta1^3 - eta3)*xi1 +
> (3*eta1^2 - eta2)*xi2), -2*(1 + alpha16)*
> (-((-2*eta1^3 + 3*eta1*eta2 - eta3)*eta4^2)/6 - eta1^2*eta4*eta5 +
> (eta1*eta4 + eta5)*eta6) + 2*(1 - alpha15)*eta1*eta7 -
> (2*eta7*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*(-(eta14*(3*eta1^2 - eta2))/2 +
> (eta13*(3*eta1*eta4 + 8*eta5))/6 -
> (eta12*(3*eta1^2*eta4 - eta2*eta4 + 4*eta6))/4 + eta1*xi4 + 2*xi5),
> 2*(1 - alpha15)*(((eta1^4 - 3*eta2^2 + 2*eta1*eta3)*eta4^2)/12 +
> ((2*eta1^3 - 9*eta1*eta2 + eta3)*eta4*eta5)/6 -
> ((3*eta1^2 - eta2)*eta5^2)/2 - ((eta1^2 - 3*eta2)*eta4*eta6)/2 +
> 4*eta1*eta5*eta6 - 2*eta6^2) - 4*(1 + alpha16)*eta4*eta7 -
> ((-(eta4*((-eta1^3/3 + (eta1*eta2)/2 - eta3/6)*eta4 + eta1^2*eta5 -
> eta1*eta6)) + 2*((eta1^3/6 - (eta1*eta2)/4 + eta3/12)*eta4^2 -
> (eta1^2*eta4*eta5)/2 + ((eta1*eta4)/2 + eta5)*eta6))*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta10*eta4*(eta1*eta4 - eta5))/2 +
> eta11*(detd*eta4 - eta1^2*eta5 + 2*eta1*eta6) -
> ((-(eta1*(eta1^2 - eta2)*eta4^2) + 2*eta1^2*eta4*eta5 - 4*eta5*eta6)*
> eta8)/4 + (eta4*((-eta1^2 - eta2)*eta4 + 4*eta6)*eta9)/4 +
> (((eta1^2 + eta2)*eta4 - 4*eta6)*xi1)/2 + (-(eta1*eta4) + eta5)*
> xi2), (-8*alpha17*eta13*eta8)/3 -
> (2*eta13*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(eta11*eta2 + (eta10*eta4)/2 + (-(eta2*eta4)/2 + eta6)*
> eta8 - 3*xi2) + 2*(1 - alpha15)*(eta1*eta13 -
> (eta12*(eta1^2 - eta2))/2 - xi3),
> -2*alpha17*(-(dwd2k*eta8) - (eta12*(eta1^2 - eta2)*eta8)/2 +
> eta13*(eta1*eta8 + (2*eta9)/3)) - 2*(1 + alpha16)*
> ((3*eta11*(-2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(eta1*eta4 + 3*eta5))/2 -
> (3*((detd*eta4)/3 + (eta1*(eta1^2 - eta2)*eta4)/2 +
> eta1*(-(eta1*eta5) + eta6))*eta8)/2 +
> (((2*eta1^2 - eta2)*eta4 + eta6)*eta9)/2 + ((-3*eta1^2 + eta2)*xi1)/
> 2) - ((epsilon + alpha14*k + (2*eta1*k)/3)*(2*eta1*eta13 -
> eta12*(eta1^2 - eta2) - 2*xi3))/k + 2*(1 - alpha15)*
> ((eta13*(eta1^2 + eta2))/2 - (eta12*(2*eta1^3 - 3*eta1*eta2 + eta3))/
> 3 - 2*eta1*xi3), -(((epsilon + alpha14*k + (2*eta1*k)/3)*
> (-2*eta11*eta2 - eta10*eta4 - (-(eta2*eta4) + 2*eta6)*eta8 +
> 6*xi2))/k) + 2*(1 - alpha15)*((eta11*(eta1^3 - 3*eta1*eta2))/2 +
> eta10*(-(eta1*eta4) + eta5) -
> (-((eta1^2 - eta2)*(-(eta1*eta4)/2 + eta5))/2 +
> eta1*(-(eta2*eta4)/2 + eta6))*eta8 - (-((eta1^2 + eta2)*eta4)/4 +
> eta6)*eta9 - ((3*eta1^2 - eta2)*xi1)/2 + 4*eta1*xi2) -
> 2*alpha17*(-(eta11*eta2*eta8) - (eta10*eta4*eta8)/2 -
> (-(eta2*eta4)/2 + eta6)*eta8^2 + 3*eta8*xi2) -
> 2*(1 + alpha16)*(-(eta13*eta4)/2 + 3*xi4)},
> {2*(1 - alpha15)*(-((eta1^3 - 2*eta1*eta2 + eta3)*eta4)/6 +
> (4*eta1^2*eta5)/9 - (eta2*eta5)/3 - (eta1*eta6)/3) -
> ((-((eta1^2 - eta2)*eta4)/2 + (5*eta1*eta5)/3 - 2*eta6)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta11*(eta1^2 - 3*eta2))/6 - (eta10*eta4)/2 -
> (((6*eta1^2 - 9*eta2)*eta4 - 14*eta1*eta5 + 18*eta6)*eta8)/18 +
> ((eta1*eta4 - eta5)*eta9)/3 - (eta1*xi1)/3 + xi2),
> 2*(1 - alpha15)*(-((eta1^4 - 2*eta1^2*eta2 - eta2^2 + 2*eta1*eta3)*
> eta4)/12 + (eta1*(3*eta1^2 - eta2)*eta5)/18 - (eta2*eta6)/3) -
> 2*(1 + alpha16)*eta7 - (((-eta1^3/3 + (eta1*eta2)/2 - eta3/6)*eta4 +
> eta1^2*eta5 - (eta2*eta5)/3 - eta1*eta6)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*((eta11*(eta1^3 - 2*eta1*eta2 - eta3))/
> 6 + (eta10*(-3*eta1*eta4 + eta5))/6 -
> ((9*eta1^3*eta4 - 15*eta1*eta2*eta4 - 18*eta1^2*eta5 + 8*eta2*eta5 +
> 18*eta1*eta6)*eta8)/36 - (((-3*eta1^2 - eta2)*eta4 + 6*eta6)*
> eta9)/12 - ((3*eta1^2 - eta2)*xi1)/6 + eta1*xi2),
> 2*(1 - alpha15)*(-((eta1^2 - eta2)*eta4^2)/12 + (4*eta1*eta4*eta5)/9 -
> eta5^2/3 - (eta4*eta6)/3) - (eta4*eta5*(epsilon + alpha14*k +
> (2*eta1*k)/3))/(6*k) - 2*alpha17*((eta11*(eta1*eta4 - 2*eta5))/6 +
> (eta4*((-3*eta1*eta4)/2 + 5*eta5)*eta8)/18 + (eta4^2*eta9)/12 -
> (eta4*xi1)/3), -2*(1 + alpha16)*(-((eta1^2 - eta2)*eta4^2)/2 +
> 2*eta1*eta4*eta5 - eta5^2 - 2*eta4*eta6) + 2*(1 - alpha15)*eta7 -
> 2*alpha17*(-(eta1*eta14) + (eta13*eta4)/2 +
> (eta12*(-3*eta1*eta4 + 2*eta5))/6 + xi4),
> -2*(1 + alpha16)*(((-2*eta1^3 + 3*eta1*eta2 - eta3)*eta4^2)/6 +
> eta1^2*eta4*eta5 - (eta1*eta4 + eta5)*eta6) +
> 2*(1 - alpha15)*eta1*eta7 - (2*eta7*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*(-(eta14*(3*eta1^2 - eta2))/2 +
> (eta13*(3*eta1*eta4 + 8*eta5))/6 -
> (eta12*(3*eta1^2*eta4 - eta2*eta4 + 4*eta6))/4 + eta1*xi4 + 2*xi5),
> ((1 - alpha15)*eta5*((3*(eta1^2 - eta2)*eta4)/2 - 4*eta1*eta5 +
> 6*eta6))/9 - ((((eta1^2 - eta2)*eta4^2)/4 - eta5^2/3 +
> eta5*(-(eta1*eta4) + eta5) + eta4*eta6)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*
> ((eta11*(3*(eta1^2 - eta2)*eta4 - 8*eta1*eta5 + 12*eta6))/12 +
> (eta5*(-3*eta1*eta4 + 8*eta5)*eta8)/18 + (eta4*eta5*eta9)/6 -
> (2*eta5*xi1)/3), (-2*alpha17*eta14*eta5)/3 -
> 2*(1 + alpha16)*((eta4*eta5^2)/2 + eta4^2*(((eta1^2 - eta2)*eta4)/4 -
> eta1*eta5 + eta6)), 2*(1 + alpha16)*eta14 -
> ((epsilon + alpha14*k + (2*eta1*k)/3)*(eta1*eta11 - (eta5*eta8)/3 +
> (-(eta1*eta4)/2 + eta5)*eta8 + (eta4*eta9)/2 - 2*xi1))/k -
> 2*alpha17*((2*eta1*eta11*eta8)/3 + ((-3*eta1*eta4 + 4*eta5)*eta8^2)/
> 9 + (eta4*eta8*eta9)/3 - (4*eta8*xi1)/3) +
> 2*(1 - alpha15)*((eta11*(eta1^2 - eta2))/2 - (eta10*eta4)/2 -
> ((3*(eta1^2 - eta2)*eta4 - 4*eta1*eta5 + 3*eta6)*eta8)/9 -
> ((-3*eta1*eta4 + 2*eta5)*eta9)/6 - eta1*xi1 + xi2),
> -(((epsilon + alpha14*k + (2*eta1*k)/3)*(-(eta11*(-eta1^2 + eta2)) -
> 2*(((eta1^2 - eta2)*eta4)/2 - eta1*eta5 + eta6)*eta8 -
> (2*eta5*eta9)/3 - eta4*(eta10 - eta1*eta9) - 2*eta1*xi1 + 2*xi2))/
> k) + 2*(1 - alpha15)*((eta11*(2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(3*eta1*eta4 + eta5))/6 -
> ((detd*eta4 + (eta1*(eta1^2 - eta2)*eta4)/2 - eta1^2*eta5 +
> eta1*eta6)*eta8)/2 + ((3*(2*eta1^2 + eta2)*eta4 - 2*eta1*eta5 -
> 3*eta6)*eta9)/18 - ((eta1^2 + eta2)*xi1)/2 + eta1*xi2) -
> 2*alpha17*(-(eta10*eta4*eta8)/2 - (((eta1^2 - eta2)*eta4)/2 -
> eta1*eta5 + eta6)*eta8^2 + ((6*eta1*eta4 - 2*eta5)*eta8*eta9)/18 +
> (eta4*eta9^2)/6 + (eta11*(3*(eta1^2 - eta2)*eta8 + 2*eta1*eta9))/6 -
> (eta1*eta8 + (2*eta9)/3)*xi1 + eta8*xi2) - 2*(1 + alpha16)*
> (-2*eta1*eta14 + eta13*eta4 + eta12*(-(eta1*eta4) + eta5) + 2*xi4),
> -2*alpha17*eta14*eta8 - (2*eta14*(epsilon + alpha14*k + (2*eta1*k)/3))/
> k - 2*(1 + alpha16)*(eta11*(-(eta1*eta4) + eta5) -
> eta4*(-(eta1*eta4)/2 + eta5)*eta8 - (eta4^2*eta9)/2 + 2*eta4*xi1) +
> 2*(1 - alpha15)*(eta1*eta14 - (eta13*eta4)/2 -
> (eta12*(-3*eta1*eta4 + 4*eta5))/6 - xi4)},
> {-2*(1 + alpha16)*(((eta1^2 - eta2)*eta4^2)/2 - 2*eta1*eta4*eta5 +
> eta5^2 + 2*eta4*eta6) - 4*(1 - alpha15)*eta7 -
> 2*alpha17*((4*eta1*eta14)/3 - eta13*eta4 + eta12*(eta1*eta4 - eta5) -
> 2*xi4), -2*(1 + alpha16)*(((2*eta1^3 - 3*eta1*eta2 + eta3)*eta4^2)/
> 6 - eta1^2*eta4*eta5 + eta1*eta4*eta6 + eta5*eta6) -
> 4*(1 - alpha15)*eta1*eta7 + (2*eta7*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*alpha17*((2*eta14*(3*eta1^2 - eta2))/3 +
> eta1^2*eta12*eta4 - eta13*(eta1*eta4 + eta5) - 2*eta1*xi4 - 2*xi5),
> (-2*alpha17*eta14*eta4)/3,
> 2*(1 - alpha15)*(((2*eta1^3 - 3*eta1*eta2 + eta3)*eta4^2)/6 -
> eta1^2*eta4*eta5 + (eta1*eta4 + eta5)*eta6) -
> ((-(eta4*(-((eta1^2 - eta2)*eta4)/2 + 2*eta1*eta5 - 2*eta6)) +
> 2*(((eta1^2 - eta2)*eta4^2)/4 + eta5*(-(eta1*eta4) + eta5) +
> eta4*eta6))*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta10*eta4^2)/2 - eta11*(eta1*eta5 - 2*eta6) +
> (((eta1^2 - eta2)*eta4^2)/2 - (3*eta1*eta4*eta5)/2 + eta5^2 +
> eta4*eta6)*eta8 + ((-(eta1*eta4^2) + eta4*eta5)*eta9)/2 +
> (eta1*eta4 - eta5)*xi1 - eta4*xi2),
> 2*(1 - alpha15)*(((eta1^4 - 3*eta2^2 + 2*eta1*eta3)*eta4^2)/12 +
> ((2*eta1^3 - 9*eta1*eta2 + eta3)*eta4*eta5)/6 -
> ((3*eta1^2 - eta2)*eta5^2)/2 - ((eta1^2 - 3*eta2)*eta4*eta6)/2 +
> 4*eta1*eta5*eta6 - 2*eta6^2) + 4*(1 + alpha16)*eta4*eta7 -
> ((-(eta4*((-eta1^3/3 + (eta1*eta2)/2 - eta3/6)*eta4 + eta1^2*eta5 -
> eta1*eta6)) + 2*((eta1^3/6 - (eta1*eta2)/4 + eta3/12)*eta4^2 -
> (eta1^2*eta4*eta5)/2 + ((eta1*eta4)/2 + eta5)*eta6))*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*alpha17*((eta10*eta4*(eta1*eta4 - eta5))/2 +
> eta11*(detd*eta4 - eta1^2*eta5 + 2*eta1*eta6) -
> ((-(eta1*(eta1^2 - eta2)*eta4^2) + 2*eta1^2*eta4*eta5 - 4*eta5*eta6)*
> eta8)/4 + (eta4*((-eta1^2 - eta2)*eta4 + 4*eta6)*eta9)/4 +
> (((eta1^2 + eta2)*eta4 - 4*eta6)*xi1)/2 + (-(eta1*eta4) + eta5)*
> xi2), (-2*alpha17*eta14*eta5)/3 - 2*(1 + alpha16)*
> (-(eta4*eta5^2)/2 - eta4^2*(((eta1^2 - eta2)*eta4)/4 - eta1*eta5 +
> eta6)), -2*alpha17*eta11*(-((eta1^2 - eta2)*eta4^2)/2 +
> 2*eta1*eta4*eta5 - eta5^2 - 2*eta4*eta6) + 2*(1 - alpha15)*
> (-(eta1^2*eta4^2*eta5)/2 + (eta2*eta4^2*eta5)/2 +
> 2*eta1*eta4*eta5^2 - eta5^3 - 2*eta4*eta5*eta6) -
> ((-(eta4*eta5^2) + eta4^2*(-((eta1^2 - eta2)*eta4)/2 + 2*eta1*eta5 -
> 2*eta6))*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> (-8*alpha17*eta14*eta8)/3 - (2*eta14*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*(1 + alpha16)*(-(eta11*(eta1*eta4 - eta5)) -
> eta4*(-(eta1*eta4)/2 + eta5)*eta8 - (eta4^2*eta9)/2 + 2*eta4*xi1) +
> 2*(1 - alpha15)*(2*eta1*eta14 - eta13*eta4 -
> eta12*(-(eta1*eta4) + eta5) - 2*xi4),
> -2*(1 + alpha16)*(-(eta11*(eta1^2 - eta2)*eta4) + eta10*eta4^2 +
> eta4*((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8 -
> eta1*eta4^2*eta9 + 2*(eta1*eta4 + eta5)*xi1 - 2*eta4*xi2) -
> ((epsilon + alpha14*k + (2*eta1*k)/3)*(-2*eta13*eta4 +
> 4*eta1*(eta14 + (eta12*eta4)/2) - 2*eta12*eta5 - 4*xi4))/k -
> 2*alpha17*(-(eta13*eta4*eta8) + eta12*(eta1*eta4 - eta5)*eta8 +
> (2*eta14*(3*eta1*eta8 + eta9))/3 - 2*eta8*xi4) +
> 2*(1 - alpha15)*(4*eta1^2*eta14 - eta13*(2*eta1*eta4 + eta5) +
> 2*eta12*(((5*eta1^2 - eta2)*eta4)/4 - eta1*eta5 + eta6) + eta4*xi3 -
> 4*eta1*xi4 - 4*xi5), 4*(1 + alpha16)*eta14*eta4 -
> ((epsilon + alpha14*k + (2*eta1*k)/3)*(2*eta1*eta11*eta4 -
> 2*eta11*eta5 + eta4*(-(eta1*eta4) + 2*eta5)*eta8 + eta4^2*eta9 -
> 4*eta4*xi1))/k - 2*alpha17*(eta11*(eta1*eta4 - eta5)*eta8 +
> eta4*(-(eta1*eta4)/2 + eta5)*eta8^2 + (eta4^2*eta8*eta9)/2 -
> 2*eta4*eta8*xi1) + 2*(1 - alpha15)*(-(eta10*eta4^2)/2 +
> eta11*((eta1^2 - eta2)*eta4 - eta1*eta5 + 2*eta6) +
> ((-((eta1^2 - eta2)*eta4^2) + eta1*eta4*eta5 + 2*eta5^2 -
> 2*eta4*eta6)*eta8)/2 + (eta4*(eta1*eta4 + eta5)*eta9)/2 -
> (eta1*eta4 + 3*eta5)*xi1 + eta4*xi2)},
> {-2*(1 + alpha16)*eta12 + 2*(1 - alpha15)*(eta10 + (eta1^2*eta8)/9 -
> (eta2*eta8)/3 - (eta1*eta9)/3) - 2*alpha17*((-2*eta1*eta8^2)/9 +
> (2*eta8*eta9)/3) - ((-(eta1*eta8)/3 + eta9)*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k, -2*(1 + alpha16)*eta13 -
> 2*alpha17*((2*eta10*eta8)/3 - (2*eta2*eta8^2)/9) +
> 2*(1 - alpha15)*(eta1*eta10 + (eta1*(3*eta1^2 - 7*eta2)*eta8)/18 -
> ((3*eta1^2 - eta2)*eta9)/6) - ((eta10 - (eta2*eta8)/3)*
> (epsilon + alpha14*k + (2*eta1*k)/3))/k,
> -2*alpha17*((2*eta11*eta8)/3 - (2*eta4*eta8^2)/9) -
> ((eta11 - (eta4*eta8)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/k +
> 2*(1 - alpha15)*((eta1*eta4*eta8)/9 - (eta5*eta8)/3 - (eta4*eta9)/3 +
> xi1), 2*(1 - alpha15)*eta13 - (8*alpha17*eta12*eta8)/3 -
> (2*eta12*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(-(eta1*eta11) - (-(eta1*eta4)/2 + eta5)*eta8 -
> (eta4*eta9)/2 + 3*xi1), (-8*alpha17*eta13*eta8)/3 -
> (2*eta13*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(-(eta11*eta2) - (eta10*eta4)/2 -
> (-(eta2*eta4)/2 + eta6)*eta8 + 3*xi2) + 2*(1 - alpha15)*
> (eta1*eta13 - (eta12*(eta1^2 - eta2))/2 - xi3),
> -2*(1 + alpha16)*eta14 - ((epsilon + alpha14*k + (2*eta1*k)/3)*
> (eta1*eta11 - (eta5*eta8)/3 + (-(eta1*eta4)/2 + eta5)*eta8 +
> (eta4*eta9)/2 - 2*xi1))/k - 2*alpha17*((2*eta1*eta11*eta8)/3 +
> ((-3*eta1*eta4 + 4*eta5)*eta8^2)/9 + (eta4*eta8*eta9)/3 -
> (4*eta8*xi1)/3) + 2*(1 - alpha15)*((eta11*(eta1^2 - eta2))/2 -
> (eta10*eta4)/2 - ((3*(eta1^2 - eta2)*eta4 - 4*eta1*eta5 + 3*eta6)*
> eta8)/9 - ((-3*eta1*eta4 + 2*eta5)*eta9)/6 - eta1*xi1 + xi2),
> (-8*alpha17*eta14*eta8)/3 - (2*eta14*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*(1 + alpha16)*(eta11*(eta1*eta4 - eta5) +
> eta4*(-(eta1*eta4)/2 + eta5)*eta8 + (eta4^2*eta9)/2 - 2*eta4*xi1) +
> 2*(1 - alpha15)*(2*eta1*eta14 - eta13*eta4 -
> eta12*(-(eta1*eta4) + eta5) - 2*xi4), (-8*alpha17*eta8^3)/9 +
> 2*(1 - alpha15)*((eta1*eta8^2)/9 + (eta8*eta9)/3) -
> (2*eta8^2*(epsilon + alpha14*k + (2*eta1*k)/3))/(3*k),
> -2*(1 + alpha16)*eta12*eta8 - (16*alpha17*eta8^2*eta9)/9 +
> 2*(1 - alpha15)*((eta10*eta8)/3 + (2*eta1*eta8*eta9)/9 + eta9^2/3) -
> (4*eta8*eta9*(epsilon + alpha14*k + (2*eta1*k)/3))/(3*k),
> 2*(1 + alpha16)*eta11*eta8 + (2*(1 - alpha15)*eta12*eta8)/3},
> {-2*(1 + alpha16)*eta13 + 2*(1 - alpha15)*((4*eta1*eta10)/3 +
> 2*detd*eta8 - ((7*eta1^2 - 3*eta2)*eta9)/9) -
> 2*alpha17*(eta10*eta8 - (4*eta1*eta8*eta9)/9 + eta9^2/3) -
> ((2*eta10 - (2*eta1*eta9)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/k,
> 2*(1 - alpha15)*((eta10*(3*eta1^2 + eta2))/3 + 2*detd*eta1*eta8 -
> (2*eta1*(3*eta1^2 - eta2)*eta9)/9) -
> 2*alpha17*(detd*eta8^2 - ((9*eta1^2 - eta2)*eta8*eta9)/18 +
> (eta10*(3*eta1*eta8 + eta9))/3) -
> (((-2*eta2*eta9)/3 + 2*(eta1*eta10 + ((eta1^3 - 3*eta1*eta2 + 2*eta3)*
> eta8)/6 - ((eta1^2 - eta2)*eta9)/2))*(epsilon + alpha14*k +
> (2*eta1*k)/3))/k - 2*(1 + alpha16)*(eta1*eta13 -
> (eta12*(eta1^2 - eta2))/2 - xi3),
> -(((epsilon + alpha14*k + (2*eta1*k)/3)*((-2*eta4*eta9)/3 + 2*xi1))/
> k) - 2*alpha17*((eta11*eta9)/3 - (4*eta4*eta8*eta9)/9 + eta8*xi1) +
> 2*(1 - alpha15)*(-(eta11*(eta1^2 - eta2))/2 + (eta10*eta4)/3 +
> (((eta1^2 - eta2)*eta4)/2 - eta1*eta5 + eta6)*eta8 +
> ((-7*eta1*eta4 + 3*eta5)*eta9)/9 + 2*eta1*xi1 - xi2),
> -2*alpha17*(eta13*eta8 + (2*eta12*eta9)/3) -
> (2*eta13*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*((-3*eta11*(eta1^2 - eta2))/2 + 2*eta10*eta4 +
> 3*(((eta1^2 - eta2)*eta4)/2 - eta1*eta5 + eta6)*eta8 -
> (2*eta1*eta4 - eta5)*eta9 + 4*eta1*xi1 - 3*xi2) +
> 2*(1 - alpha15)*(eta1*eta13 - (eta12*(eta1^2 - eta2))/2 - xi3),
> -2*alpha17*(-(dwd2k*eta8) - (eta12*(eta1^2 - eta2)*eta8)/2 +
> eta13*(eta1*eta8 + (2*eta9)/3)) - 2*(1 + alpha16)*
> ((-3*eta11*(-2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 + 3*eta5))/2 +
> (3*((detd*eta4)/3 + (eta1*(eta1^2 - eta2)*eta4)/2 +
> eta1*(-(eta1*eta5) + eta6))*eta8)/2 -
> (((2*eta1^2 - eta2)*eta4 + eta6)*eta9)/2 - ((-3*eta1^2 + eta2)*xi1)/
> 2) - ((epsilon + alpha14*k + (2*eta1*k)/3)*(2*eta1*eta13 -
> eta12*(eta1^2 - eta2) - 2*xi3))/k + 2*(1 - alpha15)*
> ((eta13*(eta1^2 + eta2))/2 - (eta12*(2*eta1^3 - 3*eta1*eta2 + eta3))/
> 3 - 2*eta1*xi3), -(((epsilon + alpha14*k + (2*eta1*k)/3)*
> (-(eta11*(-eta1^2 + eta2)) - 2*(((eta1^2 - eta2)*eta4)/2 -
> eta1*eta5 + eta6)*eta8 - (2*eta5*eta9)/3 -
> eta4*(eta10 - eta1*eta9) - 2*eta1*xi1 + 2*xi2))/k) +
> 2*(1 - alpha15)*((eta11*(2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(3*eta1*eta4 + eta5))/6 -
> ((detd*eta4 + (eta1*(eta1^2 - eta2)*eta4)/2 - eta1^2*eta5 +
> eta1*eta6)*eta8)/2 + ((3*(2*eta1^2 + eta2)*eta4 - 2*eta1*eta5 -
> 3*eta6)*eta9)/18 - ((eta1^2 + eta2)*xi1)/2 + eta1*xi2) -
> 2*alpha17*(-(eta10*eta4*eta8)/2 - (((eta1^2 - eta2)*eta4)/2 -
> eta1*eta5 + eta6)*eta8^2 + ((6*eta1*eta4 - 2*eta5)*eta8*eta9)/18 +
> (eta4*eta9^2)/6 + (eta11*(3*(eta1^2 - eta2)*eta8 + 2*eta1*eta9))/6 -
> (eta1*eta8 + (2*eta9)/3)*xi1 + eta8*xi2) - 2*(1 + alpha16)*
> (2*eta1*eta14 - eta13*eta4 - eta12*(-(eta1*eta4) + eta5) - 2*xi4),
> -2*(1 + alpha16)*(eta11*(eta1^2 - eta2)*eta4 - eta10*eta4^2 -
> eta4*((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8 +
> eta1*eta4^2*eta9 - 2*(eta1*eta4 + eta5)*xi1 + 2*eta4*xi2) -
> ((epsilon + alpha14*k + (2*eta1*k)/3)*(-2*eta13*eta4 +
> 4*eta1*(eta14 + (eta12*eta4)/2) - 2*eta12*eta5 - 4*xi4))/k -
> 2*alpha17*(-(eta13*eta4*eta8) + eta12*(eta1*eta4 - eta5)*eta8 +
> (2*eta14*(3*eta1*eta8 + eta9))/3 - 2*eta8*xi4) +
> 2*(1 - alpha15)*(4*eta1^2*eta14 - eta13*(2*eta1*eta4 + eta5) +
> 2*eta12*(((5*eta1^2 - eta2)*eta4)/4 - eta1*eta5 + eta6) + eta4*xi3 -
> 4*eta1*xi4 - 4*xi5), 2*(1 + alpha16)*eta12*eta8 -
> (16*alpha17*eta8^2*eta9)/9 + 2*(1 - alpha15)*((eta10*eta8)/3 +
> (2*eta1*eta8*eta9)/9 + eta9^2/3) -
> (4*eta8*eta9*(epsilon + alpha14*k + (2*eta1*k)/3))/(3*k),
> 2*(1 - alpha15)*(eta1*eta10*eta8 + detd*eta8^2 + (eta10*eta9)/3 -
> ((eta1^2 - eta2)*eta8*eta9)/2 + (4*eta1*eta9^2)/9) -
> 2*alpha17*(eta10*eta8^2 + (7*eta8*eta9^2)/9) -
> ((2*eta10*eta8 + (2*eta9^2)/3)*(epsilon + alpha14*k + (2*eta1*k)/3))/
> k, -2*alpha17*eta12*eta8^2 + 2*(1 - alpha15)*
> (eta13*eta8 + (2*eta12*eta9)/3) -
> (2*eta12*eta8*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(-(eta11*eta9) - eta8*xi1)},
> {2*(1 - alpha15)*((-2*eta1*eta12)/3 + 2*eta13) - 2*alpha17*eta12*eta8 -
> (2*eta12*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(-(eta1*eta11) - (-(eta1*eta4)/2 + eta5)*eta8 -
> (eta4*eta9)/2 + 3*xi1), 2*(1 - alpha15)*(2*eta1*eta13 -
> (eta12*(3*eta1^2 - eta2))/3) - 2*alpha17*eta13*eta8 -
> (2*eta13*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(-(eta11*eta2) - (eta10*eta4)/2 -
> (-(eta2*eta4)/2 + eta6)*eta8 + 3*xi2), 2*(1 - alpha15)*
> (-eta14 - (eta12*eta4)/6), -2*(1 + alpha16)*
> (-3*eta14 + (eta12*eta4)/2) - ((epsilon + alpha14*k + (2*eta1*k)/3)*
> (-2*eta1*eta11 - (-(eta1*eta4) + 2*eta5)*eta8 - eta4*eta9 + 6*xi1))/
> k - 2*alpha17*(-(eta1*eta11*eta8) - (-(eta1*eta4)/2 + eta5)*eta8^2 -
> (eta4*eta8*eta9)/2 + 3*eta8*xi1) + 2*(1 - alpha15)*
> (-(eta11*eta2) - (eta10*eta4)/2 - (-(eta2*eta4)/2 + eta6)*eta8 +
> 3*xi2), -(((epsilon + alpha14*k + (2*eta1*k)/3)*
> (-2*eta11*eta2 - eta10*eta4 - (-(eta2*eta4) + 2*eta6)*eta8 +
> 6*xi2))/k) + 2*(1 - alpha15)*((eta11*(eta1^3 - 3*eta1*eta2))/2 +
> eta10*(-(eta1*eta4) + eta5) -
> (-((eta1^2 - eta2)*(-(eta1*eta4)/2 + eta5))/2 +
> eta1*(-(eta2*eta4)/2 + eta6))*eta8 - (-((eta1^2 + eta2)*eta4)/4 +
> eta6)*eta9 - ((3*eta1^2 - eta2)*xi1)/2 + 4*eta1*xi2) -
> 2*alpha17*(-(eta11*eta2*eta8) - (eta10*eta4*eta8)/2 -
> (-(eta2*eta4)/2 + eta6)*eta8^2 + 3*eta8*xi2) -
> 2*(1 + alpha16)*((eta13*eta4)/2 - 3*xi4), -2*alpha17*eta14*eta8 -
> (2*eta14*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(-(eta11*(-(eta1*eta4) + eta5)) +
> eta4*(-(eta1*eta4)/2 + eta5)*eta8 + (eta4^2*eta9)/2 - 2*eta4*xi1) +
> 2*(1 - alpha15)*(eta1*eta14 - (eta13*eta4)/2 -
> (eta12*(-3*eta1*eta4 + 4*eta5))/6 - xi4),
> -4*(1 + alpha16)*eta14*eta4 - ((epsilon + alpha14*k + (2*eta1*k)/3)*
> (2*eta1*eta11*eta4 - 2*eta11*eta5 + eta4*(-(eta1*eta4) + 2*eta5)*
> eta8 + eta4^2*eta9 - 4*eta4*xi1))/k -
> 2*alpha17*(eta11*(eta1*eta4 - eta5)*eta8 +
> eta4*(-(eta1*eta4)/2 + eta5)*eta8^2 + (eta4^2*eta8*eta9)/2 -
> 2*eta4*eta8*xi1) + 2*(1 - alpha15)*(-(eta10*eta4^2)/2 +
> eta11*((eta1^2 - eta2)*eta4 - eta1*eta5 + 2*eta6) +
> ((-((eta1^2 - eta2)*eta4^2) + eta1*eta4*eta5 + 2*eta5^2 -
> 2*eta4*eta6)*eta8)/2 + (eta4*(eta1*eta4 + eta5)*eta9)/2 -
> (eta1*eta4 + 3*eta5)*xi1 + eta4*xi2), -2*(1 + alpha16)*eta11*eta8 +
> (2*(1 - alpha15)*eta12*eta8)/3, -2*alpha17*eta12*eta8^2 +
> 2*(1 - alpha15)*(eta13*eta8 + (2*eta12*eta9)/3) -
> (2*eta12*eta8*(epsilon + alpha14*k + (2*eta1*k)/3))/k -
> 2*(1 + alpha16)*(eta11*eta9 + eta8*xi1), 2*alpha17*eta11*eta8^2 +
> (2*eta11*eta8*(epsilon + alpha14*k + (2*eta1*k)/3))/k +
> 2*(1 - alpha15)*(-((-(eta1*eta4)/2 + eta5)*eta8^2) -
> (eta4*eta8*eta9)/2 - eta11*(eta1*eta8 + eta9) + 2*eta8*xi1)}}
>
> heatflux = {{eta9*gamma10 + eta10*gamma11 + eta11*gamma13 + eta12*gamma14 -
> eta12*gamma15 + eta8*(gamma9 + (prod - epsilon*rho)/(2*k*rho)),
> eta10*gamma10 + (eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*
> gamma11 - eta12*gamma12 - eta13*gamma15 +
> eta9*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) + gamma13*xi1,
> (eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*gamma10 +
> ((eta10*(eta1^2 + eta2))/2 + detd*eta1*eta8 - ((eta1^3 - eta3)*eta9)/
> 3)*gamma11 - eta13*gamma12 +
> (-(eta1*eta13) + (eta12*(eta1^2 - eta2))/2)*gamma15 +
> eta10*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) + gamma13*xi2 +
> gamma14*xi3, eta12*gamma10 + eta13*gamma11 + eta11*gamma12 +
> gamma15*(eta1*eta11 - ((eta1*eta4 - 2*eta5)*eta8)/2 + (eta4*eta9)/2 -
> 2*xi1) + gamma14*xi1, (eta11*eta4*gamma13)/2 +
> (eta12*eta4*gamma14)/2 + eta14*gamma15 +
> eta11*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) + gamma10*xi1 +
> gamma11*xi2, eta13*gamma10 + (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*
> gamma11 - eta14*gamma13 + eta12*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) + gamma12*(eta1*eta11 - ((eta1*eta4 - 2*eta5)*eta8)/2 +
> (eta4*eta9)/2 - 2*xi1) + gamma15*(eta11*eta2 + (eta10*eta4)/2 -
> ((eta2*eta4 - 2*eta6)*eta8)/2 - 2*xi2) +
> gamma14*((eta11*(eta1^2 - eta2))/2 - (eta10*eta4)/2 -
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 +
> (eta1*eta4*eta9)/2 - eta1*xi1 + xi2), -(eta12*eta4*gamma13)/2 -
> eta12*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) + gamma12*xi1 +
> gamma14*(-(eta11*(eta1^2 - eta2))/2 + eta10*eta4 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 +
> (-(eta1*eta4) + eta5)*eta9 + 2*eta1*xi1 - 2*xi2) +
> gamma15*((eta11*(eta1^2 - eta2))/2 - (eta10*eta4)/2 -
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 +
> (eta1*eta4*eta9)/2 - eta1*xi1 + xi2) - gamma11*xi3},
> {eta10*gamma10 + (eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*
> gamma11 + eta12*gamma12 + eta13*gamma14 +
> eta9*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) + gamma13*xi1,
> (eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*gamma10 +
> ((eta10*(eta1^2 + eta2))/2 + detd*eta1*eta8 - ((eta1^3 - eta3)*eta9)/
> 3)*gamma11 + eta10*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma13*(-(eta11*(eta1^2 - eta2))/2 + eta10*eta4 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4 - eta5)*eta9 + 2*eta1*xi1 - 2*xi2) - gamma14*xi3 +
> gamma15*xi3, ((eta10*(eta1^2 + eta2))/2 + detd*eta1*eta8 -
> ((eta1^3 - eta3)*eta9)/3)*gamma10 +
> ((eta10*(eta1^3 + 3*eta1*eta2 + 2*eta3))/6 +
> (detd*(eta1^2 + eta2)*eta8)/2 -
> ((eta1^4 + 6*eta1^2*eta2 - 3*eta2^2 - 4*eta1*eta3)*eta9)/12)*
> gamma11 + (eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*
> (gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma13*(-(eta11*(2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 + eta5))/2 +
> ((2*detd*eta4 + eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 +
> 2*eta1*eta6)*eta8)/4 - ((eta1^2*eta4 - eta6)*eta9)/2 +
> eta1^2*xi1 - eta1*xi2) + gamma12*xi3 +
> gamma15*(detd*eta12 + eta1*xi3), eta13*gamma10 +
> (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma11 +
> (eta12*eta4*gamma13)/2 + eta12*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) + gamma12*xi1 + gamma14*xi2 +
> gamma15*((eta11*(eta1^2 - eta2))/2 - (eta10*eta4)/2 -
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 +
> (eta1*eta4*eta9)/2 - eta1*xi1 + xi2), (eta12*eta4*gamma12)/2 +
> (eta13*eta4*gamma14)/2 + (eta4*gamma13*xi1)/2 +
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi1 + gamma10*xi2 +
> gamma11*(detd*eta11 - ((eta1^2 - eta2)*xi1)/2 + eta1*xi2) +
> gamma15*(eta1*eta14 + (eta1*eta12*eta4)/2 - (eta13*eta4)/2 - xi4),
> (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma10 +
> ((eta13*(eta1^2 + eta2))/2 - (eta12*(eta1^3 - eta3))/3)*gamma11 +
> eta13*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma15*((eta11*(-2*detd + eta1^3 - eta1*eta2))/4 - (eta10*eta5)/2 -
> ((eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*eta8)/4 +
> (((eta1^2 - eta2)*eta4 + 2*eta6)*eta9)/4 - ((eta1^2 - eta2)*xi1)/
> 2) + gamma12*((eta11*(eta1^2 - eta2))/2 - (eta10*eta4)/2 -
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 +
> (eta1*eta4*eta9)/2 - eta1*xi1 + xi2) +
> gamma14*((eta11*(-2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(eta1*eta4 - eta5))/2 - (eta1*((eta1^2 - eta2)*eta4 -
> 2*eta1*eta5 + 2*eta6)*eta8)/4 + (((eta1^2 + eta2)*eta4 - 2*eta6)*
> eta9)/4 - ((eta1^2 + eta2)*xi1)/2 + eta1*xi2) +
> gamma13*(-(eta1*eta14) + (eta13*eta4)/2 - (eta12*(eta1*eta4 - 2*eta5))/
> 2 + xi4), (detd*eta11 - (detd*eta4*eta8)/2 -
> (((eta1^2 - eta2)*eta4 - 4*eta1*eta5 + 4*eta6)*eta9)/4)*gamma15 +
> gamma12*(-(eta11*(eta1^2 - eta2))/2 + eta10*eta4 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4 - eta5)*eta9 + 2*eta1*xi1 - 2*xi2) +
> gamma14*(-(eta11*(2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 + eta5))/2 +
> ((detd*eta4 + (eta1*(eta1^2 - eta2)*eta4)/2 +
> eta1*(-(eta1*eta5) + eta6))*eta8)/2 +
> ((-(eta1^2*eta4) + eta6)*eta9)/2 + eta1^2*xi1 - eta1*xi2) -
> gamma10*xi3 + gamma11*(-(detd*eta12) - eta1*xi3)},
> {(eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*gamma10 +
> ((eta10*(eta1^2 + eta2))/2 + detd*eta1*eta8 - ((eta1^3 - eta3)*eta9)/
> 3)*gamma11 + eta13*gamma12 +
> (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma14 +
> eta10*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) + gamma13*xi2 -
> gamma15*xi3, ((eta10*(eta1^2 + eta2))/2 + detd*eta1*eta8 -
> ((eta1^3 - eta3)*eta9)/3)*gamma10 +
> ((eta10*(eta1^3 + 3*eta1*eta2 + 2*eta3))/6 +
> (detd*(eta1^2 + eta2)*eta8)/2 -
> ((eta1^4 + 6*eta1^2*eta2 - 3*eta2^2 - 4*eta1*eta3)*eta9)/12)*
> gamma11 + (eta1*eta10 + detd*eta8 - ((eta1^2 - eta2)*eta9)/2)*
> (gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma13*(-(eta11*(2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 + eta5))/2 +
> ((2*detd*eta4 + eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 +
> 2*eta1*eta6)*eta8)/4 - ((eta1^2*eta4 - eta6)*eta9)/2 +
> eta1^2*xi1 - eta1*xi2) - gamma12*xi3 +
> gamma14*(-(detd*eta12) - eta1*xi3),
> ((eta10*(eta1^3 + 3*eta1*eta2 + 2*eta3))/6 +
> (detd*(eta1^2 + eta2)*eta8)/2 -
> ((eta1^4 + 6*eta1^2*eta2 - 3*eta2^2 - 4*eta1*eta3)*eta9)/12)*
> gamma10 + ((eta10*(eta1^4 + 3*eta2^2 + 8*eta1*eta3))/12 +
> (detd*(eta1^3 + 3*eta1*eta2 + 2*eta3)*eta8)/6 -
> (eta2*(eta1^3 - eta3)*eta9)/3)*gamma11 +
> ((eta10*(eta1^2 + eta2))/2 + detd*eta1*eta8 - ((eta1^3 - eta3)*eta9)/
> 3)*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma13*(-(eta11*(eta1^2 - eta2)^2)/4 +
> (eta10*((eta1^2 - eta2)*eta4 + 2*eta6))/2 +
> ((3*(eta1^2 - eta2)^2*eta4 - 4*(eta1^3 - eta3)*eta5 +
> 6*(eta1^2 - eta2)*eta6)*eta8)/12 - ((eta1^3 - eta3)*eta4*eta9)/
> 3 + (2*(eta1^3 - eta3)*xi1)/3 + (-eta1^2 + eta2)*xi2) +
> gamma14*(-(detd*eta13) - ((eta1^2 - eta2)*xi3)/2) +
> gamma15*(detd*eta13 + ((eta1^2 - eta2)*xi3)/2),
> (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma10 +
> ((eta13*(eta1^2 + eta2))/2 - (eta12*(eta1^3 - eta3))/3)*gamma11 +
> (eta13*eta4*gamma13)/2 + eta13*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) + gamma12*xi2 + gamma14*(detd*eta11 -
> ((eta1^2 - eta2)*xi1)/2 + eta1*xi2) +
> gamma15*((eta11*(-2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(eta1*eta4 - eta5))/2 - (eta1*((eta1^2 - eta2)*eta4 -
> 2*eta1*eta5 + 2*eta6)*eta8)/4 + (((eta1^2 + eta2)*eta4 - 2*eta6)*
> eta9)/4 - ((eta1^2 + eta2)*xi1)/2 + eta1*xi2),
> (eta13*eta4*gamma12)/2 - ((eta1^2*eta12 - 2*eta1*eta13 - eta12*eta2)*
> eta4*gamma14)/4 + (eta4*gamma13*xi2)/2 +
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi2 +
> gamma10*(detd*eta11 - ((eta1^2 - eta2)*xi1)/2 + eta1*xi2) +
> gamma11*(detd*eta1*eta11 - ((eta1^3 - eta3)*xi1)/3 +
> ((eta1^2 + eta2)*xi2)/2) + gamma15*(eta1^2*eta14 +
> (eta1^2*eta12*eta4)/2 + eta13*eta5 - (eta13*(eta1*eta4 + 2*eta5))/
> 2 - eta1*xi4 - xi5), ((eta13*(eta1^2 + eta2))/2 -
> (eta12*(eta1^3 - eta3))/3)*gamma10 +
> ((eta13*(eta1^3 + 3*eta1*eta2 + 2*eta3))/6 -
> (eta12*(eta1^4 + 6*eta1^2*eta2 - 3*eta2^2 - 4*eta1*eta3))/12)*
> gamma11 + (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*
> (gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma12*((eta11*(-2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(eta1*eta4 - eta5))/2 - (eta1*((eta1^2 - eta2)*eta4 -
> 2*eta1*eta5 + 2*eta6)*eta8)/4 + (((eta1^2 + eta2)*eta4 - 2*eta6)*
> eta9)/4 - ((eta1^2 + eta2)*xi1)/2 + eta1*xi2) +
> gamma15*((eta11*(-4*detd*eta1 + (eta1^2 - eta2)^2))/4 -
> (eta10*(eta1^2 - eta2)*eta4)/4 -
> ((-6*detd*eta1*eta4 + 3*(eta1^2 - eta2)^2*eta4 - 4*(eta1^3 - eta3)*
> eta5 + 6*(eta1^2 - eta2)*eta6)*eta8)/12 +
> ((eta1^3 - eta3)*eta4*eta9)/6 - ((eta1^3 - eta3)*xi1)/3 +
> ((eta1^2 - eta2)*xi2)/2) + (gamma14*(2*detd*eta1*eta11 +
> eta1^2*eta11*eta2 - eta11*eta2^2 - eta1^2*eta10*eta4 -
> eta10*eta2*eta4 + 2*eta1*eta10*eta5 - 2*detd*eta1*eta4*eta8 -
> eta1^2*eta2*eta4*eta8 + eta2^2*eta4*eta8 + 4*detd*eta5*eta8 +
> 2*eta1*eta2*eta5*eta8 - 2*eta2*eta6*eta8 + 2*detd*eta4*eta9 +
> 2*eta1*eta2*eta4*eta9 - 2*eta1*eta6*eta9 - 8*detd*xi1 -
> 4*eta1*eta2*xi1 + 2*eta1^2*xi2 + 2*eta2*xi2))/4 +
> gamma13*(-(eta1^2*eta14) - (eta1^2*eta12*eta4)/2 +
> (eta13*(eta1*eta4 + 2*eta5))/2 + eta1*xi4 + xi5),
> gamma15*(-(eta10*(((eta1^2 - eta2)*eta4)/4 - eta1*eta5 + eta6)) -
> (detd*eta4*eta9)/2 + detd*xi1) +
> gamma12*(-(eta11*(2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 + eta5))/2 +
> (((2*detd + eta1^3 - eta1*eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*
> eta8)/4 - ((eta1^2*eta4 - eta6)*eta9)/2 + eta1^2*xi1 - eta1*xi2) +
> gamma14*(-(eta11*(2*detd*eta1 + (eta1^2 - eta2)*eta2))/4 +
> (eta10*(eta2*eta4 + eta1*eta5))/2 +
> ((2*detd*eta1*eta4 + (eta1^2 - eta2)*eta2*eta4 - 2*eta1*eta2*eta5 +
> 2*eta2*eta6)*eta8)/4 - ((eta1*eta2*eta4 + (eta1^2 - eta2)*eta5 -
> eta1*eta6)*eta9)/2 + (detd + eta1*eta2)*xi1 - eta2*xi2) -
> (eta4*gamma13*xi3)/2 - (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi3 +
> gamma10*(-(detd*eta12) - eta1*xi3) + gamma11*(-(detd*eta1*eta12) -
> ((eta1^2 + eta2)*xi3)/2)}, {eta12*gamma10 + eta13*gamma11 -
> eta11*gamma12 - gamma15*xi1 + gamma14*(-(eta1*eta11) +
> ((eta1*eta4 - 2*eta5)*eta8)/2 - (eta4*eta9)/2 + 2*xi1),
> eta13*gamma10 + (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma11 +
> (eta12*eta4*gamma13)/2 + eta12*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) - gamma12*xi1 + gamma14*(-(eta11*(eta1^2 - eta2))/2 +
> (eta10*eta4)/2 + (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*
> eta8)/2 - (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) - gamma15*xi2,
> (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma10 +
> ((eta13*(eta1^2 + eta2))/2 - (eta12*(eta1^3 - eta3))/3)*gamma11 +
> (eta13*eta4*gamma13)/2 + eta13*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) - gamma12*xi2 + gamma15*(-(detd*eta11) +
> ((eta1^2 - eta2)*xi1)/2 - eta1*xi2) +
> gamma14*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 - eta5))/2 + (eta1*((eta1^2 - eta2)*eta4 -
> 2*eta1*eta5 + 2*eta6)*eta8)/4 - (((eta1^2 + eta2)*eta4 - 2*eta6)*
> eta9)/4 + ((eta1^2 + eta2)*xi1)/2 - eta1*xi2),
> -(eta11*eta4*gamma13)/2 - eta14*gamma14 + eta14*gamma15 -
> eta11*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma10*(-(eta1*eta11) + ((eta1*eta4 - 2*eta5)*eta8)/2 -
> (eta4*eta9)/2 + 2*xi1) + gamma11*(-(eta11*eta2) - (eta10*eta4)/2 +
> ((eta2*eta4 - 2*eta6)*eta8)/2 + 2*xi2), -(eta14*gamma10) -
> (eta11*eta4*gamma12)/2 + (eta4*gamma14*(-(eta1*eta11) +
> ((eta1*eta4 - 2*eta5)*eta8)/2 - (eta4*eta9)/2 + 2*xi1))/2 +
> gamma15*((eta1*eta11*eta4)/2 - eta11*eta5 +
> (eta4*(-(eta1*eta4) + 2*eta5)*eta8)/4 + (eta4^2*eta9)/4 -
> eta4*xi1) - gamma11*xi4, eta14*gamma12 +
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*(-(eta1*eta11) +
> ((eta1*eta4 - 2*eta5)*eta8)/2 - (eta4*eta9)/2 + 2*xi1) +
> gamma13*(-(eta1*eta11*eta4)/2 + (eta4*(eta1*eta4 - 2*eta5)*eta8)/4 -
> (eta4^2*eta9)/4 + eta4*xi1) + gamma10*(-(eta11*eta2) -
> (eta10*eta4)/2 + ((eta2*eta4 - 2*eta6)*eta8)/2 + 2*xi2) +
> gamma11*((eta11*(eta1^3 - 3*eta1*eta2 - eta3))/3 -
> (eta1*eta10*eta4)/2 - ((eta1*(eta1^2 - 3*eta2)*eta4 -
> 2*(eta1^2 - eta2)*eta5 + 4*eta1*eta6)*eta8)/4 +
> ((eta1^2 - eta2)*eta4*eta9)/4 - (eta1^2 - eta2)*xi1 + 2*eta1*xi2) +
> gamma15*xi4, (eta12*eta4*gamma12)/2 - (eta4*gamma13*xi1)/2 -
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi1 +
> gamma11*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 + (eta10*eta5)/2 +
> ((eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*eta8)/4 -
> (((eta1^2 - eta2)*eta4 + 2*eta6)*eta9)/4 + ((eta1^2 - eta2)*xi1)/
> 2) + gamma10*(-(eta11*(eta1^2 - eta2))/2 + (eta10*eta4)/2 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma15*(eta1*eta14 + (eta1*eta12*eta4)/2 - (eta13*eta4)/2 - xi4) +
> gamma14*(-(eta1*eta14) + (eta13*eta4)/2 - (eta12*(eta1*eta4 - 2*eta5))/
> 2 + xi4)}, {(eta11*eta4*gamma13)/2 - eta14*gamma14 -
> (eta12*eta4*gamma15)/2 + eta11*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) + gamma10*xi1 + gamma11*xi2, -(eta12*eta4*gamma12)/2 -
> (eta13*eta4*gamma15)/2 + (eta4*gamma13*xi1)/2 +
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi1 + gamma10*xi2 +
> gamma11*(detd*eta11 - ((eta1^2 - eta2)*xi1)/2 + eta1*xi2) +
> gamma14*(-(eta1*eta14) - (eta1*eta12*eta4)/2 + (eta13*eta4)/2 + xi4),
> -(eta13*eta4*gamma12)/2 + ((eta1^2*eta12 - 2*eta1*eta13 - eta12*eta2)*
> eta4*gamma15)/4 + (eta4*gamma13*xi2)/2 +
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi2 +
> gamma10*(detd*eta11 - ((eta1^2 - eta2)*xi1)/2 + eta1*xi2) +
> gamma11*(detd*eta1*eta11 - ((eta1^3 - eta3)*xi1)/3 +
> ((eta1^2 + eta2)*xi2)/2) + gamma14*(-(eta1^2*eta14) -
> (eta1^2*eta12*eta4)/2 - eta13*eta5 + (eta13*(eta1*eta4 + 2*eta5))/
> 2 + eta1*xi4 + xi5), -(eta14*gamma10) + (eta11*eta4*gamma12)/2 -
> (eta4*gamma15*(-(eta1*eta11) + ((eta1*eta4 - 2*eta5)*eta8)/2 -
> (eta4*eta9)/2 + 2*xi1))/2 + gamma14*(-(eta1*eta11*eta4)/2 +
> eta11*eta5 - (eta4*(-(eta1*eta4) + 2*eta5)*eta8)/4 -
> (eta4^2*eta9)/4 + eta4*xi1) - gamma11*xi4,
> (eta11*eta4^2*gamma13)/4 - (eta14*eta4*gamma14)/2 +
> (eta14*eta4*gamma15)/2 + (eta11*eta4*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)))/2 + gamma10*(-(eta11*(eta1*eta4 - 2*eta5))/2 +
> (eta4*(eta1*eta4 - 2*eta5)*eta8)/4 - (eta4^2*eta9)/4 + eta4*xi1) +
> gamma11*(-(eta10*eta4^2)/4 - (eta11*(eta2*eta4 - 2*eta6))/2 +
> (eta4*(eta2*eta4 - 2*eta6)*eta8)/4 + eta4*xi2),
> -(eta14*eta4*gamma13)/2 + eta11*(-((eta1^2 - eta2)*eta4)/4 +
> eta1*eta5 - eta6)*gamma14 - eta14*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) + gamma12*((eta1*eta11*eta4)/2 -
> (eta4*(eta1*eta4 - 2*eta5)*eta8)/4 + (eta4^2*eta9)/4 - eta4*xi1) -
> (eta4*gamma15*(-(eta11*eta2) - (eta10*eta4)/2 +
> ((eta2*eta4 - 2*eta6)*eta8)/2 + 2*xi2))/2 - gamma10*xi4 +
> gamma11*((eta14*(eta1^2 - eta2))/2 - eta1*xi4),
> -(eta12*eta4^2*gamma13)/4 -
> (eta12*eta4*(gamma9 + (prod - epsilon*rho)/(2*k*rho)))/2 +
> (eta4*gamma12*xi1)/2 - (eta4*gamma15*(-(eta11*(eta1^2 - eta2))/2 +
> (eta10*eta4)/2 + (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*
> eta8)/2 - (eta1*eta4*eta9)/2 + eta1*xi1 - xi2))/2 +
> gamma14*(-(eta11*(eta1^2 - eta2)*eta4)/4 + (eta10*eta4^2)/4 +
> (eta4*((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/4 -
> (eta1*eta4^2*eta9)/4 + ((eta1*eta4 + 2*eta5)*xi1)/2 -
> (eta4*xi2)/2) + gamma10*(-(eta1*eta14) - (eta1*eta12*eta4)/2 +
> (eta13*eta4)/2 + xi4) - gamma11*xi5},
> {eta13*gamma10 + (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma11 -
> eta14*gamma13 + eta12*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma12*(-(eta1*eta11) + ((eta1*eta4 - 2*eta5)*eta8)/2 -
> (eta4*eta9)/2 + 2*xi1) + gamma15*(-(eta11*(eta1^2 - eta2))/2 +
> (eta10*eta4)/2 + (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*
> eta8)/2 - (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma14*(-(eta11*eta2) - (eta10*eta4)/2 + ((eta2*eta4 - 2*eta6)*eta8)/
> 2 + 2*xi2), (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*gamma10 +
> ((eta13*(eta1^2 + eta2))/2 - (eta12*(eta1^3 - eta3))/3)*gamma11 +
> eta13*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma14*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 + (eta10*eta5)/2 +
> ((eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*eta8)/4 -
> (((eta1^2 - eta2)*eta4 + 2*eta6)*eta9)/4 + ((eta1^2 - eta2)*xi1)/
> 2) + gamma12*(-(eta11*(eta1^2 - eta2))/2 + (eta10*eta4)/2 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma15*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 - eta5))/2 + (eta1*((eta1^2 - eta2)*eta4 -
> 2*eta1*eta5 + 2*eta6)*eta8)/4 - (((eta1^2 + eta2)*eta4 - 2*eta6)*
> eta9)/4 + ((eta1^2 + eta2)*xi1)/2 - eta1*xi2) +
> gamma13*(-(eta1*eta14) + (eta13*eta4)/2 - (eta12*(eta1*eta4 - 2*eta5))/
> 2 + xi4), ((eta13*(eta1^2 + eta2))/2 - (eta12*(eta1^3 - eta3))/3)*
> gamma10 + ((eta13*(eta1^3 + 3*eta1*eta2 + 2*eta3))/6 -
> (eta12*(eta1^4 + 6*eta1^2*eta2 - 3*eta2^2 - 4*eta1*eta3))/12)*
> gamma11 + (eta1*eta13 - (eta12*(eta1^2 - eta2))/2)*
> (gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma12*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 +
> (eta10*(eta1*eta4 - eta5))/2 + (eta1*((eta1^2 - eta2)*eta4 -
> 2*eta1*eta5 + 2*eta6)*eta8)/4 - (((eta1^2 + eta2)*eta4 - 2*eta6)*
> eta9)/4 + ((eta1^2 + eta2)*xi1)/2 - eta1*xi2) +
> gamma14*(-(eta11*(-4*detd*eta1 + (eta1^2 - eta2)^2))/4 +
> (eta10*(eta1^2 - eta2)*eta4)/4 +
> ((-6*detd*eta1*eta4 + 3*(eta1^2 - eta2)^2*eta4 - 4*(eta1^3 - eta3)*
> eta5 + 6*(eta1^2 - eta2)*eta6)*eta8)/12 -
> ((eta1^3 - eta3)*eta4*eta9)/6 + ((eta1^3 - eta3)*xi1)/3 -
> ((eta1^2 - eta2)*xi2)/2) + (gamma15*(-2*detd*eta1*eta11 -
> eta1^2*eta11*eta2 + eta11*eta2^2 + eta1^2*eta10*eta4 +
> eta10*eta2*eta4 - 2*eta1*eta10*eta5 + 2*detd*eta1*eta4*eta8 +
> eta1^2*eta2*eta4*eta8 - eta2^2*eta4*eta8 - 4*detd*eta5*eta8 -
> 2*eta1*eta2*eta5*eta8 + 2*eta2*eta6*eta8 - 2*detd*eta4*eta9 -
> 2*eta1*eta2*eta4*eta9 + 2*eta1*eta6*eta9 + 8*detd*xi1 +
> 4*eta1*eta2*xi1 - 2*eta1^2*xi2 - 2*eta2*xi2))/4 +
> gamma13*(-(eta1^2*eta14) - (eta1^2*eta12*eta4)/2 +
> (eta13*(eta1*eta4 + 2*eta5))/2 + eta1*xi4 + xi5),
> -(eta14*gamma12) + (gamma9 + (prod - epsilon*rho)/(2*k*rho))*
> (-(eta1*eta11) + ((eta1*eta4 - 2*eta5)*eta8)/2 - (eta4*eta9)/2 +
> 2*xi1) + gamma13*(-(eta1*eta11*eta4)/2 +
> (eta4*(eta1*eta4 - 2*eta5)*eta8)/4 - (eta4^2*eta9)/4 + eta4*xi1) +
> gamma10*(-(eta11*eta2) - (eta10*eta4)/2 + ((eta2*eta4 - 2*eta6)*eta8)/
> 2 + 2*xi2) + gamma11*((eta11*(eta1^3 - 3*eta1*eta2 - eta3))/3 -
> (eta1*eta10*eta4)/2 - ((eta1*(eta1^2 - 3*eta2)*eta4 -
> 2*(eta1^2 - eta2)*eta5 + 4*eta1*eta6)*eta8)/4 +
> ((eta1^2 - eta2)*eta4*eta9)/4 - (eta1^2 - eta2)*xi1 + 2*eta1*xi2) -
> gamma14*xi4, -(eta14*eta4*gamma13)/2 -
> eta11*(-((eta1^2 - eta2)*eta4)/4 + eta1*eta5 - eta6)*gamma15 -
> eta14*(gamma9 + (prod - epsilon*rho)/(2*k*rho)) +
> gamma12*(-(eta1*eta11*eta4)/2 + (eta4*(eta1*eta4 - 2*eta5)*eta8)/4 -
> (eta4^2*eta9)/4 + eta4*xi1) +
> (eta4*gamma14*(-(eta11*eta2) - (eta10*eta4)/2 +
> ((eta2*eta4 - 2*eta6)*eta8)/2 + 2*xi2))/2 - gamma10*xi4 +
> gamma11*((eta14*(eta1^2 - eta2))/2 - eta1*xi4),
> gamma13*(-(eta11*((eta1^2 - eta2)*eta4 + 4*eta6))/4 +
> ((eta1*eta4 - 2*eta5)*eta5*eta8)/2 - (eta4*eta5*eta9)/2 +
> 2*eta5*xi1) + (gamma9 + (prod - epsilon*rho)/(2*k*rho))*
> (-(eta11*eta2) - (eta10*eta4)/2 + ((eta2*eta4 - 2*eta6)*eta8)/2 +
> 2*xi2) + gamma10*(-(eta11*(eta1*eta2 + (-eta1^3 + eta3)/3)) -
> (eta1*eta10*eta4)/2 - (-(eta1*eta2*eta4)/2 -
> ((eta1^2 - eta2)*(-(eta1*eta4)/2 + eta5))/2 + eta1*eta6)*eta8 +
> ((eta1^2 - eta2)*eta4*eta9)/4 - (eta1^2 - eta2)*xi1 + 2*eta1*xi2) +
> gamma11*(-(eta11*(-eta1^4 + 3*eta2^2 + 4*eta1*eta3))/6 -
> (eta10*(eta1^2 + eta2)*eta4)/4 +
> (((-2*eta1^4 + 3*eta1^2*eta2 + 3*eta2^2 + 2*eta1*eta3)*eta4 +
> 4*(eta1^3 - eta3)*eta5 - 6*(eta1^2 + eta2)*eta6)*eta8)/12 +
> ((eta1^3 - eta3)*eta4*eta9)/6 - (2*(eta1^3 - eta3)*xi1)/3 +
> (eta1^2 + eta2)*xi2) + gamma15*((eta14*(eta1^2 - eta2))/2 -
> eta13*eta5 + (eta12*((eta1^2 - eta2)*eta4 + 4*eta6))/4 - xi5) +
> gamma14*(-(eta14*(eta1^2 - eta2))/2 + eta13*eta5 -
> (eta12*((eta1^2 - eta2)*eta4 + 4*eta6))/4 + xi5),
> -(eta12*(((eta1^2 - eta2)*eta4)/4 - eta1*eta5 + eta6)*gamma15) +
> gamma10*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 + (eta10*eta5)/2 +
> ((eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*eta8)/4 -
> (((eta1^2 - eta2)*eta4 + 2*eta6)*eta9)/4 + ((eta1^2 - eta2)*xi1)/
> 2) + (gamma9 + (prod - epsilon*rho)/(2*k*rho))*
> (-(eta11*(eta1^2 - eta2))/2 + (eta10*eta4)/2 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma11*(-(eta11*(-2*detd*eta1 + eta1^2*eta2 - eta2^2))/4 -
> (eta10*((eta1^2 - eta2)*eta4 - 2*eta1*eta5))/4 +
> (((eta1^2 - eta2)*eta2*eta4 - 2*eta1*eta2*eta5 + 2*eta2*eta6)*eta8)/
> 4 - (eta1*eta6*eta9)/2 - detd*xi1 + ((eta1^2 - eta2)*xi2)/2) +
> gamma13*(-(eta11*(eta1^2 - eta2)*eta4)/4 + (eta10*eta4^2)/4 +
> (eta4*((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/4 -
> (eta1*eta4^2*eta9)/4 + (eta1*eta4*xi1)/2 - (eta4*xi2)/2) +
> gamma12*(-(eta1*eta14) + (eta13*eta4)/2 - (eta12*(eta1*eta4 - 2*eta5))/
> 2 + xi4) + gamma14*(eta12*eta6 - xi5)},
> {-(eta12*eta4*gamma13)/2 - eta12*(gamma9 + (prod - epsilon*rho)/
> (2*k*rho)) - gamma12*xi1 + gamma14*(-(eta11*(eta1^2 - eta2))/2 +
> (eta10*eta4)/2 + (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*
> eta8)/2 - (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma15*((eta11*(eta1^2 - eta2))/2 - eta10*eta4 -
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (-(eta1*eta4) + eta5)*eta9 - 2*eta1*xi1 + 2*xi2) - gamma11*xi3,
> (-(detd*eta11) + (detd*eta4*eta8)/2 +
> (((eta1^2 - eta2)*eta4 - 4*eta1*eta5 + 4*eta6)*eta9)/4)*gamma14 +
> gamma12*((eta11*(eta1^2 - eta2))/2 - eta10*eta4 -
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 +
> (eta1*eta4 - eta5)*eta9 - 2*eta1*xi1 + 2*xi2) +
> gamma15*((eta11*(2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(eta1*eta4 + eta5))/2 -
> ((detd*eta4 + (eta1*(eta1^2 - eta2)*eta4)/2 +
> eta1*(-(eta1*eta5) + eta6))*eta8)/2 -
> ((-(eta1^2*eta4) + eta6)*eta9)/2 - eta1^2*xi1 + eta1*xi2) -
> gamma10*xi3 + gamma11*(-(detd*eta12) - eta1*xi3),
> gamma14*(eta10*(((eta1^2 - eta2)*eta4)/4 - eta1*eta5 + eta6) +
> (detd*eta4*eta9)/2 - detd*xi1) +
> gamma12*((eta11*(2*detd + eta1^3 - eta1*eta2))/4 -
> (eta10*(eta1*eta4 + eta5))/2 -
> (((2*detd + eta1^3 - eta1*eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*
> eta8)/4 + ((eta1^2*eta4 - eta6)*eta9)/2 - eta1^2*xi1 + eta1*xi2) +
> gamma15*((eta11*(2*detd*eta1 + (eta1^2 - eta2)*eta2))/4 -
> (eta10*(eta2*eta4 + eta1*eta5))/2 -
> ((2*detd*eta1*eta4 + (eta1^2 - eta2)*eta2*eta4 - 2*eta1*eta2*eta5 +
> 2*eta2*eta6)*eta8)/4 + ((eta1*eta2*eta4 + (eta1^2 - eta2)*eta5 -
> eta1*eta6)*eta9)/2 - (detd + eta1*eta2)*xi1 + eta2*xi2) -
> (eta4*gamma13*xi3)/2 - (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi3 +
> gamma10*(-(detd*eta12) - eta1*xi3) + gamma11*(-(detd*eta1*eta12) -
> ((eta1^2 + eta2)*xi3)/2), -(eta12*eta4*gamma12)/2 -
> (eta4*gamma13*xi1)/2 - (gamma9 + (prod - epsilon*rho)/(2*k*rho))*xi1 +
> gamma11*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 + (eta10*eta5)/2 +
> ((eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*eta8)/4 -
> (((eta1^2 - eta2)*eta4 + 2*eta6)*eta9)/4 + ((eta1^2 - eta2)*xi1)/
> 2) + gamma10*(-(eta11*(eta1^2 - eta2))/2 + (eta10*eta4)/2 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma15*(eta1*eta14 - (eta13*eta4)/2 + (eta12*(eta1*eta4 - 2*eta5))/
> 2 - xi4) + gamma14*(-(eta1*eta14) - (eta1*eta12*eta4)/2 +
> (eta13*eta4)/2 + xi4), -(eta12*eta4^2*gamma13)/4 -
> (eta12*eta4*(gamma9 + (prod - epsilon*rho)/(2*k*rho)))/2 -
> (eta4*gamma12*xi1)/2 + (eta4*gamma14*(-(eta11*(eta1^2 - eta2))/2 +
> (eta10*eta4)/2 + (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*
> eta8)/2 - (eta1*eta4*eta9)/2 + eta1*xi1 - xi2))/2 +
> gamma15*((eta11*(eta1^2 - eta2)*eta4)/4 - (eta10*eta4^2)/4 -
> (eta4*((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/4 +
> (eta1*eta4^2*eta9)/4 - ((eta1*eta4 + 2*eta5)*xi1)/2 +
> (eta4*xi2)/2) + gamma10*(-(eta1*eta14) - (eta1*eta12*eta4)/2 +
> (eta13*eta4)/2 + xi4) - gamma11*xi5,
> eta12*(((eta1^2 - eta2)*eta4)/4 - eta1*eta5 + eta6)*gamma14 +
> gamma10*(-(eta11*(-2*detd + eta1^3 - eta1*eta2))/4 + (eta10*eta5)/2 +
> ((eta1*(eta1^2 - eta2)*eta4 - 2*eta1^2*eta5 + 2*eta1*eta6)*eta8)/4 -
> (((eta1^2 - eta2)*eta4 + 2*eta6)*eta9)/4 + ((eta1^2 - eta2)*xi1)/
> 2) + (gamma9 + (prod - epsilon*rho)/(2*k*rho))*
> (-(eta11*(eta1^2 - eta2))/2 + (eta10*eta4)/2 +
> (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/2 -
> (eta1*eta4*eta9)/2 + eta1*xi1 - xi2) +
> gamma11*(-(eta11*(-2*detd*eta1 + eta1^2*eta2 - eta2^2))/4 -
> (eta10*((eta1^2 - eta2)*eta4 - 2*eta1*eta5))/4 +
> (((eta1^2 - eta2)*eta2*eta4 - 2*eta1*eta2*eta5 + 2*eta2*eta6)*eta8)/
> 4 - (eta1*eta6*eta9)/2 - detd*xi1 + ((eta1^2 - eta2)*xi2)/2) +
> gamma13*(-(eta11*(eta1^2 - eta2)*eta4)/4 + (eta10*eta4^2)/4 +
> (eta4*((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/4 -
> (eta1*eta4^2*eta9)/4 + (eta1*eta4*xi1)/2 - (eta4*xi2)/2) +
> gamma12*(eta1*eta14 - (eta13*eta4)/2 + (eta12*(eta1*eta4 - 2*eta5))/
> 2 - xi4) + gamma15*(-(eta12*eta6) + xi5),
> (-(detd*eta11) + (detd*eta4*eta8)/2 +
> (((eta1^2 - eta2)*eta4 - 4*eta1*eta5 + 4*eta6)*eta9)/4)*gamma10 +
> gamma11*(-(detd*eta1*eta11) - (eta10*((eta1^2 - eta2)*eta4 -
> 4*eta1*eta5 + 4*eta6))/4 + (detd*eta1*eta4*eta8)/2 +
> (((eta1^3 - eta3)*eta4 - 6*eta1^2*eta5 + 6*eta1*eta6)*eta9)/6 +
> 2*detd*xi1) - (eta4*gamma13*(-(eta11*(eta1^2 - eta2))/2 +
> eta10*eta4 + (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/
> 2 - (eta1*eta4 - eta5)*eta9 + 2*eta1*xi1 - 2*xi2))/2 +
> (gamma9 + (prod - epsilon*rho)/(2*k*rho))*((eta11*(eta1^2 - eta2))/2 -
> eta10*eta4 - (((eta1^2 - eta2)*eta4 - 2*eta1*eta5 + 2*eta6)*eta8)/
> 2 + (eta1*eta4 - eta5)*eta9 - 2*eta1*xi1 + 2*xi2) +
> gamma15*(eta1^2*eta14 - (eta1*eta13*eta4)/2 +
> (eta12*((3*eta1^2 - eta2)*eta4 - 4*eta1*eta5 + 4*eta6))/4 +
> (eta4*xi3)/2 - eta1*xi4 - xi5) + gamma14*(-(eta1^2*eta14) +
> (eta1*eta13*eta4)/2 - (eta12*((3*eta1^2 - eta2)*eta4 - 4*eta1*eta5 +
> 4*eta6))/4 - (eta4*xi3)/2 + eta1*xi4 + xi5)}}
DMUG-Archiv, http://www.mathematica.ch/archiv.html