-
Notifications
You must be signed in to change notification settings - Fork 3
/
VqS.c
112 lines (107 loc) · 39.1 KB
/
VqS.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
/*
* Copyright (c) 2012, Bernard Giroux
* All rights reserved.
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions are met:
*
* * Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
* * Redistributions in binary form must reproduce the above copyright
* notice, this list of conditions and the following disclaimer in the
* documentation and/or other materials provided with the distribution.
* * Neither the name of the University of California, Berkeley nor the
* names of its contributors may be used to endorse or promote products
* derived from this software without specific prior written permission.
*
* THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND ANY
* EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
* WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
* DISCLAIMED. IN NO EVENT SHALL THE REGENTS AND CONTRIBUTORS BE LIABLE FOR ANY
* DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
* (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
* LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
* ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
* SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
#include "V.h"
double complex VqS( const double lx,
const double lz,
const double cu11,
const double cu13,
const double cu33,
const double cu55,
const double alpha1,
const double alpha3,
const double M,
const double rho,
const double rhof,
const double nk1,
const double nk3,
const double m1,
const double m3,
const double omega) {
double rhof2 = rhof*rhof;
double rhof3 = rhof2*rhof;
double rhof4 = rhof2*rhof2;
double rho2 = rho*rho;
double lx2 = lx*lx;
double lz2 = lz*lz;
double lx4 = lx2*lx2;
double lz4 = lz2*lz2;
double lx6 = lx4*lx2;
double lz6 = lz4*lz2;
double M2 = M*M;
double omega2 = omega*omega;
double alpha1_2 = alpha1*alpha1;
double alpha3_2 = alpha3*alpha3;
double cu13_2 = cu13*cu13;
return csqrt(-(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3)
/(3.*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
) + ((1 - I*csqrt(3))
*(-cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,2) + 3*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
))
/(3.*cpow(2,0.6666666666666666)*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
*cpow(-2*cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,3) + 9*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
- 27*(I*cu11*cu55*lx4*lz2*M*nk1*omega - I*cu13_2*lx2*lz4*M*nk1*omega + I*cu11*cu33*lx2*lz4*M*nk1*omega - I*2.*cu13*cu55*lx2*lz4*M*nk1*omega + I*cu33*cu55*lz6*M*nk1*omega - I*alpha1_2*cu55*lx4*lz2*M2*nk1*omega - I*alpha3_2*cu11*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu13*lx2*lz4*M2*nk1*omega - I*alpha1_2*cu33*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu55*lx2*lz4*M2*nk1*omega - I*alpha3_2*cu55*lz6*M2*nk1*omega + I*cu11*cu55*lx6*M*nk3*omega - I*cu13_2*lx4*lz2*M*nk3*omega + I*cu11*cu33*lx4*lz2*M*nk3*omega - I*2.*cu13*cu55*lx4*lz2*M*nk3*omega + I*cu33*cu55*lx2*lz4*M*nk3*omega - I*alpha1_2*cu55*lx6*M2*nk3*omega - I*alpha3_2*cu11*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu13*lx4*lz2*M2*nk3*omega - I*alpha1_2*cu33*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu55*lx4*lz2*M2*nk3*omega - I*alpha3_2*cu55*lx2*lz4*M2*nk3*omega + cu11*cu55*lx4*lz2*M*m1*omega2 - cu13_2*lx2*lz4*M*m1*omega2 + cu11*cu33*lx2*lz4*M*m1*omega2 - 2*cu13*cu55*lx2*lz4*M*m1*omega2 + cu33*cu55*lz6*M*m1*omega2 - alpha1_2*cu55*lx4*lz2*M2*m1*omega2 - alpha3_2*cu11*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu13*lx2*lz4*M2*m1*omega2 - alpha1_2*cu33*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu55*lx2*lz4*M2*m1*omega2 - alpha3_2*cu55*lz6*M2*m1*omega2 + cu11*cu55*lx6*M*m3*omega2 - cu13_2*lx4*lz2*M*m3*omega2 + cu11*cu33*lx4*lz2*M*m3*omega2 - 2*cu13*cu55*lx4*lz2*M*m3*omega2 + cu33*cu55*lx2*lz4*M*m3*omega2 - alpha1_2*cu55*lx6*M2*m3*omega2 - alpha3_2*cu11*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu13*lx4*lz2*M2*m3*omega2 - alpha1_2*cu33*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu55*lx4*lz2*M2*m3*omega2 - alpha3_2*cu55*lx2*lz4*M2*m3*omega2)
*cpow(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4,2) + csqrt(4*cpow(-cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,2) + 3*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
,3) + cpow(-2*cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,3) + 9*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
- 27*(I*cu11*cu55*lx4*lz2*M*nk1*omega - I*cu13_2*lx2*lz4*M*nk1*omega + I*cu11*cu33*lx2*lz4*M*nk1*omega - I*2.*cu13*cu55*lx2*lz4*M*nk1*omega + I*cu33*cu55*lz6*M*nk1*omega - I*alpha1_2*cu55*lx4*lz2*M2*nk1*omega - I*alpha3_2*cu11*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu13*lx2*lz4*M2*nk1*omega - I*alpha1_2*cu33*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu55*lx2*lz4*M2*nk1*omega - I*alpha3_2*cu55*lz6*M2*nk1*omega + I*cu11*cu55*lx6*M*nk3*omega - I*cu13_2*lx4*lz2*M*nk3*omega + I*cu11*cu33*lx4*lz2*M*nk3*omega - I*2.*cu13*cu55*lx4*lz2*M*nk3*omega + I*cu33*cu55*lx2*lz4*M*nk3*omega - I*alpha1_2*cu55*lx6*M2*nk3*omega - I*alpha3_2*cu11*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu13*lx4*lz2*M2*nk3*omega - I*alpha1_2*cu33*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu55*lx4*lz2*M2*nk3*omega - I*alpha3_2*cu55*lx2*lz4*M2*nk3*omega + cu11*cu55*lx4*lz2*M*m1*omega2 - cu13_2*lx2*lz4*M*m1*omega2 + cu11*cu33*lx2*lz4*M*m1*omega2 - 2*cu13*cu55*lx2*lz4*M*m1*omega2 + cu33*cu55*lz6*M*m1*omega2 - alpha1_2*cu55*lx4*lz2*M2*m1*omega2 - alpha3_2*cu11*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu13*lx2*lz4*M2*m1*omega2 - alpha1_2*cu33*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu55*lx2*lz4*M2*m1*omega2 - alpha3_2*cu55*lz6*M2*m1*omega2 + cu11*cu55*lx6*M*m3*omega2 - cu13_2*lx4*lz2*M*m3*omega2 + cu11*cu33*lx4*lz2*M*m3*omega2 - 2*cu13*cu55*lx4*lz2*M*m3*omega2 + cu33*cu55*lx2*lz4*M*m3*omega2 - alpha1_2*cu55*lx6*M2*m3*omega2 - alpha3_2*cu11*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu13*lx4*lz2*M2*m3*omega2 - alpha1_2*cu33*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu55*lx4*lz2*M2*m3*omega2 - alpha3_2*cu55*lx2*lz4*M2*m3*omega2)
*cpow(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4,2),2))
,0.3333333333333333))
- ((1 + I*csqrt(3))
*cpow(-2*cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,3) + 9*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
- 27*(I*cu11*cu55*lx4*lz2*M*nk1*omega - I*cu13_2*lx2*lz4*M*nk1*omega + I*cu11*cu33*lx2*lz4*M*nk1*omega - I*2.*cu13*cu55*lx2*lz4*M*nk1*omega + I*cu33*cu55*lz6*M*nk1*omega - I*alpha1_2*cu55*lx4*lz2*M2*nk1*omega - I*alpha3_2*cu11*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu13*lx2*lz4*M2*nk1*omega - I*alpha1_2*cu33*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu55*lx2*lz4*M2*nk1*omega - I*alpha3_2*cu55*lz6*M2*nk1*omega + I*cu11*cu55*lx6*M*nk3*omega - I*cu13_2*lx4*lz2*M*nk3*omega + I*cu11*cu33*lx4*lz2*M*nk3*omega - I*2.*cu13*cu55*lx4*lz2*M*nk3*omega + I*cu33*cu55*lx2*lz4*M*nk3*omega - I*alpha1_2*cu55*lx6*M2*nk3*omega - I*alpha3_2*cu11*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu13*lx4*lz2*M2*nk3*omega - I*alpha1_2*cu33*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu55*lx4*lz2*M2*nk3*omega - I*alpha3_2*cu55*lx2*lz4*M2*nk3*omega + cu11*cu55*lx4*lz2*M*m1*omega2 - cu13_2*lx2*lz4*M*m1*omega2 + cu11*cu33*lx2*lz4*M*m1*omega2 - 2*cu13*cu55*lx2*lz4*M*m1*omega2 + cu33*cu55*lz6*M*m1*omega2 - alpha1_2*cu55*lx4*lz2*M2*m1*omega2 - alpha3_2*cu11*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu13*lx2*lz4*M2*m1*omega2 - alpha1_2*cu33*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu55*lx2*lz4*M2*m1*omega2 - alpha3_2*cu55*lz6*M2*m1*omega2 + cu11*cu55*lx6*M*m3*omega2 - cu13_2*lx4*lz2*M*m3*omega2 + cu11*cu33*lx4*lz2*M*m3*omega2 - 2*cu13*cu55*lx4*lz2*M*m3*omega2 + cu33*cu55*lx2*lz4*M*m3*omega2 - alpha1_2*cu55*lx6*M2*m3*omega2 - alpha3_2*cu11*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu13*lx4*lz2*M2*m3*omega2 - alpha1_2*cu33*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu55*lx4*lz2*M2*m3*omega2 - alpha3_2*cu55*lx2*lz4*M2*m3*omega2)
*cpow(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4,2) + csqrt(4*cpow(-cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,2) + 3*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
,3) + cpow(-2*cpow(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3,3) + 9*(-(cu11*cu55*lx4*nk1*nk3) + cu13_2*lx2*lz2*nk1*nk3 - cu11*cu33*lx2*lz2*nk1*nk3 + 2*cu13*cu55*lx2*lz2*nk1*nk3 - cu33*cu55*lz4*nk1*nk3 + I*cu11*cu55*lx4*m3*nk1*omega - I*cu13_2*lx2*lz2*m3*nk1*omega + I*cu11*cu33*lx2*lz2*m3*nk1*omega - I*2.*cu13*cu55*lx2*lz2*m3*nk1*omega + I*cu33*cu55*lz4*m3*nk1*omega + I*cu11*cu55*lx4*m1*nk3*omega - I*cu13_2*lx2*lz2*m1*nk3*omega + I*cu11*cu33*lx2*lz2*m1*nk3*omega - I*2.*cu13*cu55*lx2*lz2*m1*nk3*omega + I*cu33*cu55*lz4*m1*nk3*omega + cu11*cu55*lx4*m1*m3*omega2 - cu13_2*lx2*lz2*m1*m3*omega2 + cu11*cu33*lx2*lz2*m1*m3*omega2 - 2*cu13*cu55*lx2*lz2*m1*m3*omega2 + cu33*cu55*lz4*m1*m3*omega2 - I*cu11*lx2*lz2*M*nk1*omega*rho - I*cu55*lx2*lz2*M*nk1*omega*rho - I*cu33*lz4*M*nk1*omega*rho - I*cu55*lz4*M*nk1*omega*rho + I*alpha1_2*lx2*lz2*M2*nk1*omega*rho + I*alpha3_2*lz4*M2*nk1*omega*rho - I*cu11*lx4*M*nk3*omega*rho - I*cu55*lx4*M*nk3*omega*rho - I*cu33*lx2*lz2*M*nk3*omega*rho - I*cu55*lx2*lz2*M*nk3*omega*rho + I*alpha1_2*lx4*M2*nk3*omega*rho + I*alpha3_2*lx2*lz2*M2*nk3*omega*rho - cu11*lx2*lz2*M*m1*omega2*rho - cu55*lx2*lz2*M*m1*omega2*rho - cu33*lz4*M*m1*omega2*rho - cu55*lz4*M*m1*omega2*rho + alpha1_2*lx2*lz2*M2*m1*omega2*rho + alpha3_2*lz4*M2*m1*omega2*rho - cu11*lx4*M*m3*omega2*rho - cu55*lx4*M*m3*omega2*rho - cu33*lx2*lz2*M*m3*omega2*rho - cu55*lx2*lz2*M*m3*omega2*rho + alpha1_2*lx4*M2*m3*omega2*rho + alpha3_2*lx2*lz2*M2*m3*omega2*rho + I*2.*alpha3*cu11*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu13*lx2*lz2*M*nk1*omega*rhof - I*2.*alpha1*cu55*lx2*lz2*M*nk1*omega*rhof + I*2.*alpha3*cu55*lz4*M*nk1*omega*rhof + I*2.*alpha1*cu55*lx4*M*nk3*omega*rhof - I*2.*alpha3*cu13*lx2*lz2*M*nk3*omega*rhof + I*2.*alpha1*cu33*lx2*lz2*M*nk3*omega*rhof - I*2.*alpha3*cu55*lx2*lz2*M*nk3*omega*rhof + 2*alpha3*cu11*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu13*lx2*lz2*M*m1*omega2*rhof - 2*alpha1*cu55*lx2*lz2*M*m1*omega2*rhof + 2*alpha3*cu55*lz4*M*m1*omega2*rhof + 2*alpha1*cu55*lx4*M*m3*omega2*rhof - 2*alpha3*cu13*lx2*lz2*M*m3*omega2*rhof + 2*alpha1*cu33*lx2*lz2*M*m3*omega2*rhof - 2*alpha3*cu55*lx2*lz2*M*m3*omega2*rhof - cu11*lx4*M*omega2*rhof2 - 2*cu13*lx2*lz2*M*omega2*rhof2 - 4*cu55*lx2*lz2*M*omega2*rhof2 - cu33*lz4*M*omega2*rhof2 + alpha1_2*lx4*M2*omega2*rhof2 + 2*alpha1*alpha3*lx2*lz2*M2*omega2*rhof2 + alpha3_2*lz4*M2*omega2*rhof2)
*(cu11*lx2*nk1*nk3*rho + cu55*lx2*nk1*nk3*rho + cu33*lz2*nk1*nk3*rho + cu55*lz2*nk1*nk3*rho - I*cu11*lx2*m3*nk1*omega*rho - I*cu55*lx2*m3*nk1*omega*rho - I*cu33*lz2*m3*nk1*omega*rho - I*cu55*lz2*m3*nk1*omega*rho - I*cu11*lx2*m1*nk3*omega*rho - I*cu55*lx2*m1*nk3*omega*rho - I*cu33*lz2*m1*nk3*omega*rho - I*cu55*lz2*m1*nk3*omega*rho - cu11*lx2*m1*m3*omega2*rho - cu55*lx2*m1*m3*omega2*rho - cu33*lz2*m1*m3*omega2*rho - cu55*lz2*m1*m3*omega2*rho + I*lz2*M*nk1*omega*rho2 + I*lx2*M*nk3*omega*rho2 + lz2*M*m1*omega2*rho2 + lx2*M*m3*omega2*rho2 - I*2.*alpha3*lz2*M*nk1*omega*rho*rhof - I*2.*alpha1*lx2*M*nk3*omega*rho*rhof - 2*alpha3*lz2*M*m1*omega2*rho*rhof - 2*alpha1*lx2*M*m3*omega2*rho*rhof - I*cu11*lx2*nk1*omega*rhof2 - I*cu55*lz2*nk1*omega*rhof2 - I*cu55*lx2*nk3*omega*rhof2 - I*cu33*lz2*nk3*omega*rhof2 - cu11*lx2*m1*omega2*rhof2 - cu55*lz2*m1*omega2*rhof2 - cu55*lx2*m3*omega2*rhof2 - cu33*lz2*m3*omega2*rhof2 + lx2*M*omega2*rho*rhof2 + lz2*M*omega2*rho*rhof2 - 2*alpha1*lx2*M*omega2*rhof3 - 2*alpha3*lz2*M*omega2*rhof3)
*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
- 27*(I*cu11*cu55*lx4*lz2*M*nk1*omega - I*cu13_2*lx2*lz4*M*nk1*omega + I*cu11*cu33*lx2*lz4*M*nk1*omega - I*2.*cu13*cu55*lx2*lz4*M*nk1*omega + I*cu33*cu55*lz6*M*nk1*omega - I*alpha1_2*cu55*lx4*lz2*M2*nk1*omega - I*alpha3_2*cu11*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu13*lx2*lz4*M2*nk1*omega - I*alpha1_2*cu33*lx2*lz4*M2*nk1*omega + I*2.*alpha1*alpha3*cu55*lx2*lz4*M2*nk1*omega - I*alpha3_2*cu55*lz6*M2*nk1*omega + I*cu11*cu55*lx6*M*nk3*omega - I*cu13_2*lx4*lz2*M*nk3*omega + I*cu11*cu33*lx4*lz2*M*nk3*omega - I*2.*cu13*cu55*lx4*lz2*M*nk3*omega + I*cu33*cu55*lx2*lz4*M*nk3*omega - I*alpha1_2*cu55*lx6*M2*nk3*omega - I*alpha3_2*cu11*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu13*lx4*lz2*M2*nk3*omega - I*alpha1_2*cu33*lx4*lz2*M2*nk3*omega + I*2.*alpha1*alpha3*cu55*lx4*lz2*M2*nk3*omega - I*alpha3_2*cu55*lx2*lz4*M2*nk3*omega + cu11*cu55*lx4*lz2*M*m1*omega2 - cu13_2*lx2*lz4*M*m1*omega2 + cu11*cu33*lx2*lz4*M*m1*omega2 - 2*cu13*cu55*lx2*lz4*M*m1*omega2 + cu33*cu55*lz6*M*m1*omega2 - alpha1_2*cu55*lx4*lz2*M2*m1*omega2 - alpha3_2*cu11*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu13*lx2*lz4*M2*m1*omega2 - alpha1_2*cu33*lx2*lz4*M2*m1*omega2 + 2*alpha1*alpha3*cu55*lx2*lz4*M2*m1*omega2 - alpha3_2*cu55*lz6*M2*m1*omega2 + cu11*cu55*lx6*M*m3*omega2 - cu13_2*lx4*lz2*M*m3*omega2 + cu11*cu33*lx4*lz2*M*m3*omega2 - 2*cu13*cu55*lx4*lz2*M*m3*omega2 + cu33*cu55*lx2*lz4*M*m3*omega2 - alpha1_2*cu55*lx6*M2*m3*omega2 - alpha3_2*cu11*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu13*lx4*lz2*M2*m3*omega2 - alpha1_2*cu33*lx4*lz2*M2*m3*omega2 + 2*alpha1*alpha3*cu55*lx4*lz2*M2*m3*omega2 - alpha3_2*cu55*lx2*lz4*M2*m3*omega2)
*cpow(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4,2),2))
,0.3333333333333333))
/(6.*cpow(2,0.3333333333333333)*(-(nk1*nk3*rho2)
+ I*m3*nk1*omega*rho2 + I*m1*nk3*omega*rho2 + m1*m3*omega2*rho2 + I*nk1*omega*rho*rhof2 + I*nk3*omega*rho*rhof2 + m1*omega2*rho*rhof2 + m3*omega2*rho*rhof2 + omega2*rhof4)
));
}