I am trying to evaluate the integral

$$ \int_0^{2\pi}\frac{D_{11}^2+2D_{12}^2+D_{22}^2+D_{33}^2}{1+e\cos x(t)}dx $$

Where $D_{ij}$ are third derivatives wrt $t$ as defined below.

a[psi_] := k1/(1 + e*Cos[psi]);
d[psi_] := p1/(1 + e*Cos[psi]);
D11 = D[a[x[t]]*(Cos[x[t]]^2 - 1/3), {t, 3}];
D12 = D[a[x[t]]*Cos[x[t]]*Sin[x[t]], {t, 3}];
D22 = D[a[x[t]]*(Sin[x[t]]^2 - 1/3), {t, 3}];
D33 = D[-a[x[t]]/3];

To do this I have written the input

Assuming[Element[k1, Reals] && Element[e, Reals] && 
x'[t] == p2/(d[x[t]]^2) && Element[p2, Reals] && 
Element[p2, Reals], 
Integrate[(D11^2 + 2*D12^2 + D22^2 + D33^2)/(1 + e*Cos[x[t]])^2, {x[
t], 0, 2*Pi}]]

This is taking a very long time to evaluate, is it OK to specify to Mathematica that I have the condition $\dot{x}=\frac{p_2}{d(x(t))^2}$ in this manner? How could this be sped up?

EDIT: I have re-written this without subscripts.

  • 2
    $\begingroup$ To warn you at the outset: using subscripts for anything other than pretty formatting will give you nothing but inconvenience; avoid them altogether. $\endgroup$ – J. M. is slightly pensive Oct 11 '18 at 16:00
  • $\begingroup$ @J.M.issomewhatokay. in the sense that it is a pain to type, certainly. It won't do anything to my calculation though, right? $\endgroup$ – tbfr416 Oct 11 '18 at 16:12
  • 2
    $\begingroup$ Subscripts will harm your calculation because they appear to the human as a variable or symbol but to the kernel as a formatting construct. Just use D11 and D23 instead. $\endgroup$ – bill s Oct 11 '18 at 16:20

Your Answer

By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Browse other questions tagged or ask your own question.