I have exponential terms of temperature in reaction rates used in parabolic heat equation. As I apply Galerkin method to find its ODE, it contains temperature in exponential form at right hand side of equation, how temperature can be updated in time by implicit Euler? if there any other numerical technique I should use to solve my ODE.