3 edited body

I don't think you'll get a closed form but you get get an expression that involves products and sums that isn't too horrible and is easy to evaluate for particular cases. It'll take me forever to type it up here with nice formatting but I'll describe the approach in a way that's easy to reproduce.

We're starting with the number 1 and repeatedly performing these two operations: integration and multiplication by $e^{nat}$ for various $n$. Switch to working with Laplace transforms. In Laplace transform space, integration (starting from zero) is simply multiplication by $1/s$ and multiplication by the exponential takes $f(s)$ to $f(s-na)$.

So the Laplace transform of $I_1$ is $1/(s(s-a))$, the Laplace transform of $I_2$ is $1/(s(s-3a)(s-4a))$ and so on. Notice how when we keep going, the denominator is made up of factors of the form $s-(r^2-i^2)a$. So it's easy to write the Laplace transform of $I_r$ using product notation. You now simply need to invert the Laplace transform. To do this, we need to write our rational function in $s$ as partial fractions. The factors in the denominator are all distinct so this is the easy case. After the inverse transform you should end up with something like

$e^{r^2t}\sum_{i=0}^r\exp(-i^2t)/\prod_{j\ne r}(i^2-j^2)$i}(i^2-j^2)$I expect that I've made a slip somewhere as I only have a few moments spare. That's why I said 'something like'. But the method is sound and if you put in more work than me you'll get a definite result. Update: Did some testing in Mathematica. I think the expression above is correct. 2 added 1 characters in body; edited body I don't think you'll get a closed form but you get get an expression that involves products and sums that isn't too horrible and is easy to evaluate for particular cases. It'll take me forever to type it up here with nice formatting but I'll describe the approach in a way that's easy to reproduce. We're staring starting with the number 1 and repeatedly performing these two operations: integration and multiplication by$e^{nat}$for various$n$. Switch to working with Laplace transforms. In Laplace transform space, integration (starting from zero) is simply multiplication by$1/s$and multiplication by the exponential takes$f(s)$to$f(s-na)$. So the Laplace transform of$I_1$is$1/(s(s-a))$, the Laplace transform of$I_2$is$1/(s(s-3a)(s-4a))$and so on. Notice how when we keep going, the denominator is made up of factors of the form$s-(r^2-i^2)a$. So it's easy to write the Laplace transform of$I_r$using product notation. You now simply need to invert the Laplace transform. To do this, we need to write our rational function in$s$as partial fractions. The factors in the denominator are all distinct so this is the easy case. After the inverse transform you should end up with something like$e^{rt^2}\sum_{i=0}^r\exp(-i^2t)/\prod_{j\ne e^{r^2t}\sum_{i=0}^r\exp(-i^2t)/\prod_{j\ne r}(i^2-j^2)$I expect that I've made a slip somewhere as I only have a few moments spare. That's why I said 'something like'. But the method is sound and if you put in more work than me you'll get a definite result. Update: Did some testing in Mathematica. I think the expression above is correct. 1 I don't think you'll get a closed form but you get get an expression that involves products and sums that isn't too horrible and is easy to evaluate for particular cases. It'll take me forever to type it up here with nice formatting but I'll describe the approach in a way that's easy to reproduce. We're staring with the number 1 and repeatedly performing these two operations: integration and multiplication by$e^{nat}$for various$n$. Switch to working with Laplace transforms. In Laplace transform space, integration (starting from zero) is simply multiplication by$1/s$and multiplication by the exponential takes$f(s)$to$f(s-na)$. So the Laplace transform of$I_1$is$1/(s(s-a))$, the Laplace transform of$I_2$is$1/(s(s-3a)(s-4a))$and so on. Notice how when we keep going, the denominator is made up of factors of the form$s-(r^2-i^2)a$. So it's easy to write the Laplace transform of$I_r$using product notation. You now simply need to invert the Laplace transform. To do this, we need to write our rational function in$s$as partial fractions. The factors in the denominator are all distinct so this is the easy case. After the inverse transform you should end up with something like$e^{rt^2}\sum_{i=0}^r\exp(-i^2t)/\prod_{j\ne r}(i^2-j^2)\$

I expect that I've made a slip somewhere as I only have a few moments spare. That's why I said 'something like'. But the method is sound and if you put in more work than me you'll get a definite result.

Update: Did some testing in Mathematica. I think the expression above is correct.