I am trying to numericaly solve the heat equation in Mathematica.
One of the boundary conditions:
Lim (dTe/dr) = 0 when r -> Infinity
How can I handle this condition in Mathematica using NDSolve????
Attached is a file with what I have got so far\
Attachment: EarthTempMetric.nb, URL: ,