I have the following function: $$\frac{2e^x}{e^{2x}+1+2x}=\sum_{n=0}^\infty \varepsilon_n\frac{x^n}{n!}$$
I would like to find a closed form for the $\varepsilon_k$. One thing that I do know is that the $\varepsilon_n$ satisfy the following recurrence relation:
$$\varepsilon_n=1-2n\varepsilon_{n-1}-\sum_{k=0}^{n-2}\binom{n}{k}2^{n-k-1}\varepsilon_k$$
The first 11 numbers are below:
$$1,-1,3,-15,93,-725,6815, -74627, 933849, -13148361, 205690779$$
I've tried slowly plugging in small values of $n$, expanding the recurrence relation and then condensing it back to something smaller. I try to find patterns then amongst the summands and products but this method is very long and tedious, and the size of the summands expands quickly. What are some of the methods the have been used for closed forms similar to this and how can I better approach the problem?