I am trying to solve the following PDE:
\begin{equation} \partial_t u = \partial_x^2u + u - f(x,t) \end{equation}
with the initial condition: $u(x,t=0) = 0$. ($f(x,t)$ is known).
The first thing to think about is to separate the $u$ into $u = u_h + u_{nh}$, where $u_h$ solves the homogeneous version,and $u_{nh}$ solves the nonhomogeneous one. However, this does not work, as we will return to the initial problem for the latter.
So, what is the method to solve this kind of nonhomogeneous PDE?