hules.mws

>    deq1:=c1*m1*diff(T1(t),t)=-C*A*((T1(t)/100)^4-(T2(t)/100)^4); deq2:=c2*m2*diff(T2(t),t)=-C*A*((T2(t)/100)^4-(T1(t)/100)^4);

deq1 := c1*m1*diff(T1(t),t) = -C*A*(1/100000000*T1(t)^4-1/100000000*T2(t)^4)

deq2 := c2*m2*diff(T2(t),t) = -C*A*(1/100000000*T2(t)^4-1/100000000*T1(t)^4)

>    with(DEtools):

>    c1:=1; c2:=1.2; m1:=1; m2:=1; C:=1; A:=1;

c1 := 1

c2 := 1.2

m1 := 1

m2 := 1

C := 1

A := 1

>    DEplot({deq1,deq2},{T1(t),T2(t)},t=0..3,
[[T1(0)=273,T2(0)=373]],scene=[t,T1(t)],stepsize=0.1,linecolor=t); DEplot({deq1,deq2},{T1(t),T2(t)},t=0..3,
[[T1(0)=273,T2(0)=373]],scene=[t,T2(t)],stepsize=0.1,linecolor=t);

[Maple Plot]

[Maple Plot]

>