Alguém tem algum código ou algoritmo de cálculo da equação diferencial de calor 1D??

Estou fazendo um trabalho que é pra fazer um programa em SCILAB, usando o método de diferenças finitas. Problema: Suponha uma barra de metal isolada termicamente, com as suas duas extremidades em contato com blocos de gelo a 0º C, e aquecida instantaneamente em seu ponto médio por um maçarico. Qual a temperatura da barra após um certo tempo? Lembrado que sua temperatura até atingir um instante t deverá cair até 0º C.