So I have this PDE: d2T/dr2 + 1/r dT/dr + d2T/dθ2 = 0. How do I implement dT/dr || [r = 0] = 0? Also, what should I do about 1/r? This is actually the first time I am going to attack FDF in polar/cylindrical coordinates. I can finite-difference the base equation fairly decently; I am just having a hard time in implementing the derivative boundary condition at r = 0. Can someone give me an idea what to do?