Hi Adokiye,
I’ll try to write something for you, but it might be a week or two before I have time. Feel free to e-mail me the details of your problem. From your description, it sounds like it will be best to use polar coordinates. In the meantime, here are two things to consider:
1) In general, parabolic PDEs like the diffusion equation are easier to work with than elliptic PDEs. Mathematica’s built-in routines in NDSolve can probably do a lot of the work for you, saving you the time of setting up a finite difference method.
2) If you do want to get a start right away, Igor Mandric and Ecanterina Bunduchi have a demonstration project here: . You can download the author code and take a look at what they’ve done.

Hi, Ben I appreciate your response. If the NDSolve function can solve the problem then i will be very happy with that than the FDM, here, i need to take boundary
conditions that will create hole in the middle of the plot.
Please assist me.
cheers,