9 .......... N:     número de volumes de controle
     10 .......... M:     número total de passos de tempo
    20.0d+0 ..... tf:    instante final para solução de T(x) [s]
     0.5d+0 ..... theta: formulação empregada (0.0 = explícita; 0.5 = Crank-Nicolson; 1.0 = implícita)
     1.17d-4 .... alpha: difusividade térmica [m2/s]
     0.1d+0 ..... L:     comprimento total do domínio [m]
     0.0d+0 ..... Ta:    condição de contorno na face oeste (x=0) [ºC]
     0.0d+0 ..... Tb:    condição de contorno na face leste (x=L) [ºC]