tags:

views:

59

answers:

1

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.

+1  A: 
Conspicuous Compiler