% define parâmetros como globais para acessar de dentro da função
global A L c k m1 m2
m1 = 200;
m2 = 200;
c = 466;
k = 400;
A = 8000e-6;
L = 0.5;
T01 = 800;
T02 = 300;

% instantes de tempo
t = [0:100:100000];

% testa função
gradienteTermicoDois_Octave([T01;T02],0)

% integração numérica por Octave

lsode_options('integration method','adams');
T = lsode("gradienteTermicoDois_Octave",[T01;T02],t);

% registra gráfico
plot(t,T)
grid
xlabel('tempo - t / s')
ylabel('Temperatura - T / K')
shg