function ydot = gradienteTermicoDois_Octave(y,t)
  global k A c L m1 m2;
  M = k*A/(c*L) * [-1/m1 1/m1; 1/m2 -1/m2];
  
  ydot = M * y;
endfunction