The analysis of effectiveness of the gradient algorithm for the two-dimension steady state heat transfer problems is being performed. The three gradient algorithms - the BCG (biconjugate gradient algorithm), the BICGSTAB (biconjugate gradient stabilized algorithm), and the CGS (conjugate gradient squared algorithm) are implemented in a computer code. Because the first type boundary conditions are imposed, it is possible to compare the results with the analytical solution. Computations are carried out for different numerical grid densities. Therefore it is possible to investigate how the grid density influences the efficiency of the gradient algorithms. The total computational time, residual drop and the iteration time for the gradient algorithms are additionally compared with the performance of the SOR (successive over-relaxation) method.