Numerical Method for Three-Dimensional Heat Conduction in Cylindrical and Spherical Coordinates