Mathematica 8 отлично работает здесь.Предположим, что я интерпретирую Resi следующим образом (не знаю, правильно ли это, но возможно):
Resi = 1754.2*Sin[(Pi*x)/200]*A[1, 0][t] + 6.64531*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t] + 3.01996*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*
A[1, 0][t]*A[1, 7][t] - 861.531*Cos[7*\[Theta]]*Sin[(Pi*x)/200]^2*
A[1, 0][t]*A[1, 7][t] + 0.533983*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 0][t]^2*A[1, 7][t] - 0.369945*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 0][t]^2*A[1, 7][t] - 107.444*A[1, 7][t]^2 +
107.444*Cos[(Pi*x)/100]*A[1, 7][t]^2 - 3.95111*Cos[7*\[Theta]]^2*
Sin[(Pi*x)/200]^2*A[1, 7][t]^2 + 3.95111*Cos[(Pi*x)/200]^2*
Sin[7*\[Theta]]^2*A[1, 7][t]^2 + 0.533983*Cos[7*\[Theta]]^2*Sin[(Pi*x)/200]*
A[1, 0][t]*A[1, 7][t]^2 - 1.84972*Cos[(Pi*x)/100]*Cos[7*\[Theta]]^2*
Sin[(Pi*x)/200]*A[1, 0][t]*A[1, 7][t]^2 + 0.133496*Cos[14*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 0][t]*A[1, 7][t]^2 - 1.47978*Cos[(Pi*x)/200]*
Sin[(Pi*x)/100]*Sin[7*\[Theta]]^2*A[1, 0][t]*A[1, 7][t]^2 +
52.6477*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]^3 -
52.6477*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]^3 +
0.133496*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]^3 +
39.1352*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t] +
0.987777*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*A[1, 0][t]*A[1, 14][t] -
3438.75*Cos[14*\[Theta]]*Sin[(Pi*x)/200]^2*A[1, 0][t]*A[1, 14][t] +
0.533983*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]^2*A[1, 14][t] -
0.484011*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]^2*
A[1, 14][t] + 13.5898*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*A[1, 7][t]*
A[1, 14][t] + 0.0257425*Cos[(Pi*x)/100]*Cos[21*\[Theta]]*A[1, 7][t]*
A[1, 14][t] - 10.542*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]^2*
A[1, 7][t]*A[1, 14][t] + 8.43363*Cos[(Pi*x)/200]^2*Sin[7*\[Theta]]*
Sin[14*\[Theta]]*A[1, 7][t]*A[1, 14][t] + 0.266991*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 0][t]*A[1, 7][t]*A[1, 14][t] -
1.66475*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*
A[1, 7][t]*A[1, 14][t] + 1.06797*Cos[7*\[Theta]]*Cos[14*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 0][t]*A[1, 7][t]*A[1, 14][t] -
7.25708*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 0][t]*A[1, 7][t]*A[1, 14][t] + 0.266991*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 0][t]*A[1, 7][t]*A[1, 14][t] -
0.0283811*Cos[(Pi*x)/100]*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*
A[1, 7][t]*A[1, 14][t] - 3.92758*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*
Sin[7*\[Theta]]*Sin[14*\[Theta]]*A[1, 0][t]*A[1, 7][t]*A[1, 14][t] +
0.266991*Cos[7*\[Theta]]^2*Sin[(Pi*x)/200]*A[1, 7][t]^2*A[1, 14][t] -
8.32376*Cos[(Pi*x)/100]*Cos[7*\[Theta]]^2*Sin[(Pi*x)/200]*A[1, 7][t]^2*
A[1, 14][t] + 210.591*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]^2*
A[1, 14][t] - 210.591*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 7][t]^2*A[1, 14][t] + 0.133496*Cos[14*\[Theta]]^2*Sin[(Pi*x)/200]*
A[1, 7][t]^2*A[1, 14][t] + 0.266991*Cos[7*\[Theta]]*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t]^2*A[1, 14][t] -
0.040995*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 7][t]^2*A[1, 14][t] - 6.65901*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*
Sin[7*\[Theta]]^2*A[1, 7][t]^2*A[1, 14][t] - 0.0378415*Cos[(Pi*x)/200]*
Sin[(Pi*x)/100]*Sin[7*\[Theta]]*Sin[21*\[Theta]]*A[1, 7][t]^2*A[1, 14][t] -
429.778*A[1, 14][t]^2 + 429.778*Cos[(Pi*x)/100]*A[1, 14][t]^2 -
1.06283*Cos[14*\[Theta]]^2*Sin[(Pi*x)/200]^2*A[1, 14][t]^2 +
1.06283*Cos[(Pi*x)/200]^2*Sin[14*\[Theta]]^2*A[1, 14][t]^2 +
0.533983*Cos[14*\[Theta]]^2*Sin[(Pi*x)/200]*A[1, 0][t]*A[1, 14][t]^2 -
2.42005*Cos[(Pi*x)/100]*Cos[14*\[Theta]]^2*Sin[(Pi*x)/200]*A[1, 0][t]*
A[1, 14][t]^2 + 0.133496*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*
A[1, 14][t]^2 - 1.93604*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*
Sin[14*\[Theta]]^2*A[1, 0][t]*A[1, 14][t]^2 + 210.591*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t]*A[1, 14][t]^2 -
210.591*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*
A[1, 14][t]^2 + 0.266991*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 7][t]*A[1, 14][t]^2 - 28.3008*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*
Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*A[1, 14][t]^2 +
0.266991*Cos[14*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*
A[1, 14][t]^2 - 0.0788365*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t]*A[1, 14][t]^2 +
0.133496*Cos[7*\[Theta]]*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*
A[1, 14][t]^2 - 13.318*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*Sin[7*\[Theta]]*
Sin[14*\[Theta]]*A[1, 7][t]*A[1, 14][t]^2 - 0.075683*Cos[(Pi*x)/200]*
Sin[(Pi*x)/100]*Sin[14*\[Theta]]*Sin[21*\[Theta]]*A[1, 7][t]*A[1, 14][t]^2 +
842.364*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]^3 -
842.364*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]^3 +
0.133496*Cos[14*\[Theta]]*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]^3 +
6.64531*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t] +
3.01996*Cos[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 0][t]*B[1, 7][t] -
861.531*Sin[(Pi*x)/200]^2*Sin[7*\[Theta]]*A[1, 0][t]*B[1, 7][t] +
0.533983*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 0][t]^2*B[1, 7][t] -
0.369945*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 0][t]^2*
B[1, 7][t] - 7.90221*Cos[(Pi*x)/200]^2*Cos[7*\[Theta]]*Sin[7*\[Theta]]*
A[1, 7][t]*B[1, 7][t] - 7.90221*Cos[7*\[Theta]]*Sin[(Pi*x)/200]^2*
Sin[7*\[Theta]]*A[1, 7][t]*B[1, 7][t] + 1.06797*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 0][t]*A[1, 7][t]*B[1, 7][t] -
3.69945*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
A[1, 0][t]*A[1, 7][t]*B[1, 7][t] + 2.95956*Cos[(Pi*x)/200]*
Cos[7*\[Theta]]*Sin[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 0][t]*A[1, 7][t]*
B[1, 7][t] + 0.266991*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 0][t]*
A[1, 7][t]*B[1, 7][t] + 52.6477*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
A[1, 7][t]^2*B[1, 7][t] - 52.6477*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*A[1, 7][t]^2*B[1, 7][t] + 0.133496*Cos[14*\[Theta]]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 7][t]^2*B[1, 7][t] +
0.266991*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]^2*
B[1, 7][t] - 13.5898*Cos[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 14][t]*
B[1, 7][t] - 10.542*Cos[14*\[Theta]]*Sin[(Pi*x)/200]^2*Sin[7*\[Theta]]*
A[1, 14][t]*B[1, 7][t] - 8.43363*Cos[(Pi*x)/200]^2*Cos[7*\[Theta]]*
Sin[14*\[Theta]]*A[1, 14][t]*B[1, 7][t] + 0.0257425*Cos[(Pi*x)/100]*
Sin[21*\[Theta]]*A[1, 14][t]*B[1, 7][t] - 0.266991*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*A[1, 0][t]*A[1, 14][t]*B[1, 7][t] +
1.66475*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 0][t]*
A[1, 14][t]*B[1, 7][t] + 1.06797*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*A[1, 0][t]*A[1, 14][t]*B[1, 7][t] -
7.25708*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
A[1, 0][t]*A[1, 14][t]*B[1, 7][t] + 3.92758*Cos[(Pi*x)/200]*
Cos[7*\[Theta]]*Sin[(Pi*x)/100]*Sin[14*\[Theta]]*A[1, 0][t]*A[1, 14][t]*
B[1, 7][t] + 0.266991*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*A[1, 0][t]*
A[1, 14][t]*B[1, 7][t] - 0.0283811*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[21*\[Theta]]*A[1, 0][t]*A[1, 14][t]*B[1, 7][t] +
0.*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 7][t]*A[1, 14][t]*
B[1, 7][t] + 0.*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 7][t] +
0.266991*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 7][t]*
A[1, 14][t]*B[1, 7][t] - 0.040995*Cos[(Pi*x)/100]*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 7][t] +
0.*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Sin[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 7][t]*
A[1, 14][t]*B[1, 7][t] + 0.0378415*Cos[(Pi*x)/200]*Cos[21*\[Theta]]*
Sin[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 7][t] +
0.266991*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]*
A[1, 14][t]*B[1, 7][t] + 0.266991*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*
Sin[21*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 7][t] -
0.040995*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*
A[1, 7][t]*A[1, 14][t]*B[1, 7][t] + 0.0378415*Cos[(Pi*x)/200]*
Cos[7*\[Theta]]*Sin[(Pi*x)/100]*Sin[21*\[Theta]]*A[1, 7][t]*A[1, 14][t]*
B[1, 7][t] + 210.591*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 14][t]^2*
B[1, 7][t] - 210.591*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
A[1, 14][t]^2*B[1, 7][t] - 0.266991*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*A[1, 14][t]^2*B[1, 7][t] + 28.3008*Cos[(Pi*x)/100]*
Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 14][t]^2*B[1, 7][t] +
0.133496*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 14][t]^2*
B[1, 7][t] - 13.318*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Sin[(Pi*x)/100]*
Sin[14*\[Theta]]*A[1, 14][t]^2*B[1, 7][t] + 0.075683*Cos[(Pi*x)/200]*
Cos[21*\[Theta]]*Sin[(Pi*x)/100]*Sin[14*\[Theta]]*A[1, 14][t]^2*B[1, 7][t] +
0.266991*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*A[1, 14][t]^2*
B[1, 7][t] - 0.0788365*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
Sin[21*\[Theta]]*A[1, 14][t]^2*B[1, 7][t] - 107.444*B[1, 7][t]^2 +
107.444*Cos[(Pi*x)/100]*B[1, 7][t]^2 + 3.95111*Cos[(Pi*x)/200]^2*
Cos[7*\[Theta]]^2*B[1, 7][t]^2 - 3.95111*Sin[(Pi*x)/200]^2*Sin[7*\[Theta]]^2*
B[1, 7][t]^2 - 0.133496*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*
B[1, 7][t]^2 - 1.47978*Cos[(Pi*x)/200]*Cos[7*\[Theta]]^2*Sin[(Pi*x)/100]*
A[1, 0][t]*B[1, 7][t]^2 + 0.533983*Sin[(Pi*x)/200]*Sin[7*\[Theta]]^2*
A[1, 0][t]*B[1, 7][t]^2 - 1.84972*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]^2*A[1, 0][t]*B[1, 7][t]^2 + 52.6477*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t]*B[1, 7][t]^2 - 52.6477*Cos[(Pi*x)/100]*
Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*B[1, 7][t]^2 -
0.133496*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*
B[1, 7][t]^2 + 0.266991*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[14*\[Theta]]*
A[1, 7][t]*B[1, 7][t]^2 + 210.591*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 14][t]*B[1, 7][t]^2 - 210.591*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 14][t]*B[1, 7][t]^2 -
0.133496*Cos[14*\[Theta]]^2*Sin[(Pi*x)/200]*A[1, 14][t]*B[1, 7][t]^2 +
6.65901*Cos[(Pi*x)/200]*Cos[7*\[Theta]]^2*Sin[(Pi*x)/100]*A[1, 14][t]*
B[1, 7][t]^2 - 0.0378415*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Cos[21*\[Theta]]*
Sin[(Pi*x)/100]*A[1, 14][t]*B[1, 7][t]^2 -
0.266991*Sin[(Pi*x)/200]*Sin[7*\[Theta]]^2*A[1, 14][t]*B[1, 7][t]^2 +
8.32376*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]^2*A[1, 14][t]*
B[1, 7][t]^2 + 0.266991*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[21*\[Theta]]*
A[1, 14][t]*B[1, 7][t]^2 - 0.040995*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*Sin[21*\[Theta]]*A[1, 14][t]*B[1, 7][t]^2 +
52.6477*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]^3 -
52.6477*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]^3 -
0.133496*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]^3 +
39.1352*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*B[1, 14][t] +
0.987777*Cos[(Pi*x)/100]*Sin[14*\[Theta]]*A[1, 0][t]*B[1, 14][t] -
3438.75*Sin[(Pi*x)/200]^2*Sin[14*\[Theta]]*A[1, 0][t]*B[1, 14][t] +
0.533983*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 0][t]^2*B[1, 14][t] -
0.484011*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 0][t]^2*
B[1, 14][t] + 13.5898*Cos[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 7][t]*
B[1, 14][t] - 8.43363*Cos[(Pi*x)/200]^2*Cos[14*\[Theta]]*Sin[7*\[Theta]]*
A[1, 7][t]*B[1, 14][t] - 10.542*Cos[7*\[Theta]]*Sin[(Pi*x)/200]^2*
Sin[14*\[Theta]]*A[1, 7][t]*B[1, 14][t] + 0.0257425*Cos[(Pi*x)/100]*
Sin[21*\[Theta]]*A[1, 7][t]*B[1, 14][t] + 0.266991*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*A[1, 0][t]*A[1, 7][t]*B[1, 14][t] -
1.66475*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 0][t]*
A[1, 7][t]*B[1, 14][t] + 3.92758*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*
Sin[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 0][t]*A[1, 7][t]*B[1, 14][t] +
1.06797*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 0][t]*A[1, 7][t]*
B[1, 14][t] - 7.25708*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]*A[1, 0][t]*A[1, 7][t]*B[1, 14][t] +
0.266991*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*A[1, 0][t]*A[1, 7][t]*
B[1, 14][t] - 0.0283811*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*
A[1, 0][t]*A[1, 7][t]*B[1, 14][t] + 0.266991*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 7][t]^2*B[1, 14][t] -
8.32376*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
A[1, 7][t]^2*B[1, 14][t] + 6.65901*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*
Sin[(Pi*x)/100]*Sin[7*\[Theta]]*A[1, 7][t]^2*B[1, 14][t] +
0.0378415*Cos[(Pi*x)/200]*Cos[21*\[Theta]]*Sin[(Pi*x)/100]*Sin[7*\[Theta]]*
A[1, 7][t]^2*B[1, 14][t] + 210.591*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
A[1, 7][t]^2*B[1, 14][t] - 210.591*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]*A[1, 7][t]^2*B[1, 14][t] + 0.133496*Cos[14*\[Theta]]*
Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]^2*B[1, 14][t] +
0.266991*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*A[1, 7][t]^2*
B[1, 14][t] - 0.040995*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*
Sin[21*\[Theta]]*A[1, 7][t]^2*B[1, 14][t] - 2.12567*Cos[(Pi*x)/200]^2*
Cos[14*\[Theta]]*Sin[14*\[Theta]]*A[1, 14][t]*B[1, 14][t] -
2.12567*Cos[14*\[Theta]]*Sin[(Pi*x)/200]^2*Sin[14*\[Theta]]*A[1, 14][t]*
B[1, 14][t] + 1.06797*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
A[1, 0][t]*A[1, 14][t]*B[1, 14][t] - 4.84011*Cos[(Pi*x)/100]*
Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 0][t]*A[1, 14][t]*
B[1, 14][t] + 3.87208*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*
Sin[14*\[Theta]]*A[1, 0][t]*A[1, 14][t]*B[1, 14][t] +
0.266991*Sin[(Pi*x)/200]*Sin[28*\[Theta]]*A[1, 0][t]*A[1, 14][t]*
B[1, 14][t] + 0.266991*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
A[1, 7][t]*A[1, 14][t]*B[1, 14][t] - 28.3008*Cos[(Pi*x)/100]*
Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*A[1, 7][t]*A[1, 14][t]*
B[1, 14][t] + 13.318*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*
Sin[7*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 14][t] +
0.266991*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]*
A[1, 14][t]*B[1, 14][t] - 28.3008*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 14][t] +
0.266991*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]*
A[1, 14][t]*B[1, 14][t] - 0.0788365*Cos[(Pi*x)/100]*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 14][t] +
13.318*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Sin[(Pi*x)/100]*Sin[14*\[Theta]]*
A[1, 7][t]*A[1, 14][t]*B[1, 14][t] + 0.075683*Cos[(Pi*x)/200]*
Cos[21*\[Theta]]*Sin[(Pi*x)/100]*Sin[14*\[Theta]]*A[1, 7][t]*A[1, 14][t]*
B[1, 14][t] + 0.266991*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*
A[1, 7][t]*A[1, 14][t]*B[1, 14][t] - 0.0788365*Cos[(Pi*x)/100]*
Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[21*\[Theta]]*A[1, 7][t]*A[1, 14][t]*
B[1, 14][t] + 0.075683*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*
Sin[21*\[Theta]]*A[1, 7][t]*A[1, 14][t]*B[1, 14][t] +
0.266991*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[28*\[Theta]]*A[1, 7][t]*
A[1, 14][t]*B[1, 14][t] + 842.364*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
A[1, 14][t]^2*B[1, 14][t] - 842.364*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]*A[1, 14][t]^2*B[1, 14][t] + 0.133496*Cos[28*\[Theta]]*
Sin[(Pi*x)/200]*Sin[14*\[Theta]]*A[1, 14][t]^2*B[1, 14][t] +
0.266991*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[28*\[Theta]]*A[1, 14][t]^2*
B[1, 14][t] + 13.5898*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*B[1, 7][t]*
B[1, 14][t] + 8.43363*Cos[(Pi*x)/200]^2*Cos[7*\[Theta]]*Cos[14*\[Theta]]*
B[1, 7][t]*B[1, 14][t] - 0.0257425*Cos[(Pi*x)/100]*Cos[21*\[Theta]]*
B[1, 7][t]*B[1, 14][t] - 10.542*Sin[(Pi*x)/200]^2*Sin[7*\[Theta]]*
Sin[14*\[Theta]]*B[1, 7][t]*B[1, 14][t] + 0.266991*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 0][t]*B[1, 7][t]*B[1, 14][t] -
1.66475*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*
B[1, 7][t]*B[1, 14][t] - 0.266991*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 0][t]*B[1, 7][t]*B[1, 14][t] + 0.0283811*Cos[(Pi*x)/100]*
Cos[21*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*B[1, 7][t]*B[1, 14][t] -
3.92758*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*
A[1, 0][t]*B[1, 7][t]*B[1, 14][t] + 1.06797*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*Sin[14*\[Theta]]*A[1, 0][t]*B[1, 7][t]*B[1, 14][t] -
7.25708*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[14*\[Theta]]*
A[1, 0][t]*B[1, 7][t]*B[1, 14][t] + 0.266991*Cos[7*\[Theta]]^2*
Sin[(Pi*x)/200]*A[1, 7][t]*B[1, 7][t]*B[1, 14][t] -
8.32376*Cos[(Pi*x)/100]*Cos[7*\[Theta]]^2*Sin[(Pi*x)/200]*A[1, 7][t]*
B[1, 7][t]*B[1, 14][t] - 0.266991*Cos[7*\[Theta]]*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t]*B[1, 7][t]*B[1, 14][t] +
0.040995*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 7][t]*B[1, 7][t]*B[1, 14][t] - 6.65901*Cos[(Pi*x)/200]*
Cos[7*\[Theta]]^2*Sin[(Pi*x)/100]*A[1, 7][t]*B[1, 7][t]*B[1, 14][t] -
0.0378415*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/100]*
A[1, 7][t]*B[1, 7][t]*B[1, 14][t] + 0.266991*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]^2*A[1, 7][t]*B[1, 7][t]*B[1, 14][t] -
8.32376*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]^2*A[1, 7][t]*
B[1, 7][t]*B[1, 14][t] - 6.65901*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*
Sin[7*\[Theta]]^2*A[1, 7][t]*B[1, 7][t]*B[1, 14][t] +
0.266991*Sin[(Pi*x)/200]*Sin[14*\[Theta]]^2*A[1, 7][t]*B[1, 7][t]*
B[1, 14][t] + 0.266991*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[21*\[Theta]]*
A[1, 7][t]*B[1, 7][t]*B[1, 14][t] - 0.040995*Cos[(Pi*x)/100]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[21*\[Theta]]*A[1, 7][t]*B[1, 7][t]*
B[1, 14][t] + 0.0378415*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*Sin[7*\[Theta]]*
Sin[21*\[Theta]]*A[1, 7][t]*B[1, 7][t]*B[1, 14][t] +
0.266991*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]*
B[1, 7][t]*B[1, 14][t] - 28.3008*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*
Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]*B[1, 7][t]*B[1, 14][t] -
0.266991*Cos[14*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]*
B[1, 7][t]*B[1, 14][t] + 0.0788365*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*
Cos[21*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]*B[1, 7][t]*B[1, 14][t] +
13.318*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*
A[1, 14][t]*B[1, 7][t]*B[1, 14][t] - 0.075683*Cos[(Pi*x)/200]*
Cos[14*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/100]*A[1, 14][t]*B[1, 7][t]*
B[1, 14][t] - 0.266991*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[14*\[Theta]]*
A[1, 14][t]*B[1, 7][t]*B[1, 14][t] + 28.3008*Cos[(Pi*x)/100]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[14*\[Theta]]*A[1, 14][t]*B[1, 7][t]*
B[1, 14][t] - 13.318*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*Sin[7*\[Theta]]*
Sin[14*\[Theta]]*A[1, 14][t]*B[1, 7][t]*B[1, 14][t] +
0.266991*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*Sin[21*\[Theta]]*A[1, 14][t]*
B[1, 7][t]*B[1, 14][t] - 0.0788365*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]*Sin[21*\[Theta]]*A[1, 14][t]*B[1, 7][t]*B[1, 14][t] +
0.075683*Cos[(Pi*x)/200]*Sin[(Pi*x)/100]*Sin[14*\[Theta]]*Sin[21*\[Theta]]*
A[1, 14][t]*B[1, 7][t]*B[1, 14][t] + 0.266991*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*Sin[28*\[Theta]]*A[1, 14][t]*B[1, 7][t]*B[1, 14][t] +
0.266991*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]^2*
B[1, 14][t] - 8.32376*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*B[1, 7][t]^2*B[1, 14][t] - 0.266991*Cos[21*\[Theta]]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]^2*B[1, 14][t] +
0.040995*Cos[(Pi*x)/100]*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
B[1, 7][t]^2*B[1, 14][t] + 6.65901*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*
Sin[(Pi*x)/100]*Sin[7*\[Theta]]*B[1, 7][t]^2*B[1, 14][t] +
210.591*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*B[1, 7][t]^2*B[1, 14][t] -
210.591*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*B[1, 7][t]^2*
B[1, 14][t] - 0.133496*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
B[1, 7][t]^2*B[1, 14][t] - 0.0378415*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*
Sin[(Pi*x)/100]*Sin[21*\[Theta]]*B[1, 7][t]^2*B[1, 14][t] -
429.778*B[1, 14][t]^2 + 429.778*Cos[(Pi*x)/100]*B[1, 14][t]^2 +
1.06283*Cos[(Pi*x)/200]^2*Cos[14*\[Theta]]^2*B[1, 14][t]^2 -
1.06283*Sin[(Pi*x)/200]^2*Sin[14*\[Theta]]^2*B[1, 14][t]^2 -
0.133496*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*A[1, 0][t]*B[1, 14][t]^2 -
1.93604*Cos[(Pi*x)/200]*Cos[14*\[Theta]]^2*Sin[(Pi*x)/100]*A[1, 0][t]*
B[1, 14][t]^2 + 0.533983*Sin[(Pi*x)/200]*Sin[14*\[Theta]]^2*A[1, 0][t]*
B[1, 14][t]^2 - 2.42005*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]^2*A[1, 0][t]*B[1, 14][t]^2 + 210.591*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*A[1, 7][t]*B[1, 14][t]^2 -
210.591*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*A[1, 7][t]*
B[1, 14][t]^2 - 0.133496*Cos[7*\[Theta]]*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 7][t]*B[1, 14][t]^2 - 13.318*Cos[(Pi*x)/200]*Cos[7*\[Theta]]*
Cos[14*\[Theta]]*Sin[(Pi*x)/100]*A[1, 7][t]*B[1, 14][t]^2 -
0.075683*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*Cos[21*\[Theta]]*Sin[(Pi*x)/100]*
A[1, 7][t]*B[1, 14][t]^2 + 0.266991*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*
Sin[14*\[Theta]]*A[1, 7][t]*B[1, 14][t]^2 - 28.3008*Cos[(Pi*x)/100]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*Sin[14*\[Theta]]*A[1, 7][t]*B[1, 14][t]^2 +
0.266991*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*Sin[21*\[Theta]]*A[1, 7][t]*
B[1, 14][t]^2 - 0.0788365*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]*Sin[21*\[Theta]]*A[1, 7][t]*B[1, 14][t]^2 +
842.364*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]*B[1, 14][t]^2 -
842.364*Cos[(Pi*x)/100]*Cos[14*\[Theta]]*Sin[(Pi*x)/200]*A[1, 14][t]*
B[1, 14][t]^2 - 0.133496*Cos[14*\[Theta]]*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*
A[1, 14][t]*B[1, 14][t]^2 + 0.266991*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
Sin[28*\[Theta]]*A[1, 14][t]*B[1, 14][t]^2 + 210.591*Sin[(Pi*x)/200]*
Sin[7*\[Theta]]*B[1, 7][t]*B[1, 14][t]^2 - 210.591*Cos[(Pi*x)/100]*
Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]*B[1, 14][t]^2 -
0.133496*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*Sin[7*\[Theta]]*B[1, 7][t]*
B[1, 14][t]^2 + 13.318*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*
Sin[7*\[Theta]]*B[1, 7][t]*B[1, 14][t]^2 + 0.266991*Cos[7*\[Theta]]*
Sin[(Pi*x)/200]*Sin[14*\[Theta]]*B[1, 7][t]*B[1, 14][t]^2 -
28.3008*Cos[(Pi*x)/100]*Cos[7*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
B[1, 7][t]*B[1, 14][t]^2 - 0.266991*Cos[21*\[Theta]]*Sin[(Pi*x)/200]*
Sin[14*\[Theta]]*B[1, 7][t]*B[1, 14][t]^2 + 0.0788365*Cos[(Pi*x)/100]*
Cos[21*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*B[1, 7][t]*B[1, 14][t]^2 -
0.075683*Cos[(Pi*x)/200]*Cos[14*\[Theta]]*Sin[(Pi*x)/100]*Sin[21*\[Theta]]*
B[1, 7][t]*B[1, 14][t]^2 + 842.364*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
B[1, 14][t]^3 - 842.364*Cos[(Pi*x)/100]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
B[1, 14][t]^3 - 0.133496*Cos[28*\[Theta]]*Sin[(Pi*x)/200]*Sin[14*\[Theta]]*
B[1, 14][t]^3 + ((6.90612*Sin[(Pi*x)/200])/10^7)*
Derivative[2][A[1, 0]][t] + ((6.90612*Cos[7*\[Theta]]*Sin[(Pi*x)/200])/
10^7)*Derivative[2][A[1, 7]][t] +
((6.90612*Cos[14*\[Theta]]*Sin[(Pi*x)/200])/10^7)*Derivative[2][A[1, 14]][
t] + ((6.90612*Sin[(Pi*x)/200]*Sin[7*\[Theta]])/10^7)*
Derivative[2][B[1, 7]][t] + ((6.90612*Sin[(Pi*x)/200]*Sin[14*\[Theta]])/
10^7)*Derivative[2][B[1, 14]][t];
Тогда:
With[{L = 200},
Print[AbsoluteTiming[r1 = ExpandAll[Chop[Integrate[Resi*Cos[7*1*\[Theta]]*Sin[(1*Pi*x)/L], {\[Theta], 0, 2*Pi},
{x, 0, L}], 10^(-6)]]]];
Print[AbsoluteTiming[r2 = ExpandAll[Chop[Integrate[Resi*Cos[7*1*\[Theta]]*Sin[(1*Pi*x)/L], {x, 0, L},
{\[Theta], 0, 2*Pi}], 10^(-6)]]]]; r1 === r2]
дает, подтверждая йоду,:
{36.7641028,
2087.69 A[1, 7][t] - 230144. A[1, 0][t] A[1, 7][t] +
225.867 A[1, 0][t]^2 A[1, 7][t] + 24830.6 A[1, 7][t]^3 -
2655.33 A[1, 7][t] A[1, 14][t] +
774.631 A[1, 0][t] A[1, 7][t] A[1, 14][t] +
100500. A[1, 7][t] A[1, 14][t]^2 +
24830.6 A[1, 7][t] B[1, 7][t]^2 - 2655.33 B[1, 7][t] B[1, 14][t] +
774.631 A[1, 0][t] B[1, 7][t] B[1, 14][t] +
100500. A[1, 7][t] B[1, 14][t]^2 +
0.000216962 (A[1, 7]^\[Prime]\[Prime])[t]}
{23.8343632,
2087.69 A[1, 7][t] - 230144. A[1, 0][t] A[1, 7][t] +
225.867 A[1, 0][t]^2 A[1, 7][t] + 24830.6 A[1, 7][t]^3 -
2655.33 A[1, 7][t] A[1, 14][t] +
774.631 A[1, 0][t] A[1, 7][t] A[1, 14][t] +
100500. A[1, 7][t] A[1, 14][t]^2 +
24830.6 A[1, 7][t] B[1, 7][t]^2 - 2655.33 B[1, 7][t] B[1, 14][t] +
774.631 A[1, 0][t] B[1, 7][t] B[1, 14][t] +
100500. A[1, 7][t] B[1, 14][t]^2 +
0.000216962 (A[1, 7]^\[Prime]\[Prime])[t]}
Конечно, это действительно выглядит слишком медленно.Поэтому, если вы поможете Интеграции, немного переставив интегральную функцию, вы сможете значительно ускорить процесс:
With[{L = 200},
Print[AbsoluteTiming[rint0 = Expand[Collect[Expand[Resi*Cos[7*\[Theta]]*Sin[(Pi*x)/L], x] /.
p_Times /; !FreeQ[p, x] :> Hold[Integrate][Select[p, !FreeQ[#1, x] & ], {x, 0, L}]*
Select[p, FreeQ[#1, x] & ], Hold[Integrate][_, _]] /. Hold[Integrate] -> Integrate, \[Theta]];
rint1 = Chop[Expand[Collect[rint0 /. p_Times /; !FreeQ[p, \[Theta]] :> Hold[Integrate][Select[p,
!FreeQ[#1, \[Theta]] & ], {\[Theta], 0, 2*Pi}]*Select[p, FreeQ[#1, \[Theta]] & ],
Hold[Integrate][_, _]] /. Hold[Integrate] -> Integrate], 1/10^6]]];
Print[AbsoluteTiming[tint1 = Expand[Collect[Expand[Resi*Cos[7*\[Theta]]*Sin[(Pi*x)/L], \[Theta]] /.
p_Times /; !FreeQ[p, \[Theta]] :> Hold[Integrate][Select[p, !FreeQ[#1, \[Theta]] & ], {\[Theta], 0, 2*Pi}]*
Select[p, FreeQ[#1, \[Theta]] & ], Hold[Integrate][_, _]] /. Hold[Integrate] -> Integrate];
tint2 = Chop[Expand[Collect[tint1 /. p_Times /; !FreeQ[p, x] :> Hold[Integrate][Select[p,
!FreeQ[#1, x] & ], {x, 0, L}]*Select[p, FreeQ[#1, x] & ], Hold[Integrate][_, _]] /.
Hold[Integrate] -> Integrate], 1/10^6]]]; Chop[rint1 - tint2, 1/10^6]]
В результате
{1.5800904,2087.69 A[1,7][t]-230144. A[1,0][t] A[1,7][t]+225.867 A[1,0][t]^2 A[1,7][t]+24830.6 A[1,7][t]^3-2655.33 A[1,7][t] A[1,14][t]+774.631 A[1,0][t] A[1,7][t] A[1,14][t]+100500. A[1,7][t] A[1,14][t]^2+24830.6 A[1,7][t] B[1,7][t]^2-2655.33 B[1,7][t] B[1,14][t]+774.631 A[1,0][t] B[1,7][t] B[1,14][t]+100500. A[1,7][t] B[1,14][t]^2+0.000216962 (A[1,7]^\[Prime]\[Prime])[t]}
{1.5820904,2087.69 A[1,7][t]-230144. A[1,0][t] A[1,7][t]+225.867 A[1,0][t]^2 A[1,7][t]+24830.6 A[1,7][t]^3-2655.33 A[1,7][t] A[1,14][t]+774.631 A[1,0][t] A[1,7][t] A[1,14][t]+100500. A[1,7][t] A[1,14][t]^2+24830.6 A[1,7][t] B[1,7][t]^2-2655.33 B[1,7][t] B[1,14][t]+774.631 A[1,0][t] B[1,7][t] B[1,14][t]+100500. A[1,7][t] B[1,14][t]^2+0.000216962 (A[1,7]^\[Prime]\[Prime])[t]}