function ydot = gradienteTermicoDois(t,y)
    k = 400;
    A = 8000e-6;
    L = 0.5;
    m1 = 200;
    m2 = 200;
    c = 466;
    AA = k*A/(c*L)*[-1/m1 +1/m1; 1/m2 -1/m2];
    ydot = AA * y;
endfunction