@RyderRude yes, but derivation the article is doing is for lagrangians that could be anything, so they still use only the first order. because e is infinetisemal, approximation of only first order will be pretty good, but it won't be 100% perfect which would make the later calculations incorrect.
$a \frac{\partial L}{\partial a} = \frac{d\Lambda}{dt}$ (option 1)
$a \frac{\partial L}{\partial a} + a^2 \frac{\partial^2 L}{\partial a^2} = \frac{d\Lambda}{dt}$ (option 2)
for all lagrangians, option 1 might not give the exact calculations which would cause the conservation law to be derived in…