-
Notifications
You must be signed in to change notification settings - Fork 1
/
Lfh.m
140 lines (137 loc) · 15.6 KB
/
Lfh.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
function [Lfh] = Lfh(Y,alpha,thm,thp)
Lfh(1,1)=Y(7)*((6*alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*...
alpha(1,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) + (6*alpha(1,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 - (6*alpha(1,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 + (30*...
alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 - (30*...
alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (30*...
alpha(1,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (15*...
alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 -...
(60*alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 + (60*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 + (60*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 - (60*alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + 1) - Y(6)*((6*alpha(1,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) -...
(6*alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(1,6)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(1,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(1,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4) - Y(9)*((3*alpha(1,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (3*...
alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (3*alpha(1,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 + (3*alpha(1,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (15*...
alpha(1,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (15*...
alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (15*...
alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (15*...
alpha(1,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (30*...
alpha(1,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^3 -...
(30*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (30*alpha(1,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (30*alpha(1,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4);
Lfh(2,1)=Y(8) - Y(9)*((3*alpha(2,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm -...
thp) - (3*alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (3*alpha(2,6)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 + (3*alpha(2,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (15*...
alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (15*...
alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (15*...
alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (15*...
alpha(2,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (30*...
alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^3 -...
(30*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (30*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (30*alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4) - Y(6)*((6*alpha(2,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*...
alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(2,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(2,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(2,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4) - Y(7)*((6*alpha(2,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*...
alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(2,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(2,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(2,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(2,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(2,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(2,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(2,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4);
Lfh(3,1)=Y(9)*((3*alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (3*...
alpha(3,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) + (3*alpha(3,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 - (3*alpha(3,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 + (15*...
alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 - (15*...
alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 - (15*...
alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(3,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^3 +...
(30*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 + (30*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 - (30*alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + 1) - Y(6)*((6*alpha(3,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) -...
(6*alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(3,6)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(3,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(3,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4) - Y(7)*((6*alpha(3,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*...
alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(3,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(3,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(3,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(3,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(3,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(3,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(3,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4);
Lfh(4,1)=Y(10) - Y(9)*((3*alpha(4,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm -...
thp) - (3*alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (3*alpha(4,6)*(Y(1) +...
Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 + (3*alpha(4,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (15*...
alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (15*...
alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (15*...
alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (15*...
alpha(4,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (30*...
alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm - thp)^3 -...
(30*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (30*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (30*alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4) - Y(6)*((6*alpha(4,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*...
alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(4,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(4,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(4,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4) - Y(7)*((6*alpha(4,1)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*...
alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^5)/(thm - thp) - (6*alpha(4,6)*(Y(1) + Y(2) +...
Y(4)/2 - thp)^5)/(thm - thp)^6 + (6*alpha(4,7)*(Y(1) + Y(2) + Y(4)/2 - thp)^5)/(thm - thp)^6 - (30*...
alpha(4,2)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(Y(1) + Y(2) + Y(4)/2 - thp))/(thm - thp)^2 + (30*...
alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 - (30*...
alpha(4,6)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)*(Y(1) + Y(2) + Y(4)/2 - thp)^4)/(thm - thp)^5 + (15*...
alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^4*(2*Y(1) + 2*Y(2) + Y(4) - 2*thp))/(thm - thp)^2 +...
(60*alpha(4,3)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 - thp)^2)/(thm -...
thp)^3 - (60*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^3*(Y(1) + Y(2) + Y(4)/2 -...
thp)^2)/(thm - thp)^3 - (60*alpha(4,4)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4 + (60*alpha(4,5)*((Y(1) + Y(2) + Y(4)/2 - thp)/(thm - thp) - 1)^2*(Y(1) + Y(2) + Y(4)/2 -...
thp)^3)/(thm - thp)^4);