Abstract
Errors of numerical solutions for unsteady heat conduction problems by two different methods were evalued, compareing with analytical solutions, two methods were implicit and Crank-Nicolson method. Differences of numerical solutions and analytical solutions were very small and increased with dimensions of problems in three dimensional problems, errors of numerial solution by Crank-Nicolson method were about 0.7% in central portion of analytical region.