I would like to numerically find the solution to
u_t - u_xx - u_yy = f
on the square y ∈ [0,1], x ∈ [0,1]
where f=1 if in the unit circle and f=0 otherwise. The boundary conditions are u=0 on all four edges.
I have been trying to use ND-solve for ages but I keep getting error messages. I'm new to Mathematica so I don't know how to define f before ND-Solve.
Thank you in advance.