function Tdot = gradienteTermico(t,T)
    k = 400;
    A = 8000e-6;
    L = 0.5;
    m1 = 200;
    c = 466;
    T2 = 300;
    Tdot = -k*A/(m1*c*L)*( T - T2)
endfunction