Determine the closed-form solution for the non-homogeneous recurrence relation $a_n = 3a_{n-1} - 2a_{n-2} + 2^n$ with initial conditions $a_0 = 1$ and $a_1 = 3$.