13 ......... Nx:    número de volumes de controle (direção x)
    13 ......... Ny:    número de volumes de controle (direção y)
     1.0d+0 .... Lx:    comprimento total do domínio, direção x [m]
     1.0d+0 .... Ly:    comprimento total do domínio, direção y [m]
     1.0d+0 .... k:     condutividade térmica [W/m.K]
     0.0d+0 .... To:    estimativa inicial para temperaturas
   500 ......... itmax: número máximo de iterações (método de Gauss-Seidel)
    -1.0d-12 ... tol:   tolerância