-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathdkN2_orthotropic.cpp
102 lines (101 loc) · 127 KB
/
dkN2_orthotropic.cpp
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
#include <math.h>
#include "GreenAnisotropic2D.hpp"
double medium::dkN2_orthotropic(double t, int k, int i, int j) {
if ((k==0)&&(i==1)&&(j==1)) {
return (pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1)/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==0)&&(i==2)&&(j==2)) {
return (pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1)/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==0)&&(i==1)&&(j==2)) {
return (pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==0)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,0,1,2);
}
else if ((k==1)&&(i==1)&&(j==1)) {
return -(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==1)&&(i==2)&&(j==2)) {
return -(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==1)&&(i==1)&&(j==2)) {
return (2*pow(-1, K)*R0*cos(4*Th + 2*t) - 2*T1*cos(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) - (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2);
}
else if ((k==1)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,1,1,2);
}
else if ((k==2)&&(i==1)&&(j==1)) {
return (2*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1) - (-16*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 2*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-4*pow(-1, K)*R0*cos(4*Th + 2*t) + 8*R1*cos(2*Th + 2*t) - 4*T1*cos(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==2)&&(i==2)&&(j==2)) {
return (2*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1) - (-16*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 2*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-4*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 8*R1*cos(2*Th + 2*t) + 4*T1*cos(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==2)&&(i==1)&&(j==2)) {
return (2*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t)) + (-4*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*T1*sin(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) - (4*pow(-1, K)*R0*cos(4*Th + 2*t) - 4*T1*cos(2*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2);
}
else if ((k==2)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,2,1,2);
}
else if ((k==3)&&(i==1)&&(j==1)) {
return (6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 3*(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t)) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-12*pow(-1, K)*R0*cos(4*Th + 2*t) + 24*R1*cos(2*Th + 2*t) - 12*T1*cos(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1) + (8*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) + 8*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==3)&&(i==2)&&(j==2)) {
return (6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 3*(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t)) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-12*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 24*R1*cos(2*Th + 2*t) + 12*T1*cos(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1) + (8*pow(-1, K + 1)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) - 8*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t));
}
else if ((k==3)&&(i==1)&&(j==2)) {
return (2*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(6*pow(-1, K)*R0*cos(4*Th + 2*t) - 6*T1*cos(2*t)) - (-4*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*T1*sin(2*t))*(24*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + 3*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-8*pow(-1, K)*R0*cos(4*Th + 2*t) + 8*T1*cos(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 6*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3))*(pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t));
}
else if ((k==3)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,3,1,2);
}
else if ((k==4)&&(i==1)&&(j==1)) {
return (2*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-24*pow(-1, K)*R0*cos(4*Th + 2*t) + 48*R1*cos(2*Th + 2*t) - 24*T1*cos(2*t)) - (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 4*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(8*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) + 8*T1*sin(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-24*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 4*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 4*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t)) + (16*pow(-1, K)*R0*cos(4*Th + 2*t) - 32*R1*cos(2*Th + 2*t) + 16*T1*cos(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + ((-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 6*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1);
}
else if ((k==4)&&(i==2)&&(j==2)) {
return (2*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-24*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 48*R1*cos(2*Th + 2*t) + 24*T1*cos(2*t)) - (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 4*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(8*pow(-1, K + 1)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) - 8*T1*sin(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-24*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 4*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 4*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t)) + (16*pow(-1, K + 1)*R0*cos(4*Th + 2*t) - 32*R1*cos(2*Th + 2*t) - 16*T1*cos(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + ((-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 6*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1);
}
else if ((k==4)&&(i==1)&&(j==2)) {
return (12*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 6*(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-4*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*T1*sin(2*t)) + (16*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*T1*sin(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) - (-32*pow(-1, K)*R0*cos(4*Th + 2*t) + 32*T1*cos(2*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (8*pow(-1, K)*R0*cos(4*Th + 2*t) - 8*T1*cos(2*t))*(-(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 6*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t))*((-256*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + 6*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3));
}
else if ((k==4)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,4,1,2);
}
else if ((k==5)&&(i==1)&&(j==1)) {
return (20*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 10*(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(8*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) + 8*T1*sin(2*t)) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(80*pow(-1, K)*R0*cos(4*Th + 2*t) - 160*R1*cos(2*Th + 2*t) + 80*T1*cos(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-40*pow(-1, K)*R0*cos(4*Th + 2*t) + 80*R1*cos(2*Th + 2*t) - 40*T1*cos(2*t)) + (-32*pow(-1, K)*R0*sin(4*Th + 2*t) + 64*R1*sin(2*Th + 2*t) - 32*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t))*(5*(-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 5*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 30*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 5*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1)*(-(-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 60*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 90*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3));
}
else if ((k==5)&&(i==2)&&(j==2)) {
return (20*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 10*(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(8*pow(-1, K + 1)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) - 8*T1*sin(2*t)) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(80*pow(-1, K + 1)*R0*cos(4*Th + 2*t) - 160*R1*cos(2*Th + 2*t) - 80*T1*cos(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-40*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 80*R1*cos(2*Th + 2*t) + 40*T1*cos(2*t)) + (-32*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 64*R1*sin(2*Th + 2*t) + 32*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t))*(5*(-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 5*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 30*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 5*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1)*(-(-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 60*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 90*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3));
}
else if ((k==5)&&(i==1)&&(j==2)) {
return (2*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-80*pow(-1, K)*R0*cos(4*Th + 2*t) + 80*T1*cos(2*t)) + (-4*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*T1*sin(2*t))*(-10*(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 60*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 10*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) - (16*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*T1*sin(2*t))*(40*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + 5*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (10*pow(-1, K)*R0*cos(4*Th + 2*t) - 10*T1*cos(2*t))*((-256*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + 6*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (32*pow(-1, K)*R0*cos(4*Th + 2*t) - 32*T1*cos(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t))*(-(-1920*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 60*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 120*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 90*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2));
}
else if ((k==5)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,5,1,2);
}
else if ((k==6)&&(i==1)&&(j==1)) {
return (2*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(240*pow(-1, K)*R0*cos(4*Th + 2*t) - 480*R1*cos(2*Th + 2*t) + 240*T1*cos(2*t)) - (-48*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 6*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-32*pow(-1, K)*R0*sin(4*Th + 2*t) + 64*R1*sin(2*Th + 2*t) - 32*T1*sin(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 20*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 20*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(8*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) + 8*T1*sin(2*t)) + (-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t))*(-6*(-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 720*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 360*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 540*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 6*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 6*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (-64*pow(-1, K)*R0*cos(4*Th + 2*t) + 128*R1*cos(2*Th + 2*t) - 64*T1*cos(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-60*pow(-1, K)*R0*cos(4*Th + 2*t) + 120*R1*cos(2*Th + 2*t) - 60*T1*cos(2*t))*((-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 6*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1)*((-1536*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 12*(2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 720*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 6)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(-28800*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 1800*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 480*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 1080*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 90*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 360*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 20*pow(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (-480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 30*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-256*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 90*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4));
}
else if ((k==6)&&(i==2)&&(j==2)) {
return (2*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(240*pow(-1, K + 1)*R0*cos(4*Th + 2*t) - 480*R1*cos(2*Th + 2*t) - 240*T1*cos(2*t)) - (-48*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 6*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-32*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 64*R1*sin(2*Th + 2*t) + 32*T1*sin(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 20*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 20*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(8*pow(-1, K + 1)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) - 8*T1*sin(2*t)) + (-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t))*(-6*(-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 720*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 360*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 540*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 6*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 6*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (-64*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 128*R1*cos(2*Th + 2*t) + 64*T1*cos(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-60*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 120*R1*cos(2*Th + 2*t) + 60*T1*cos(2*t))*((-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 6*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1)*((-1536*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 12*(2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 720*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 6)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(-28800*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 1800*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 480*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 1080*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 90*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 360*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 20*pow(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (-480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 30*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-256*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 90*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4));
}
else if ((k==6)&&(i==1)&&(j==2)) {
return (30*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 15*(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(16*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*T1*sin(2*t)) + (-64*pow(-1, K)*R0*sin(4*Th + 2*t) + 64*T1*sin(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-4*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*T1*sin(2*t))*(15*(-256*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 360*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 15*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 15*(-64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + 90*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (-160*pow(-1, K)*R0*cos(4*Th + 2*t) + 160*T1*cos(2*t))*(-(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 6*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (12*pow(-1, K)*R0*cos(4*Th + 2*t) - 12*T1*cos(2*t))*(-(-1920*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 60*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 120*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 90*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) - (192*pow(-1, K)*R0*cos(4*Th + 2*t) - 192*T1*cos(2*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t))*((-15360*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 480*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 20*pow(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 360*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 720*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 6)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(28800*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 1800*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 90*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 1080*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + (96*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + 12*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (-1920*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 30*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))*(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 90*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (256*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2));
}
else if ((k==6)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,6,1,2);
}
else if ((k==7)&&(i==1)&&(j==1)) {
return (42*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 21*(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-32*pow(-1, K)*R0*sin(4*Th + 2*t) + 64*R1*sin(2*Th + 2*t) - 32*T1*sin(2*t)) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-448*pow(-1, K)*R0*cos(4*Th + 2*t) + 896*R1*cos(2*Th + 2*t) - 448*T1*cos(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(560*pow(-1, K)*R0*cos(4*Th + 2*t) - 1120*R1*cos(2*Th + 2*t) + 560*T1*cos(2*t)) + (-2*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) - 2*T1*sin(2*t))*(7*(-1536*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 12*(2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 5040*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 6)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - 7*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(-28800*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 1800*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 3360*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 7560*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 7*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 90*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 7*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 360*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 140*pow(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 7*(-480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 30*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 7*(-256*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 630*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4)) + (8*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) + 8*T1*sin(2*t))*(35*(-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 840*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 35*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 210*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 35*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (128*pow(-1, K)*R0*sin(4*Th + 2*t) - 256*R1*sin(2*Th + 2*t) + 128*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-84*pow(-1, K)*R0*cos(4*Th + 2*t) + 168*R1*cos(2*Th + 2*t) - 84*T1*cos(2*t))*(-(-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 60*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 90*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (pow(-1, K)*R0*cos(4*Th + 2*t) - 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 + T1*cos(2*t) + T1)*(-(-16128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 126*(2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-3360*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 420*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(-672*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 42*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 5040*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 7)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 8) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)*(-241920*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 15120*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - 4200*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) - 12600*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(53760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 840*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-80640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 5040*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-10080*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 630*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-3584*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 14*(8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 2520*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 630*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(4480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 70*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (512*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (32768*pow(-1, K + 1)*R0*T1 + 32768*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2));
}
else if ((k==7)&&(i==2)&&(j==2)) {
return (42*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 21*(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(-32*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 64*R1*sin(2*Th + 2*t) + 32*T1*sin(2*t)) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-448*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 896*R1*cos(2*Th + 2*t) + 448*T1*cos(2*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-6*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(560*pow(-1, K + 1)*R0*cos(4*Th + 2*t) - 1120*R1*cos(2*Th + 2*t) - 560*T1*cos(2*t)) + (-2*pow(-1, K + 1)*R0*sin(4*Th + 2*t) + 4*R1*sin(2*Th + 2*t) + 2*T1*sin(2*t))*(7*(-1536*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 12*(2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 5040*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 6)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - 7*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(-28800*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 1800*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 3360*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 7560*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 7*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 90*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 7*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 360*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 140*pow(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 7*(-480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 30*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 7*(-256*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 630*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4)) + (8*pow(-1, K + 1)*R0*sin(4*Th + 2*t) - 16*R1*sin(2*Th + 2*t) - 8*T1*sin(2*t))*(35*(-64*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 840*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 35*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 210*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 35*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (128*pow(-1, K + 1)*R0*sin(4*Th + 2*t) - 256*R1*sin(2*Th + 2*t) - 128*T1*sin(2*t))/(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (-84*pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 168*R1*cos(2*Th + 2*t) + 84*T1*cos(2*t))*(-(-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 120*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 60*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 90*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (pow(-1, K + 1)*R0*cos(4*Th + 2*t) + 2*R1*cos(2*Th) - 2*R1*cos(2*Th + 2*t) + T0 - T1*cos(2*t) + T1)*(-(-16128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 126*(2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-3360*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 420*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(-672*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 42*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 5040*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 7)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 8) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)*(-241920*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 15120*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - 4200*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) - 12600*pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(53760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 840*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + pow(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-80640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 5040*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-10080*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 630*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))*(64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-3584*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 14*(8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 2520*(-8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 630*(32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(-16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(4480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 70*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (512*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (32768*pow(-1, K + 1)*R0*T1 + 32768*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K + 1)*pow(R0, 2) - 2*pow(R1, 2) + 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + pow(T0, 2) + 2*T0*T1 + (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2));
}
else if ((k==7)&&(i==1)&&(j==2)) {
return (2*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2))*(672*pow(-1, K)*R0*cos(4*Th + 2*t) - 672*T1*cos(2*t)) - (-64*pow(-1, K)*R0*sin(4*Th + 2*t) + 64*T1*sin(2*t))*(56*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + 7*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + (-4*pow(-1, K)*R0*sin(4*Th + 2*t) + 4*T1*sin(2*t))*(-21*(-1920*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 60*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 21*(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 20*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 2520*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) + 21*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(3840*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 240*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 21*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 10*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 1890*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 21*(128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (16*pow(-1, K)*R0*sin(4*Th + 2*t) - 16*T1*sin(2*t))*(-35*(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 210*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 35*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(96*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 6*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (-280*pow(-1, K)*R0*cos(4*Th + 2*t) + 280*T1*cos(2*t))*((-256*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 8*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + 24*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(576*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 36*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (-64*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - (512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) + 6*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3)) + (-128*pow(-1, K)*R0*cos(4*Th + 2*t) + 128*T1*cos(2*t))/(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t)) + (14*pow(-1, K)*R0*cos(4*Th + 2*t) - 14*T1*cos(2*t))*((-15360*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 480*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + 20*pow(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 360*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 720*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 6)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)*(28800*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 1800*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) - pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(-5760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 90*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 1080*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + (96*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + 12*(8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (-1920*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 30*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))*(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 90*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (256*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2)) + (pow(-1, K)*R0*sin(4*Th + 2*t) - 2*R1*sin(2*Th) - T1*sin(2*t))*(-(-134400*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - 4200*(128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 4)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) - (-512*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (32768*pow(-1, K + 1)*R0*T1 + 32768*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 2) - 420*pow(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t), 2)*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(80640*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 5040*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + (-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*(-4480*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 70*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) - 630*(-32*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) - (128*pow(-1, K + 1)*R0*T1 + 128*pow(R1, 2))*sin(4*Th + 4*t))*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - 5040*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 7)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 8) + pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 5)*(241920*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 15120*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 7) + pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*(-53760*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 840*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) - 12600*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 3)*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 2)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 6) - 126*pow(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t), 2)*(128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) - (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(-40320*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - 630*(512*pow(-1, K + 1)*R0*T1 + 512*pow(R1, 2))*cos(4*Th + 4*t))*(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 4) + 2520*(8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*pow(16*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + (32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t), 3)/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 5) + (8*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (8*pow(-1, K + 1)*R0*T1 + 8*pow(R1, 2))*sin(4*Th + 4*t))*(3584*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 14*(8192*pow(-1, K + 1)*R0*T1 + 8192*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3) + (128*R1*(pow(-1, K)*R0 - T0)*sin(2*Th + 2*t) + (2048*pow(-1, K + 1)*R0*T1 + 2048*pow(R1, 2))*sin(4*Th + 4*t))*(672*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) + 42*(32*pow(-1, K + 1)*R0*T1 + 32*pow(R1, 2))*cos(4*Th + 4*t))/pow(pow(-1, 2*K)*pow(R0, 2) + 2*pow(R1, 2) - 4*R1*(pow(-1, K)*R0 - T0)*cos(2*Th + 2*t) - pow(T0, 2) - 2*T0*T1 - (2*pow(-1, K + 1)*R0*T1 + 2*pow(R1, 2))*cos(4*Th + 4*t), 3));
}
else if ((k==7)&&(i==2)&&(j==1)) {
return dkN2_orthotropic(t,7,1,2);
}
return 0;
};