I have the following PDE

such that A is a constant

I've set up my ODEs using seperation of variables to get

where k_1 and k_2 are constants.

How do I account for my source term (f(x,y,t))? I'm reading up on eigenfunction expansion but so far it's only for the dimensional case.