13 ......... Nx:    número de volumes de controle (direção x)
   13 ......... Ny:    número de volumes de controle (direção y)
     1.0d-2 .... gama:  difusividade térmica [m2/s]
     0.0d+0 .... To:    estimativa inicial para temperaturas
  1000 ......... itmax: número máximo de iterações (método de Gauss-Seidel)
     -1.0d-10 ... tol:   tolerância
     1.0d+0 ..... beta:  fator de mistura (funções de interpolação): 0 = UDS; 1 = CDS