As laid out on Wikipedia, the function $$f(x):=\begin{cases} \exp(-1/x) & x>0\\ 0 & x\le 0 \end{cases}$$
has the expression for derivatives at $x>0$, $$ f^{(n)}(x) = \frac{p_n(x)}{x^{2n}}f(x) $$ where the polynomials $p_n$ of degree $n-1$ are recursively defined by $p_1=1$, $$ p_{n+1} = (1-2nx)p_{n} + x^2 p_{n}'$$
I am interested in bounding $f^{(n)}$ on sets like $[0,a]$. Since its well known that all derivatives vanish at $x=0$, and $f(x)/x^{2n}$ never vanishes for $x>0$, the local extrema of $f^{(n)}$ occurs at a zero of $p_{n+1}$, so bounding $f^{(n)}(x)$ depends on value of its zeros.
Alternatively, perhaps a crude bound on $p_n$ on these sets $[0,a]$ is obtainable from the recurrence. This sounds easier, but I have yet to make progress here as well.
Question
How do I get a bound on $f^{(n)}$? (Preferrably one that looks like the inequality below)
Motivation
to test my understanding of Gevrey functions, I have decided to check what which space (if any) this function belongs to. This requires getting some bounds on the derivative. My suspicion is that $f$ is Gevrey of order 2, i.e. for any compact set $K$(and in particular any compact set containing $0$) there is a constant $C$ such that $$ \sup_{x\in K} |f^{(n)}(x)| \le C^{n+1} n!^2 $$
Notable "similar question"(currently unanswered)
Bounding $n$-th derivative of $x \mapsto \exp\left({-\frac{1}{x^2}}\right)$
There is a comment there that reaches somewhat the same idea (the associated recursive polynomial is important) but doesn't say any more.