% 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(0,[T01;T02])

% integração numérica
help ode45

[t,T] = ode45(@gradienteTermicoDois, t, [T01,T02]);

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