I calculated this integral:
$${\int_{0}^{\infty}e^{-cx}x^{n}\arctan(ax)\mathrm{d}x=\frac{n!}{c^{n+1}}\left(\text{Ci}(b)\sum_{k=0}^{n}\sin\left(b-\frac{k\pi}{2}\right)\frac{b^k}{k!}+\left(\frac{\pi}{2}-\text{Si}(b)\right)\sum_{k=0}^{n}\cos\left(b-\frac{k\pi}{2}\right)\frac{b^k}{k!}+\sum_{k=1}^{n}\frac{1}{k}\sum_{j=0}^{k-1}\sum_{l=0}^{n-k}\frac{\sin\left(\pi\cdot\frac{l-j}{2}\right)}{j!l!}b^{j+l}\right)}$$
Where $b:=\dfrac{c}{a}$
I wanted to ask for some help with manipulating the last summations
$$\sum_{k=1}^{n}\frac{1}{k}\sum_{j=0}^{k-1}\sum_{l=0}^{n-k}\frac{\sin\left(\pi\cdot\frac{l-j}{2}\right)}{j!l!}b^{j+l}$$
since the only values that sine can have are $1$, $-1$ and $0$ I imagine that at least one of the summations can be removed.
In general I noticed that:
$$\sin\left(\frac{l-j}{2}\pi\right)=\begin{cases}
0&\text{if }j-l=2n_0\\
+1&\text{if }j-l=4n_0+1\\
-1&\text{if }j-l=4n_0-1\end{cases}$$
Edit
I noticed that the polynomial is odd, so technically it is possible to substitute $j+l=2s-1$ where $1\leq s\leq \left\lfloor\dfrac{n}{2}\right\rfloor$?
Thanks in advance for any suggestions