> | 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); |
> | with(DEtools): |
> | 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); |
> |