@@ -43,11 +43,7 @@ def t0_eta():
4343 - (delta_rho * k ** 2 * mu + k ** 2 * mu * rhog ) * D * cosh (D * k ) ** 2
4444 - (delta_rho * k * mu + k * mu * rhog ) * sinh (D * k ) * cosh (D * k )
4545 - sqrt (
46- (
47- delta_rho ** 2 * k ** 2
48- + 2 * delta_rho * k ** 2 * rhog
49- + k ** 2 * rhog ** 2
50- )
46+ (delta_rho ** 2 * k ** 2 + 2 * delta_rho * k ** 2 * rhog + k ** 2 * rhog ** 2 )
5147 * D ** 2
5248 * cosh (D * k ) ** 4
5349 - 2
@@ -61,11 +57,7 @@ def t0_eta():
6157 * sinh (D * k )
6258 * cosh (D * k ) ** 3
6359 + (
64- (
65- delta_rho ** 2 * k ** 2
66- - 2 * delta_rho * k ** 2 * rhog
67- + k ** 2 * rhog ** 2
68- )
60+ (delta_rho ** 2 * k ** 2 - 2 * delta_rho * k ** 2 * rhog + k ** 2 * rhog ** 2 )
6961 * D ** 2
7062 - 4 * delta_rho * rhog
7163 )
@@ -97,11 +89,7 @@ def t0_xi():
9789 - (delta_rho * k ** 2 * mu + k ** 2 * mu * rhog ) * D * cosh (D * k ) ** 2
9890 - (delta_rho * k * mu + k * mu * rhog ) * sinh (D * k ) * cosh (D * k )
9991 + sqrt (
100- (
101- delta_rho ** 2 * k ** 2
102- + 2 * delta_rho * k ** 2 * rhog
103- + k ** 2 * rhog ** 2
104- )
92+ (delta_rho ** 2 * k ** 2 + 2 * delta_rho * k ** 2 * rhog + k ** 2 * rhog ** 2 )
10593 * D ** 2
10694 * cosh (D * k ) ** 4
10795 - 2
@@ -115,11 +103,7 @@ def t0_xi():
115103 * sinh (D * k )
116104 * cosh (D * k ) ** 3
117105 + (
118- (
119- delta_rho ** 2 * k ** 2
120- - 2 * delta_rho * k ** 2 * rhog
121- + k ** 2 * rhog ** 2
122- )
106+ (delta_rho ** 2 * k ** 2 - 2 * delta_rho * k ** 2 * rhog + k ** 2 * rhog ** 2 )
123107 * D ** 2
124108 - 4 * delta_rho * rhog
125109 )
@@ -175,17 +159,12 @@ def nond_F(x, t):
175159 + delta_rho * k ** 2
176160 + k ** 2 * rhog
177161 - sqrt (
178- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
179- * sinh (k ) ** 4
162+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
180163 + delta_rho ** 2 * k ** 2
181164 + 2 * delta_rho * k ** 2 * rhog
182165 + k ** 2 * rhog ** 2
183166 + 2
184- * (
185- delta_rho ** 2 * k
186- + 2 * delta_rho * k * rhog
187- + k * rhog ** 2
188- )
167+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
189168 * sinh (k )
190169 * cosh (k )
191170 + (
@@ -273,17 +252,12 @@ def nond_F(x, t):
273252 + delta_rho * k ** 2
274253 + k ** 2 * rhog
275254 + sqrt (
276- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
277- * sinh (k ) ** 4
255+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
278256 + delta_rho ** 2 * k ** 2
279257 + 2 * delta_rho * k ** 2 * rhog
280258 + k ** 2 * rhog ** 2
281259 + 2
282- * (
283- delta_rho ** 2 * k
284- + 2 * delta_rho * k * rhog
285- + k * rhog ** 2
286- )
260+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
287261 * sinh (k )
288262 * cosh (k )
289263 + (
@@ -385,11 +359,7 @@ def nond_F(x, t):
385359 + (delta_rho * k - k * rhog ) * cosh (k ) ** 2
386360 - sqrt (
387361 delta_rho ** 2 * k ** 2 * sinh (k ) ** 4
388- - 2
389- * delta_rho ** 2
390- * k ** 2
391- * sinh (k ) ** 2
392- * cosh (k ) ** 2
362+ - 2 * delta_rho ** 2 * k ** 2 * sinh (k ) ** 2 * cosh (k ) ** 2
393363 + delta_rho ** 2 * k ** 2 * cosh (k ) ** 4
394364 - 2 * delta_rho * k ** 2 * rhog * sinh (k ) ** 4
395365 + 2 * delta_rho * k ** 2 * rhog * cosh (k ) ** 4
@@ -420,11 +390,7 @@ def nond_F(x, t):
420390 + (delta_rho * k - k * rhog ) * cosh (k ) ** 2
421391 + sqrt (
422392 delta_rho ** 2 * k ** 2 * sinh (k ) ** 4
423- - 2
424- * delta_rho ** 2
425- * k ** 2
426- * sinh (k ) ** 2
427- * cosh (k ) ** 2
393+ - 2 * delta_rho ** 2 * k ** 2 * sinh (k ) ** 2 * cosh (k ) ** 2
428394 + delta_rho ** 2 * k ** 2 * cosh (k ) ** 4
429395 - 2 * delta_rho * k ** 2 * rhog * sinh (k ) ** 4
430396 + 2 * delta_rho * k ** 2 * rhog * cosh (k ) ** 4
@@ -482,17 +448,12 @@ def nond_F(x, t):
482448 + delta_rho * k ** 2
483449 + k ** 2 * rhog
484450 + sqrt (
485- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
486- * sinh (k ) ** 4
451+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
487452 + delta_rho ** 2 * k ** 2
488453 + 2 * delta_rho * k ** 2 * rhog
489454 + k ** 2 * rhog ** 2
490455 + 2
491- * (
492- delta_rho ** 2 * k
493- + 2 * delta_rho * k * rhog
494- + k * rhog ** 2
495- )
456+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
496457 * sinh (k )
497458 * cosh (k )
498459 + (
@@ -522,17 +483,12 @@ def nond_F(x, t):
522483 + delta_rho * k ** 2
523484 + k ** 2 * rhog
524485 - sqrt (
525- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
526- * sinh (k ) ** 4
486+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
527487 + delta_rho ** 2 * k ** 2
528488 + 2 * delta_rho * k ** 2 * rhog
529489 + k ** 2 * rhog ** 2
530490 + 2
531- * (
532- delta_rho ** 2 * k
533- + 2 * delta_rho * k * rhog
534- + k * rhog ** 2
535- )
491+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
536492 * sinh (k )
537493 * cosh (k )
538494 + (
@@ -666,17 +622,12 @@ def nond_G(x, t):
666622 + delta_rho * k ** 2
667623 + k ** 2 * rhog
668624 - sqrt (
669- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
670- * sinh (k ) ** 4
625+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
671626 + delta_rho ** 2 * k ** 2
672627 + 2 * delta_rho * k ** 2 * rhog
673628 + k ** 2 * rhog ** 2
674629 + 2
675- * (
676- delta_rho ** 2 * k
677- + 2 * delta_rho * k * rhog
678- + k * rhog ** 2
679- )
630+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
680631 * sinh (k )
681632 * cosh (k )
682633 + (
@@ -797,17 +748,12 @@ def nond_G(x, t):
797748 + delta_rho * k ** 2
798749 + k ** 2 * rhog
799750 + sqrt (
800- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
801- * sinh (k ) ** 4
751+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
802752 + delta_rho ** 2 * k ** 2
803753 + 2 * delta_rho * k ** 2 * rhog
804754 + k ** 2 * rhog ** 2
805755 + 2
806- * (
807- delta_rho ** 2 * k
808- + 2 * delta_rho * k * rhog
809- + k * rhog ** 2
810- )
756+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
811757 * sinh (k )
812758 * cosh (k )
813759 + (
@@ -937,11 +883,7 @@ def nond_G(x, t):
937883 + (delta_rho * k - k * rhog ) * cosh (k ) ** 2
938884 - sqrt (
939885 delta_rho ** 2 * k ** 2 * sinh (k ) ** 4
940- - 2
941- * delta_rho ** 2
942- * k ** 2
943- * sinh (k ) ** 2
944- * cosh (k ) ** 2
886+ - 2 * delta_rho ** 2 * k ** 2 * sinh (k ) ** 2 * cosh (k ) ** 2
945887 + delta_rho ** 2 * k ** 2 * cosh (k ) ** 4
946888 - 2 * delta_rho * k ** 2 * rhog * sinh (k ) ** 4
947889 + 2 * delta_rho * k ** 2 * rhog * cosh (k ) ** 4
@@ -972,11 +914,7 @@ def nond_G(x, t):
972914 + (delta_rho * k - k * rhog ) * cosh (k ) ** 2
973915 + sqrt (
974916 delta_rho ** 2 * k ** 2 * sinh (k ) ** 4
975- - 2
976- * delta_rho ** 2
977- * k ** 2
978- * sinh (k ) ** 2
979- * cosh (k ) ** 2
917+ - 2 * delta_rho ** 2 * k ** 2 * sinh (k ) ** 2 * cosh (k ) ** 2
980918 + delta_rho ** 2 * k ** 2 * cosh (k ) ** 4
981919 - 2 * delta_rho * k ** 2 * rhog * sinh (k ) ** 4
982920 + 2 * delta_rho * k ** 2 * rhog * cosh (k ) ** 4
@@ -1040,17 +978,12 @@ def nond_G(x, t):
1040978 + delta_rho * k ** 2
1041979 + k ** 2 * rhog
1042980 + sqrt (
1043- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
1044- * sinh (k ) ** 4
981+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
1045982 + delta_rho ** 2 * k ** 2
1046983 + 2 * delta_rho * k ** 2 * rhog
1047984 + k ** 2 * rhog ** 2
1048985 + 2
1049- * (
1050- delta_rho ** 2 * k
1051- + 2 * delta_rho * k * rhog
1052- + k * rhog ** 2
1053- )
986+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
1054987 * sinh (k )
1055988 * cosh (k )
1056989 + (
@@ -1106,17 +1039,12 @@ def nond_G(x, t):
11061039 + delta_rho * k ** 2
11071040 + k ** 2 * rhog
11081041 - sqrt (
1109- (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 )
1110- * sinh (k ) ** 4
1042+ (delta_rho ** 2 - 2 * delta_rho * rhog + rhog ** 2 ) * sinh (k ) ** 4
11111043 + delta_rho ** 2 * k ** 2
11121044 + 2 * delta_rho * k ** 2 * rhog
11131045 + k ** 2 * rhog ** 2
11141046 + 2
1115- * (
1116- delta_rho ** 2 * k
1117- + 2 * delta_rho * k * rhog
1118- + k * rhog ** 2
1119- )
1047+ * (delta_rho ** 2 * k + 2 * delta_rho * k * rhog + k * rhog ** 2 )
11201048 * sinh (k )
11211049 * cosh (k )
11221050 + (
0 commit comments